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Self-excited rotor whirl represents a serious hazard in the j 

operation of turbomachines . The reported investigation has, j 

therefore, the objective to measure the lateral forces acting j 
I on the rotor and to determine the characteristic pressure 1 

distribution in the rotor clearance area. A description is I 

presented of an approach for calculating the leakage flow in 

i the case of an eccentric rotor position on the basis of empirics] 
loss coefficients. The results are reported of an experimental 
investigation with a turbine stage, taking into account a varia- 
tion of the clearance characteristics. 
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In an ©ngin©©ring ©ffos^t t© gjignifieantly rais© tho output of 
turlsordackin© imitiSj, in iiigk pressure tunbinss th©r© a,pp©az“? to an 
©TCr ineneasing ©stents flezuml vibration® unrelated to tSi© 
rotational ®p©©ds plaeing sever© operating restrictions on rotors 
tliat are theriiiodynamically properly designed. Shis output^dependent 
©sssitation of self “generated vibrations is not caused by bearing 
instability s but by forces produced as a consequence of the 
clearance flow generated betv/een rotor and housing She 

d©velox»ent of these forces is generally suriimariaed under the 
concent of clearance es^citation. 


A theoretical description of the vibrational system is already 
available [9 s 10] p for multi-=supported shafts of any form. However, 
the treatment assumes Isnowledge of the sup|)ort oharaeteristies and 
the clearance asscitation forces. Recently, several important papers 
have been published (for instance, [ 6 ]) in the area of friction 
bearing research 5 x^dth their help, bea,ring instability (oil whip) 
can be substantially avoided by means of constructive measures. In 
addition, it is possible to determine system damping, \vhioh in a 
vibrating turbine shaft is predominantly caused by the bearings. In 
publications to date on rotor insta.bility due to clearance 
esscitation, the essciting forces are determined almost exclusively 
via theoretica .1 sta.tem®nts based on Thomas “ [ 1 ] fundamental 
considerations. Only recently have measurements become available 
[53 that allow a, reliable estimation of the limiting output. 


She clearance excitation forces originate In the sealing clearance 
= variable around the perimeter » that ot;cu£*jg whari a deflection 
exists between turbine rotor and cs-sing. Sue t/i the changing 
clearance loss, the rotor blades are subjected t© differing 
peripheral forces, the resultant of ?/hieh has an excitina 


•••■ numbers 


in the margin indicate foreign pa,gination 
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js’ a sliaft iriM'=atiag ia ‘feSie sara© s©e©© as tii© diree'&ioa 


of £>otation« Sine© ia tim'^Ma© stages the 2?©t©5?°s peripheral 
veloeity is var^r lars©s th© flow through tli© soaling olearaaoe 
is of spiral type* For aa ecceatrie rotor position, the 
conseouQBs© is a pressi^r© distrihiition that ¥ari©s along the 
periioeter - due to the different flow crosses© et ions •= whieh in 
the ease of handed "huekets" considerahly magnifies the eseiting 
f orees « 


She goal of this study is th© Bi©a,surea©nt of the transirers© 
foress aeting on th© rotor, as X7©11 as the determination of the 
characteristic pressure distrihution in th© rotor clearance q In 
a.ddition, a procedure is provided with lAich the clearance flow, 
affected by torsional forces, can he calculated for eccentric 
rotor positioning, by means of empirica.1 loss coefficients. She 
result contains the variable clearance throughput and the 
pressure distribution at the sealing clearance under 
consideration. 2his provides the tivo characteristics of the 
clearance flow through v/hich transverse forces proportionsil to 
the lateral deflection act on the rotor. 
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2 . 1 . 6 ©n©z*al definition based on a rotos* model 

^ 'Fv " IS’ • 

She oomplex oonfigumtion of a turbine rotor oan b® represented 



Figure 2 a 1 Rotor aiodel 


by sections of constant cro s s= sect ion » provided the subdivisions 
are sufficiently saiall. She dynaaiics of each section is 
described by a system of differential equations which can be 
solved closed 0 in the linear casso Matrix transfer procedures 
are best suited to the description of such elastomechanic 
problems » They allow a particularly clear expression of the 
quantities describing the rotor elements » The knowledge of the 
external forces acting through the flo?/ processes is of 
essentia,! importance «. They shall be defined below in a generally 
valid mannero 


Figure 2olj above 9 shows a simple vibration models suitable for 
fundamental studies. For the spatial coordinates x, and the 
g^assis describes the shaft's stationary rest position. If we 
represent the shaft's deformed center^line in a top and side pro- 
jection, according to Figure 2.2j, belov/s then v/e have, for the 
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Figure 2„2 Pore® definition on rotor aiodel 

As usual 5 the eirbernal forces acting on the rotor are defined as 
positiTe in the direction of the coordinates.. The same applies 
silso for external moments s in the right-handed system choseno 

Assuming small dynamic displacements c the external forces and 
moments occurring in addition to the stationary load can he 
assumed to he linear e, In terms of this linear theorys the load 
vector as a function of the motion vector •!!> can he 

described as follows b in matrix notations 

cS ^ • ( 3.1 ) 

If we now consider only the forces F and moments M acting from 
the outside on the turbine shaft, then the load vector ^ and 
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•feho motion vaotor ^ as*© defined fey 


’R. ■ 



Pv 
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tf 


? 

‘f 

Mvl 

t • c] 


D 

'pt/ 1 

^ « 


( 2.2 3 


G 0 n@ides*ing th© @5yram©ts*y ooMitions ° aiaplioafel© "bQisaus© of /J_ 
the ps*©^!©®*® isots^opy- th© forces and moments caused fey tfee 
flo?/ V7ill in g©n©s*al fee deserifeed fey th© coefficients of the 
deflection' amtriss and the velocity Timtriss oU » 
coefficients of this mat2*i3s have the diaiensions of a spring or 
z*espectivQlyj daaiping constant o For a shaft rotating clockv/isej 
their sign uill fee given fey the coordinate system to fee 
introduced later® which rotates with the viferational motion. 
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Q/J = 0 (,, Qg 

“bg b,j • C 2.3 ) 
b^_ cbt;. feg . 


The hearings” dynamic restoring forces can fee represented in the 
spatial coordinate system x^y in a similar manner. In generals 
the symmetry conditions expressed afeov© are not satisfied here. 
On the other hands assuming point=feearingSs the external moments 
vanish. In additions forces due to the hearings” tilt as a rule 
can fee neglected. The deflection matrix and the velocity 
matrix 


for the hearing forces then are 
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gh© l)©ariag©* ehs^aetoristle laS*© eompos©d of fom* 

©pi’iag and damping oonal5an'fe.'3s wIiosj© sign is ©stablisli©:?. "bj tlm 
coordinate system abov©. S'os* tii© usual friction Isearings a large 
volume of tests results is already avaHafel© from Glieniek© [63« 
It should h© pointed out that in those papers the direction of 
rotation is defined inversely j with the eonsequeno© that the 
signs Biust lb© cha,ng©d during the transfer to the coordinate 
system chosen hares 


2s2 SiFijple rotor stability behavior 

In turbine ©ngineeringo the foroe Qgo described by the 
coefficient q^s is of decisive signif icance o It acts 
perpendicularly to the direction of deviation and for a. 
codirect ionally vibrating shaft acts against the damping force a 
This can incite the systera to self “gene rated vibrationo Figure 
2e3, below shov/Sp in a defined coordinate system Zt, y the 
intersection point o displaced by the qua,ntity ^ ivith the angle 
ijij for a shaft rotating with angular velocity Wo The temporal 
derivative of the angle \p represents the circular frequency 
of the vibrational movement o 



Figure 2o3 Yector diagram for circumpolar vibrations 
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Ssi 'feliG 02^0't5Gsa qs'OIMcI %ho p@Giti©j;ia tlio 

e©at2*i£^gal l*@2*@© siow| ae'feo ia doTiatioa di£*©6ljl©si "i” osid 




'fellG 0pE“iSlG fOS^GG -GQ, t©gGtllG5i' witli 2?Gcrl3©2?jing 
fos^eo °0,.j 5 dQtGi?FiilaGd "bj tSiG fl0?/e SSiG damping 'Sq'Sqq S* 
feal©aG©d at tii© limit ©f ©tability W tSi© ^© 2 *©© ^ 2 * against 

tk© viMmtlGsml BGtioa. Ukil© th© £>©0t©.ring ^’©£’0© ©an bo 

a 0 gl 0 et©d{i in tiirbinGG •= boeauGG of it© ©Biall aiagnitwd© •= in 
eordpaj?!©©!! to tli© spring f©re©0i tb© clotGraiination of tbG 
©sseitation fere© Q« is th© main goal of this study. 


Building on Shoams* fimdamental work C1 ]j s©'^©ral papers /£ 
eontain theoretioal eonsiderations on th© ©ff©ot of ©ssternal 
forees on th© st@,bility behavior of sisapl© rotor aiod©l 0 . In 
analogy to Shoaiass Sasoh E?] investigatss = for rigid bearings - 
th© ©ffeet of only th© ©zeitation constant ®» o©ntral 

arrangement of th© aiass, I'/hil© Pilta [133 takss into s,eount all 
th© eosffioients of th© matrioes P-.,. and % in an asymmstrieal 
rotor® Sh© stability limits for th© simply engaged Laval shaft® 
taking th© effect of friction bearings into account through 
eonsideration of th© cl©ar©,ne© ©sseitation foro© Qg® 
ealeulated by ICra©ii©r [2]® Pollmann [33 and ?og©l [93 » 


Ml ©ssampl© of such calculations 
sts,bility value was forntsd v/ith 
© 2 seitation constant v/hich is 


is shoivn in Pigur® 2o4® b©lov7» 
th© shaft rigidity q and th© 
plotted against the shaft's 


(S’ 


% 


/s 


ajustment Besides the bearing type® the relative shaft 


iriQ 't 
c 2frv' 


elasticity ims varied® xThich 
under its own weight® to tha 


indicates th© ratio of static flesture 


friction bearing's diaraeter play. She 



SemsiQS’iQM coef fie Ion'S’ 





Adjustment to/W[^ 


Figure 2 ,^ Stability chart for the symriietrical imit^raass 
vibrator with friction bearings. Soj^ = 0.2§ b/d = 0,5 


and the bearing's width ratio b/d is kept eonstanto Below the 
characteristic limiting curve the shaft is stable ® while above it 
instability begins. The steep drop of the rotational rate - which in- 
creases v/ith ” characterizes instability due the friction 

bearing's so-=called oil whip. 


Essential information regarding the construction of a turbine shaft 
can be derived from stability charts of this kind. However, a 
quantitatively accurate vibration calculation is possible only when 
the rotor's geometry has been defined more accurately, by 
subdivision into individual fields, and once the shafting's 
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filurijlpl© support can "Jj© taken into aeooirnt# Ski© wa© ©olv©d» £ 
iastaaae, hj aatniz t5?an©S©E* pi’ooedm'*©© , ky Ssliirmes? [10] and V 
[9], a© well a© hj Gasek with finite elements [8]« 


)g©l 


2.3 Flow=-g©n©rated foreea 


leseareli 


status 


/I 


When turhin© rotors are mounted in casings, at th© non-contact 
sealing clsaranees of th© stator and th© rotor, of necessity 
clearance widths will differ along the ciremiif erenoe • Shis leads 
to clearance loss and tangential force differences at the rotor 
and in turn, to a resulting transverse force ®,t the shaft. In 
addition, an asymmetrical pressure distribution appears a.t the 
sealing clearances, which generates a force that acts on the 
rotor. For a vibration calculation, both components - the 
transverse fore© from clearance loss Qg and that from the 


pressure distribution Q 


D 


Biust be added I 


Q “ Qs Qd • C 2,S ) 

To the extent that those forces Q are linear and a function of 
the deflection £, the constants q = Q/e of the deflection matrix 
of equation ( 2 o 3 ) are easily determined. 


The flow velocity caused by these forces in the usual turbines 
are two orders of magnitude higher than th© rotor's vibration 


velocity. It is for this reason that forces from the matrix in 
equation (2.3) that are proportional to the velocity, can % 
usually be neglected. 


2 . 3 fl Forces from the variable tangential force at the rotor 
grid 


2.3«l«le Eccentric rotor position 
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A piijsjiaal Qssplanation of an © 22 ©itati©n f©5?©G [i] s'Sas’ts 

til© pz*©mi©© tlmt a tii 2 ?'binQ z^otor as siiown in Figus?© 2,5, 
■beloiv, will ©sspe 2 ?i©na© imevoa tangential forees, fos* a do* flection 
©. For a asiall loeal sealing eleamnoe (i!^ ^ 0), fiu© t© tii© low 
eleamne© loss tii© p©ripii© 2 *al foree will "b© l©,rg©n, whil© it 
vdll 1m eoi^respondingly smaller at th© diametrioally opposed 
aid© . 



Figure 2,5 Derivation of the clearance excitation forces 


The peripheral or tangential force U. of a turbine stage is 
obtained from the specific work, a. , the throughput, m, and the 
tangential velocity, us 


U; 


•O; fil 


C 2.0 ) 


The dependence of the internal tangential force on the 
sealing clearance width can be represented either by the minimum 
work, Cl4]u or as in [1], by a quantity-loss, m_„. Both 
effects can be combined in the tangential velocity lost due to 
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C 2.7 ) 


tiiG «3i©amae©, 


11“ y 




U‘, 


5?e U is til© taagesitial fora© witlioiit olearane© losses. 
y. 

define tlm isentropio tangential for go, U_, that would "b© 
attained v/itli a loss^fne© flow, then th© internal and the 
peripheral effioieney, as well as the olearanc© loss can be 
represented as fore© ratios 


IJ, 

IT^' 


Us” 


r : . 

U. 


C 2.S 


W© 


/H 


For th© local tangential foro© v;e thereby obtain the simple 
©qua, t ion 


ciU-, 


C3 



C 2.9 ) 


The integral of this force along the rotor perimeter yields “ 
using the coordinates of Figure 2,5 ■“ the forces acting on the 

rotors the constant tangential efficiency, n , cancels out. 

1/1 



2-^ 

( r. 

U_s 

7,^ 

1 

=> ” 1 -“Jl; S'fi'-P “ 

2tr 

J Csp ? 

* 

V 


} i 

Zxr 

jr 


=- J fUil cos If 

Us 

V icp ® 

i 

0 


0 



( 2.10 ) 


lEhis integral can be solved only if the local clearance loss 
along the periaieter of the variable sealing clearance is Imovm. 

There are equations for the dependence of clearance loss on 
clearance width in several of the papers comparatively reviewed 
by Winter [11]. According to Traupel [lA], clearance loss can be 
expressed as the ratio of lost work, a , caused by the sealing 
clearance to the available isentropie heat gradient, Ah 
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( 2.11 3 


If ‘fell© ieentropio tangen'feial fore© is defined 


IL 


m£A 


( 2.12 ) 


/M 


then the clearance loss, t, , 

sp 


is identical 


where the tangential force, reduced by the 
takes the va,lu© of equation (2.6). 


1 0 equation (2.8), 
clearance effect 


She stage clearance loss is composed of a loss at the stator 
blades, , and a loss at the rotor blades, 

Sp bp 


Up 




^S? 


a; j<: , -Cp x" 

T ^ ' A"' 'fi ' 


It depends primarily on the ratio of the clearance area, A , to 
the grid transverse area. A, and on a factor K, whicn essentially 
depends on the design and the inclination of the seal. 


I 

I 

i<; K" 



6a , flp 


Figure 2.6 Clearance loss coefficient 
("buckets" without shroud band) 
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"'biielsQte'' ivithoiit eiiroud tends aeeoMing to Figure 
IS' OS's respect i¥ely, IC", depend only on tlie angle 


For "blades 
E»6, abore, 

Ms or respectively s A3s on f loi7 deflection. For buckets with 
shroud bands and stator bases a dependency results on the 


already mentioned inclination of the stator, 2Ah' or the 
rotor, respectively 5 2Ah" /w|,-, according to Figure 2.?, belov/. 
Here 



Figure 2»7 Clearance loss coefficient 
("buckets" with shroud band) 


we must appropriately use the inclination a.t the base of the 
blade for the sta.tor, and for the rotor that at the blade head. 
For labyrinth seals, the clearance loss is reduced depending on 
the niffiiber of clearance peaks, z*, or respectively, z". For the 
areas, in the general case of a ehaaiber step according to 
Figure 2. 5? we shall have 

A* = ir I sin 04^ 1 A = ir dm ^ S.nAa J 
A'sp - s' 1 ° ' 
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W© tilus obtain, Sot tli© cl©as*ano© loss©, 


\[^ s'tnot^ 


4 


J 1 

■;:.i/3. 


V? 


ik 

dm 


1 “ 


( 2.J3 ) 


Consequently, the clearance loss at the stator and at the rotor 
depend linearly on the elearanee v/idth. If v/© 
neglecting compensating flov/s " that the dependence is 
applicable a,lso to local elearanee widths 


I' = s' - Q cos <(» , 
l" o s" ~ 6 COC ^ 


■ 

• { a.1'1 1 


according to Figure 2.5, then we obtain, from equation (2.10) 

/U 


Ois “ 
0X5 “ 


0 , 

<1 u 1 f JlL da 1 

I t' I ?? sin«, dm i: 


+ 




( 2.1S ) 


Because of the direct proportionality between clearance loss and 
clearance v/idth, the force vanishes in the direction of 

deflection. 


The clearance excitation force, Q 2 g» <i®P®^ds linearly on the 
eccentricity e/l" mentioned, the isentropic tangential force, 
and the coefficient ^£ 3 * which basically describes the 
construction form of the seal. 


£23 

u; 

K: 


t< 


2.S 


' 2 a 2 


e 

? > 

.J<L. dull 

\ig'sino<4 dm t* 




cJtn 


( 2.16 ) 
( 2.17 ) 


In the dimensionless representation of equation (2.16), 
indicates the slope of the excitation force mentioned, /U , 
over the relative eccentricity e/l" and is therefore referred to 
as clearance excitation coefficient K 2 s’ 
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la paineiplOf aon-linoas* eleamneo loss equations eouM also "be 
usod for th© iatogpatioa of the loeal tangential fore©, such as 
those oontaiaed in th© 2nd edition of [14], for instance* Shis 
effect is disoussed in section 4*4o2, in eomieetion with a 
measured efficiency distribution* It is shorn there that even 
with a non-linear approach for forces can develop that 

depend linearly on th© eccentricity £, just as h©rat 


2.3»l»2e Rotor-to-housing inclination 

Forces acting on th® rotor can also be caused by an inclination 
of the rotor ivith respect to the housing, if the a^tial sealing 
effect of the clearance is significante In general, according to 
Figure 2«5« the inclination of the rotor along the shaft’s 
bending line is coupled to a certain eccentricity. If w© start 
from the premise that the radial clearance v/idth is large in 
comparison to the assial sealing clearance, then the clearance 
loss will be determined only by the axial clearance. For this 
type of construction clearance loss equations exist, as in [14], 
conforming to Figure 2.7 that make it possible to caleiilats the 
forces a,eting on the rotor, which in this case will be due only 
to the inclination, The local axial clearance of the rotor 


with 




( 2. VS ) 


shroud band towards the housing changes along the perimeter 
according to Figure 2.5s for an inclination of the disc equal to th© 
bending angle a® The ms,gnitude of the non=unif©rmity a depends 
on the angle of inclination <x, and on the diai'ieter d^ of th© 
seal. The clearance loss is again taken as linear with respect 
to the ratio of the seal area and the rotor areap 




I" 


( 2.19 ) 
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Ab wasj vtio ©aso f@r tlio mdial seals tliQ elea<?aE©Q loss 

IL for til© assial sfotos* seal depends on the already 


eoQifiaio: 


ge^erred to rotor inclination 2Ah^/Wg^ and can he talsen from 

j?igBrQ 2 ofo 


In analogy to an eqoentrie rotor position, a local tangential 
fore© aan "oe fom©d that depends on the elearano© loss, ? 
and til© isentropic tangential fore© Ug (equation (2.12)) 




Inta/p-at: 


along the perimoter yields the transverse forees 
aeting on the rotor. 


9 




Qss I ^4 f Sm/Sa. dm ‘ 


( 2.20 ) 


/ 16 

In acoordance v/itli its definition in the deflection matrix, the 
for©© has a, destabilizinq effect on the rotor, it depends 
linearly on relative deflection a/l" and on th© isentropie 
tangential fore©, and it ©an he represented hy th© ©o©ffi©ient 
as the slope of the dimensionless esseiting force, over the 
relative deflection. 
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( 2.21 ) 
C 2.22 ) 


She exciting force due to the inclination of th© rotor disc can 
he calculated only for a purely axial seal effect, with th© 
clearance loss equations given. Hov/ever, it vrould h® on th© saf© 
side to consider this effect, in addition to th© clearance 
excitation coefficient for a radial seal. For a clearance that 


shows both a radial and an axial seal effect, in a rotor 
position as shown in Figure 2.5, due to the inclination a 
non-uniform clearance loss will occur along the perimeter o For 


20 



this r©ss©si 


inlet turbine stages x«7ill foe 


sufojQot t© g^’eatos’ 


Qssoitatioa iOX’ees, foeeaiis© ©f tli© 


,3»2 F©5?ees du© to th© pressure distribution at 
elearanoo 


sealing 


Srutnoivslsy [153 provides a eomprehensiv© revieiv of the 
ealeulations for non°eontaet seals. While Most of the procedures 
deserihed there determine only the throughput of a elearanc© 
that is uniform along the entire perimeter, vm are interested 
here primarily in the effeets by meajis of whieh transverse 


forces act 


on the rotor due to ele^.'^ance flowi 


Per an eccentric position of the rotor with respect to the 
housing, not only the throughput but also the pressure gradient 
a.long the direction of flov/ changes with the local clearance 
iiridth. Assuaiing that at the exit the flow ivill be subject to 
atmospheric pressure, a characteristic pressure distribution / 
will develop along the perimeter, v/hose maximurii value will 
coincide with the na.rrov/est clearance. Shis effect was first 


described by Loms-kin [I?] using a sreiooth clearance with purely 
axial flow. A similar derivation can be found in [18]. Because 
of the assumption of a purely axial flov; xvithout compensation at 
the perimeter, the pressure drop is linear- If we apply this to 
the shroud band of a turbine rotor, as in Figure 2.8, belov/, we 



Figure 2.8 Rotor clearance 



©■b'baisi “ fOE" tlio ® 

pE'^Qsgm'’© vaE^la'i^ioa p(s,i!0 


f eoordiaafeop of Hgui’© 3.5 ° tjh© 
ia assial and tangonijial diz^o&'i^ion 


pis,’?)-p^ 



1 fZ ' 

b<\. 


( 2.23 ) 


'SliQ frietion oo©ffiei©nt X pip® hydaaulies is assuaied 
eonsfeant, h©i"®. In addition, ©ntx^ana© losses and a possilsl© 
pE'essur© reGovsE^y at tii® ©zit have been n©gl©etad. 


The integration of a variable pressure vas^iation yields, first, 
the foi*ees acting on the rotor in general, 

t p b 2i? 

Q'lo “ “ J 1 ” ” j 1 P(®'‘P) *" ^ 


0 0 


0 0 
b 2 tt 


Aid " ~J dP » -JJ p(s,<^) s'lHc^ T ej«|> da . 

0 0 0 0 


( 2.24 ) 


Using equa,tion (2.23), w© obtain the forces 


/l8 


( Pi- Pi\ f (t) /TTxTv^ 5 

Qie ” Q * 


The force Q|q» 


because of its negative sign, 


acts against the 


deflection £ and a,t larger pressure differences " ?2 ” they 
occur, for instance, with boiler feed pumps ■= have an essential 


effect on the system's vibration beha^vior. Because of a 
simplified integration of the pressure variation, this force 
depends linearly on the eccentricity, while the dimensions of 
clearance d^ , b and s. occur in sortie non-linear terms. Since the 
pressure variation is symmetrical with respect to the 
deviation,, the force Q^jj, perpendicular to it vanishes. 
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Duo t© tlio s^ta^losiy a tsagoatiaX velocity Is goaemtefi 

at tliG olocrmneOf wliicSi in gas l3©a5?lnes, for instanoo ° 
tli 02 ?©iiglily iavostlgated fey Stingolia [g?] •* geaemtes a ps^ssuro 
BiassimiBi ahead of the aa3?3?o\7©st cloas^aaeoe Hero, too, Qp^ forees 
are developed that could initiate a viferatloml system. However, 
they fee com© important only for sufficiently long olearanees with 
SI very asiall radial elearanee v/idtli» wliieh in general ea.n n© 
longer fee implemented inthefelad© channel of a turfeoaaehin© . 


In turbine stages a large velocity component in a tangential 
direetionoecurs ©specially ahead of the rotor’s seal olearanee, 
with the eons©g,ii©nc© that the clearance flow is no longer axial 
but develops a diagonal flow. With eccentric rotor 

positioning this leads to variable cross-sections along a stream 
tufe@. A problem that can fee adapted to ours was treated fey 
Kiesewetter and Welter [193* who calculated transverse forces 
acting on a conical piston with longitudinal flov/ against it, as 
in Figure 2.9s 





V77m77?777. 

Figure 2»9 Conical piston with longitudinal flov/ 


With the flow as shown, forces act in the deflection direction 
that can cause hydraulic lock in valves, for instance^ 

If Vi?e apply the basic geometrical relations from [193 to a 
cylinder eccentrically placed with respect to the housing, with 
a, diagonal flow ©.gainst it, then the lo?;est flov/ xvill occur 
already before the narrov/est clearance = 0), causing a 
pressure masdmum according to the qimlitative representation in 
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Figure 2.10 CI 3 ), beloxi7 



Figure 

135 


2.10 Qualitative pressure variation 
a 5 Piirely assial flo?/ 

Flow affected by torsional forces 


Integration of this 
yields, besides the 


pressure distribution along the perimeter 
restoring force, also an essciting force. 


Similar 


considerations for the 


flow affected by torsional 


forces 


and for purely a^sial flow can also be transferred to the 
labyrinth seals regularly used in turbines, if on© assumes that 


the pressure drop caused by the seal peaks depend on the local 
radial clearance width (ef. also section 3»3«1)» 

ICostyuk [21], ho\70ver, starts from the premise that at a sesil 
peak the kinetic energy is completely turbulent. Without the 
effect of an arriving torsional flow, under this assuiiiption 


there are no transverse forces to act on the rotor. 


ij 


it 


IS 


displaced v/ith regard to the housing in a direction parallel to 
the a3sis« If the shaft is inclined with respect to the housing, 
vibratioii'-causing forces are generated, due to the non-uniform 
clearance width in the flov/ direction. Assuming a spiral flow 
pattern through the clearance, which could be caused either by 
shaft rotation or an incoming rotational flow, Rosenberg [28] 
studied a labyrinth seal ivith two peaks. 


He obtained the pressure distribution (b) in Figure 2.10, vjhich 
agreed well with his own measurements. But since here the flow 
line pattern is fixed, only qualitative statements 
may be made about the effect of torsional flow. 


24 



sE’dlaf:? ‘5© /il: 


fos* Q, rjoal ia vM%oli Ijlio 

ont^cskQO Qloas’asioo width is Gsallos? than tSiat at tho Qsit» 
tmiia^Qs^so f©2»©0Q dovolop duo to tho displaooaont offoet of tho 
xfiM'^atiag j?@tei 2?8 whieh ©an onoit© th© S 2 f 0 t©si to SGlf-{|’©nGmt©d 
vihmtioa. tJhilG AlfoM admits a© eompeasatiag flov/ ia the 
lah^iath ohaa'bQ^Ss uMo^ these eoaditions» Spu 3 g*k aM Keipea? 
l21] oljfeaisi til© opposit© result "by talsiag this ©ffeat into 
eerisiderat ion « 


Heehreuther [293 p©rforai©4 a fundaaiental study on elearane© 
flow, v/hieh starts hy solving the desorihing differential 
©ouation at a plsdn el©aranG©i using a diff©r©na© aisthod. Under 
th© ©ffeets of rotations boundary conditions can b© satisfied 
for both laminar and turbulent flow,, according to Lomaisin. She 
proeedur© ims then widened to include labyrinth olea,ranees; 
here 5 complete turbulence of th© kinetic energy v/as assimed at 
th© seal 
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Processes at th© meridian channel and the clearance 
entrance 


If the clearance vjidth is variable along th© perimeter, the 
local throughput will vary, as will the tangential efficiency at 
th© rotor blades of a turbine stage,. Shis was taken into 
consideration in equation (2.7) » in section 2e3*l* ^or constant 
pressure blaming without shroud band, Pilts [12] investigated 
both effects, by determining pressure changes in th© bla.des, in 
addition to changes in the local triangle of velocities. Her© » 
a.t least for blading without shroud band •= pressure^-eaused 
forces can act on th© rotor. ¥ith e&\ occasionally severe 
simplification of th© three-dimensional flow, Pilta calculated 
all th© coefficients for the deflection and velocity matrices in 
equation (2.3)® However, he showed by means of vibration /21 
calculations tha.t at least for a central rotor arrangement 
between th© bearings, the additional effects determined by him 


25 



QS?Q saeillf ifi eoiapa2*iGoa to eoasidoriag oaljr th© ©ssoitatioa 
foE>©G Qq aoooMlag to [1]. 


If we etart from th© pi?©miee that despite local velocity 
variations the pere©nta/|© reaction is constant along the 
periphery, even for aneccentrie rotor position, then »> in 
eontra,st to [12] ~ by introdneing esipirical elearanc© loss 
coefficients, all effects that can cause local variations of the 
peripheral force have been taken into consideration. 
Disregarding balancing flows, this is valid ©specially also for 
an efficiency distribution measured as a function of the 
clearance width, ivhieh was used to calculate local clearance 
losses. 


At constant pressure ahead of and behind the turbine step, a 
variation in the percent reaction along the periphery » 
according to section A.A.l. - caused by the local clearance 
loss, is conceivable. Here the turbine stage percent reaction 
is calculated from the pressure gradient and the turbine's 
throughput. The only quantity variable along the periphery to be 
considered here is the stator's cleara,nee throughput according 
to equa,tion (4.13). It turns out that it can very substantially 
affect the result, if one disregards balancing flo?/s ahead of 
the stator. Thus, for a large local stator clearance, a larger 
pressure gradient should be observed at the rotor. This would 
displace the local stator to rotor clearance loss ratio, ivhieh 
could affect the forces from the variable tangential or 
peripheral force to very differing degrees. The pressure 
distribution in the sealing clearances could be affected also. 
Hov/ever, a variable percent reaction in the blading's flov/ 
channel does not necessarily have to cause compressive forces, 
since varying pressure distributions could cancel each other, in 
the ease of banded "bucket "-channels. 


Since the effect of a percent reaction varying along the 
periphery is not yet sufficiently established theoretically, nor 
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doaiMGSTtJGd "bjf aQa!3m'“©a©nts, w© oliall liav© t© a©gl©et it* at 
p5?GS©sit. 2t aaii in addition feo asauraed timt prsasyne difi*QS*enoes 
in tlio raonidian eliamiQl will eanool out raueh raore roadily du© to 
oorapensiating flov/s in tangential din© ct ion, tlian in radial /22 
sealing el©amn©®s, sine© th© nelativ© eccantrioity - referred 
to eliannel height - will h© smaller. Only in ©ztrem© eases, with 
high tangential velocities in the meridian channel, ooapressiv© 
iorees will he ©z©rt©d on long, eyliMrieal rotor parts, such as 


those ohserved during 
standing hlading (i.e 


flov; control in test turbines [53 v/ith 
, large flow^off velocities). 


Sh© influx conditions at th© clearance entrance are of great 
iraportane© to the development of a pressure distribution in a 
radial clearance. La Roche [233 investigated th© effect of a 
point of irregularity in the side-wall of th© bladed flow 


channel, between the stator and the rotor, both theoretically 
and experimentally. To begin with, he showed that the off-set 
height u (see Figure 2,8) is of substantial significance to the 
initial pressure at the clearance and hence also to its 
throughput 0 


The calculations and basic essperimentation were performed by La 
Roche for a bidimensional model in v/hich the flow wa.s 
perpendicular to the off-set, Apparently an optimura off-set 
height could be found, here, for which the losses ivould be 
minimiged. In contrast to the usual construction, it was 
characterized by a negative value for the off-set and a 
ivell-rounded entrance edge, Ho’wever, transfer to an off-set with 
an oblique influx still appears to be somewhat of a problem, 
even though La Roche obtained good agreement with the 
bidiraensional model, with a uniform- pres sure turbine. 


Since the overlapping u of the rotor "bucket" height, as 
compared to the stator blading, varies as the local radial 
clearance width, for an eccentric rotor position. La. Roche's 
equations should be included in the calculations of local 
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olearana© Ioessqs. It timis out, hov/ever, that this loss is v©2*3r 
st3?oagly depeadeat oa @@coada£^ ©fleets (displae©ra©at ©ffeet of 
the rotor "blades* wall boundary layer at the stator* fora of the 
shroud "band edg©), v/liieh raak® a r©lia"bl© applioation iMpossi"bl©« 
Whil© th© off- set oartainly has an ©ff©et, for the a"boT© reasons 
La Roche's results can not b© axaplied to this study. 
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3 CMiOUMSIOI PEOCSDURS FOR SHE SRARSVSRSS FORCES AT 
TUR3IRS GRIDS, CAUSED BY THE CLSARARCS FLOl? 


/II 


By introducing 
to r©pr©s©nt a 


the appropriat© control surfaces, it is possiUl© 
seal clearance of arbitrary gsoBietry as a series 


of 


con’siguous 


stream tubes iirith variable oross” 


•sections. 


If v/e 


consider average velocities, 


at certain reference points in the 


elearane®, then they can be described by the continuity 
equa/tion, as a function of the local cross-section. The 


throughput and the pressure drop along a stream tube can be 


calculated by means of the energy equation, using empirical loss 
coefficients? the flow directions are determined v/ith the aid of 
the theorem of momentum, under special consideration for the 
torsional effect at the entrance. The application of energy and 
momentum equations to such stream tubes thus does not require a 
Imowledge of the processes inside the flow domain under 
consideration. It does, however, assume loss coefficients, to be 
determined from Imown empirical laws. 


Prom this view, the separation of clearance flow from the flov 
along the meridian channel must be possible? this will occur 
only if the blade ends are fitted with shroud bands or similar 
features. Free-standing "buckets"' or blades, coupled to the loss 
of volume, also shoiv losses at the blade ends, caused by the 
energy exchange between the clearance flow and the main flow. 

For this reason, such constructions shall be precluded, here. 

For the same reasons we must establish restrictions for blading 
xvith discontinuous shroud bands, or for stator bottoms with 
balance holes, since the continuity equation of a stream tube 
can not be applied in the manner described, due to a pressure 
equalizing flow. 


In the calculation method below, we study the clearance flow for 
a stationary displacement of the rotor ?dth respect to the 
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houising} whieh would allov; tiie detoraiination of all tli® 
eoeffieient® for th© deviation matnis (2.3)® Hov;©veri "beoause of 
theis* Biuaii greater signifioane© t© tJi© vibration behavior, v® 
shall only deterriiin© her© th© foroes relevant to an eccentrie ^ 
rotor position# Sinee th© ratio of the clearano© flow to th© 
raain flow in th© meridian channel is alr/ays small, vj® shall 
consider th© latter independent of th© rotor's ecc©ntricity, as 
a first approximation. In addition, the calculation siethod shall 
be limited to eases in which the clearance flow aiay be 
considered ineomoressible. 


3 . 1 . Definition of the control spaces at the seal clearance 
The non-contact seals predominantly used for turbines are shovm 


Figure 3«1 Seal clearance at a turbine step 



in Figure 3*1 » above# For reasons of operating safety the 
clearance at the radial entrance E and exit A is usually much 
larger than the radial clea,mee S, which in the newer snaohines 
is fashioned in labyrinth form. 

Sorae constructions show a plain clearance at the rotor, in which 
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ease th© sealing ©ffeot oan also "b© aGeomplishsd by an ©ntrano© 
©dg© (E). ¥© shall eonsid©^*, as a gsnera,! ease, the rotor 
olearano© in Fig-ur© 3,1, v/here either a plain or a labyrinth 
olearanc© oan be used at section So However, the caloulation 


procedure ohoseri will also be applicable to the seal clearance 
at the stator, if the corresponding radii are used. By /2$ 

appropriate modifications to the boundary conditions, it is readily 
possible to eliminate the radial entrance or exit, which allows for 
other areas of application, such as shaft seals of the housing. 


The flow processes at a seal clearance may be considered 
unidimensional, provided the control spaces can be varied according 
to the course of the flov/ lines. These are described by support 
reference points which in the entrance and exit lie on constant 
radii? in the ra,dial clearance they lie on planes perpendicular to 
the axis of rotation. In the tangential direction the support 
points are variable and are determined by the local angle of the 
floW“line tangent to the reference line. The variable 
eross“Sections of a stream tube are then given by the distances 
between neighboring support points on the perimeter and by th© 
local clearance width. In this procedure, the floii? angles are 
obtained by iteration of the basic equations for a stream tube? 
they are assumed Imown, in the following sections. 


3 . 1 . 1 . Location of the support points 


The cleara,nce is divided into stream tubes in th© peripheral 
direction, as shown in Figure 3*2, below, where the selection of 
their reference points is arbitrary, sines the stream tube's xvidth 


is also defined, thereby. While for certain calculation cases it 
may be appropriate to arrange the support points locally closer to 
the perimeter, here we ha.ve esta.blished a uniform subdivision on 


the radius rE of the entrance. Th© support points located on this 
radius have the subscript , being characterized by th© subscript 
k in the periphral direction. With the constant step ?/idth 
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t!:o siappos'j points at the ©ntrance ar© determined by the angle at 

il:o O0nt©r£, 


0» 




starting with these 


support; 


point 


locations th© 


ra'iiaining 


point? 


in th© ilow dirsetion ar© determined using only the looal flov/ 
angles 0 For this reason, together with th© angle at the center, 
in th© flo?/ direction the corresponding distances, Ai|> •- more 
precisely deterrtiined in section J,l,2o » also change. 


She clearance (1, S and A) is divided respectively into nE, nS 
and nA support points, in the flmv direction, as shown in Figure 
3o3» above, characterized by the subscript 4* Using the 
peripheral angle as an exasiple, the complete indexing is shown 
for the radial clearance S. Me have left it out whenever there 
is no possibility of confusion. If we require the support points 
to lie on a flow line, the tangents to the flow lines are given 
together with the corresponding flov/ angles. These can be 
approximated by straight line segments, if appropriate 
assumptions are made regarding the intersection of two 
eonse cut ive t angent s . 


If we subdivide the radial entrance E of radius r„ to r« into nE 

ill D 

support points, as in Figure 3*^* belov/, then at constant step 
v/idth, 


the radii 
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Figure 3°^- Support point location at radial entrance 

'3 + 2AreSign{r^-fj.l j. 


are fished. With Imown flovj^ angles a^, the location of the 
support points is then determined by the peripheral angles 








EW 

• 2 1*** ( k * ^ ) 


9“Sn,k 


As 




El.k *> 


( 3.2 ) 


For the angular composition shovm in Figure 3*^!-s this equation 
is valid only as an approximation, assuming small step width, 
Ar. 


In radial clearances, according to Figure 3*5» below, the same 
quations are used for plain clearances. However, in order to 
Obtain a better comparison with measurements, a variable 
subdivision of the support points was chosen. If the peripheral 
angles remain unchanged as the flow moves around the corners, 
then we have, for all support points k at the beginning of the 
radial clearance 
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For a polats aan bo defiaGfi onl^ isjsiodiatQly 

at the poais» nineo m useful assasaptioas a ‘^oloelty 

distribution eaa bo oado ia tho eSiaiaborso Boeauoe of this faet, 
th© r©latioi::'''.ip betv/soa th© siipport point looation and th© flow 
angles at tho a©q,uire partieular ©igaif ieaae© • Starting 
vdth th© 3 »©>C ^nl k entrane©, it is assuBiQd, 

according to l.ijnr© 3.6, above, that the flow in th© first 
chamber boooEc ' turbulent , following a path at th© perimeter of 
the lab^rrinldi anasibQr that is proportional to th© height h of 
th© ohamber end t© th© floxv angles at th© end of th© rsidial 
entrance . 






Sl,k • 


/30 


According to this, as in th© oas© of the simotli clear®,noe, the 


flow lin© is composed of straight 
of intersection of th© flow lines 


line segiiients, where the point 
can be changed by means of the 


th© weighting factor 0 <g^ <1, 
gg = 0 e 5 0 


which in Figure 3,6 is drav/n for 


fSi,K= 






.4 ) 


In the radial exit, th© peripheral ^gles can be calculated in 
the same manner as in the entrane© 0 If a sraooth radial clearance 
precedes the exit, than as a first approximation there will be 
no change in the peripheral angle at the corner to the exit 

f “ H'^nSjk S 

v;hile in a labyrinth the width of the last chamber must still be 
talcen into ebnsideration: 




Ag-nS 

ft 
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la aimlogj 'bo 


IPOS' tliQ 2’Qs'smiaiag p©3?ipli©aal aagl©o v;© sliall liav©, 

Figm’Q 3 ,^^- 


V%. - T/tHk * • ( s. s ) 


Shae, th© path followed hy the fluid from the olearano© entrance 
to its 032it is deterralned only hy local flow angles and the 
distances to the reference lines# Obviously, eoaiplicated 
clearance forms can also he studied in this manner, provided an 


adecmate relationship can he found between the flow angles at 
the support points and the peripheral angles, lor© precise 
predictions = ©specially in the case of labyrinth seals •=■ can be 
made only I'i/hen it becoaies possible to obtain eaipirieal 
information from flow lines rendered visible. 


3,1.2, Calculation of the control surfaces /31 

In order to determine the local flov; cross-sections , some 
assumptions have to be made regarding the corresponding width of 
the stream tube, rA\]j, If, in a^esment with Figure 3*?j below. 



Figure 3«? Determination of channel width 


we place the lateral limit of the control space between two 
support points, then we shall ha.ve, taking into consideration 
the definition of the peripheral angle as sho?;n in Figure 3,2 
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( S.C ) 


%i,k •* I hil(,ii “ ‘Pi.H-'sl ^ ik‘*2.y..Y^] 


Assuming a rotor displacement ivith regard to the housing 

parallel to the assis, according to Figure 3»1 the local 

clearance widths remain constant in the radial entrance and 

n# 

exit s^, while the radial clearance is dependent on both the 
peripheral angle xp and the eccentricity » £. To reduce the formal 
structures w© desisted at this point from also including an 
inclination of the rotor with respect to the housings Beeause of 


the bending line of the vibrating turbine shaft, the 
inclinations in the domain of the steps ar© small, in any event o 
For this reason we may neglect the effect, at this point. 


The local radial clearance v/idth according to Figure 3.2 can be 
satisfactorily approximated by 

gag tos (p 

where it is assumed that at a support point v|j. that width v/ill 
correspond to the average value of Aip. . . For a sufficiently 
narrovj subdivision at the perimeter, this approisimation has no 
effect on the final result, in comparison to an exa,ct 
integration. Thus we have, for the local surfaces of all stream 
tubes, 

^Ei.k ‘ /'^A.k % ^6 » 

* /"AX.k \ ‘ 
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SiiG eGsitme'Sjion ©©offiGionts u iatK*©clm©oc1. Iiq^q moko it possilsl© 
to talso iiito e@sisiclQ 2 ?ati©n o 2 ®ss«sQGtioa ©oa 0 t 2 *iati©ns that 
©©©US’ diiQ to doacl mtos* ©s? 0ts*oaffi eontmetiona at tli© poalsa. For 
plaia eloamneGS asid turbulent i“l0v/ vrl'Ui u " 1, \70 shall bogia 
b 2 f assuBiiag that th© r©f©r©a©© ©r©ss-s©eti©a is oosaplotely 
fillod b^ th© streams for labyrinth soalSi W©umami«*s [24] 
measured values ©an be used. At eoastant seal peals v/idth A, the 
©ffeativ© fl©i7 or©sa=s©©tioa d©er©as©s with iaoreasiag radial 
elearano© v/idth as Bliovm. in Figur© J»B, 



Figure 3' 


Relative * ridge thickness s/A 
Contraction eoeffioients (from [24]) 


According to itj 


taking th© stream contraction into account, 


a 


clearance that is variable along th© periraeter becoiaes somewhat 
more uniform. It rema,ins to be established, however, whether 
this oontrsction coefficient is irrestrietedly applicable also 
to a seal peak with a diagonal flov; through it. 


3.2 Basic stream tube equations 


An essential advantage of the decomposition into individual 
stream tubes is tha.t the basic equations of fluid flov/ mechanics 
can be ssitisfied unidimensionally for every control space. 
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Shis iBQfees it possible to applgr loss eooffioionts from th© 
lahsrriath theory » as v/oll as tho laws of tiihe hyfiraulias. In 
prinolple, one eouM also ea,lonlate eompressihly along the 
gtreaai tube, using proeediares already ImQvm [153» However, th© 
eont@zt of this study assumes ineompressihl© fluids, v/hieh seems 
particularly permissible , considering the relatively low 
pressure drop @,nd Ma.eh numbers Ivfe,< 0.5 at the clearances 
eonsidered, espeeially that of the rotor. 


3 . 2 . 1 . Boundary conditions at clearance entry and exit 

Th© pressures and velocities caused in the blade strea.® must be 
Imov/n immediately at the clearance entrance and ©sit. Starting 
from these conditions, the throughput of each stream tube will 
result from th© losses in the clearance . Once the average 
crosS“Section calculations have been performed for a turbine 
step (of. section 4.4.1.), it can be generally assumed that the 
flow \vill obey the potential vortex law. Through it v/e obtain, 
from a tangential component e, _ at velocity e at a radius r , 
the tangential component e at radius rs 

If vjQ further assume that the axial velocity e sin a is constant 
along the radius, then w© have 


C sin w. - Cm sin • 


Let the density be uniform throughout the channel 1 we can then 
apply the energy equation 



+ 


•? 

f 


1 2. 

J 


Ss 

? 


It is thus possible to calculate the velocity 
and the angle of flow a at a radius r. 


£* 



the pressure p 


40 



c 3.8 3 


P Pm 



a ® 






A® a iirst approisimation tixis yield® the preesures and 
veloeitie® "before the elearanoe, at the radii r,^ of the entrance 


"E 


and r^ of the essit. 


tJ© may no\7 assume that at the clearance entrance, according to 
Figure 3 *3 / 21 flow line exists that separates the mass stream 
flov/ing through the stream tube from the raa,in flow. Thereby the 
clearance throughput is no longer determined only by the static 
pressure drop, but also by the energy of the incoming 
velocities, whose impulse determined the flow direction inside 
the clearance 0 By the same token, shearing forces due to mixed 
friction at the separating flov/ line could act on the clearance 
flow, which can be neglected here, hov;ever, in comparison to the 
incoming impulse. The effect of an overlapping in La Roche's 
sense (cf. section 2.3»3») shall be neglected, at this point. 

But it could be described generally by mesons of the loss 
coefficients to be introduced later. 


At the exit, the clearance flow becomes mixed with the meridian 
channel flo ?/5 here, because of the differences in velocity, 
mixed friction forces from the main flow will act on the 
clearance's stream tubes. It is conceivable that during this 
process part of the kinetic energy of either the clearance flow 
or main flow is converted into pressure. However, according to 
the transverse pressure eq,uation of flow mechanics (ef. [16]), 
it can be assumed « for most seal constructions =■ that the 
exiting clearance flov/ v/ill be subject to the static pressure 
behind the clearance. 
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®i 0 pE’oeeeaes! cl©so5?il3©ci afeo"^© for th© elearano© ontmsia© aM /33_ 
© 2 zit are expressed below in energy and impuls© eq.uations. 
this end we must first introduce formal coefficients for the 


pressure loss and ineoaiing impulse, whose magnitude will be 
©sjamined more closely only in section 3«3« 


3.2.2. Cont inui ty equation 

For each support point i^ the continuity equation can be 
satisfied with the throughput of stream tub© k. It riiay be 
assuBied that in the cross-section, perpendicularly to the flov/ 
direction A. ,^sina. ,, there is an average velocity, since in 

X @ X 9 Ik 

equation (3 "7) we already introduced contraction coefficients 
p. , that take into account a variable velocity distribution 

X 9 Ik 

v/ithin the effective flov/ cross- sect ion. Therefore the 
velocities 








( 3.9 ) 


can be calculated for all support points in the clearances E, S 
and A, from the throughput m^^ of the stream tube and at constant 
density . 


3.2.3* Energy equation 

In each stream tube, the energy equations can be set up froai one 
support point to the nesct. According to section 3*2. ?/e must 
distinguish here between three basic types of control spaces, 
which are summarized in Figure 3*9» belov;. The pressure losses 
caused by friction against the channel walls or by velocity 
vortexing can be assumed to be proportional to the kinetic 
energy of the flov/, with a loss coefficient 
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c 3.10 ) 




Soponding on definition in th© energy equation, this 

pr©S3E5ur© loss oan oGour ahead of or behind the refereao© point 
at nhieli the average velocity is 


/M 



Figure 3«9 Control spaces in a stream tube 


She control spaces were selected such that the velocities at 
locations I and II could be determined with the continuity 
equation, from the local surfaces and flow angles. The total 
pressure loss occurring in the control space is composed of tv/o 
portions, proportional to the kinetic energy of the entrance and 
eisit velocity. Hence the energy equation is 

f I wj * I + 

/JZ 
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Sh© loss co®ffioi©nts ? d©psM ©ssentialljr on th© looal 
clearance form and the length of the stream tube under 


consideration • fhey ar© discussed in more detail in section 3«3» 

For a plain radial elearanee, according to Figure 3»5s t?jo 
neighboring control spaces vjill have loss coefficients) of th© 
same raagnltud© ahead and behind a support point, because the 
flo?/ paths are of equal length due to the composition of the 
flow lineo Hence we have 


Sai " §u+i j (*• J, 


( 3.12 


This equation can also be used as ^ approximation for th© 
clearance in the radial entrance and exit, as in Figure J,k, In 
a labyrinth, th© loss coefficient in energy equation (3«11) 
describes the portion of kinetic energy that is conserved due to 
incomplete vortexing. describes the entrance loss caused by 
the flow towards the paako Hov/ever, since losses are smaller 
during flow acceleration than they are for flov/ retarda,tion, in 
a labyrinth we can generally set = Oo 


With a diagonal or transverse flow aroimd a corner, as in Figure 
3 . 8 , the determination of the loss coefficients is particularly 
problematic. Here we can define coefficients and for the 
pressure losses, caused by the normal, or respectively, the 
tangential velocity components. The energy equation then becomes 


S t 1 Wj f sin\ (l“Cf) eos^otjd" ] 

^ 2. * l' ( 3.13 ) 

. + 1 Wj [ ^ ^ • 

The above description makes a simpler estimate of the pressure 
loss at a corner possible j for v/e obtain the same 

relationship of equation in v/hieh the entire pressure 

loss in the control room is divided into inflow and outflow 
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'fell© g©a©ml ©nergy ©qiia'feionss (3, 
ci©^eE"il 3 © 'fell© pr©s®iir© d^op (e©© Safel© 




and 3»i3) w© ean 
1* b©lov7) along the 


SAILI 3.1 


Energy equations for the stream tub© 
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QntijfQ soal ©loamnaog ia aeeoadanae vrloli Pigias*© 3.3» 
Sliaiiaa'^iag oil© "^reloeiti©® riieana ©f tb© e©atiau:l'&sr ©cpafeioa 
(3»9)» obtaia ■Sli© pressui*© differeaeee b©twe©a aeigiiboaiag 
suppos“‘6 poia'bQ as a function of the staeara tube's throughput, 
th© loaal flot? oao 0 s°s©etions aad angles, and th© ©iapirieally 
d©tersiin©d pr©ssur© loss ©oeffieients. For th© angles w© 
introduoed th© transformation 




vm shall need later. She relationships of the loss coefficients 
for smooth clearances in accordance v/ith equation (3«12) is also 
talcen into consideration. Introducing the control rtiagnitud© 

/ 

fj, • 1 Plain clearance 

® Labyrinth (?n=0) - 


( 3.14 ) 


the equations given can be applied also to a labyrinth v/ith a 
loss-free afflu 2 s to the seal peak. In general the loss 
coefficients contain th© subscript of the support point under 
consideration; a complete compilation is presented in Table 3.5« 


By addition of all the energy equations in Table 3 . 1 , we obtain 
a single equation that allows us to determine the throughput of 
a stream tube as a fimction of the total pressure differential 
available, provided the areas, angles and loss coefficients are 
known. Since the kinetic energy (velocity energy) of the 
tangential components before and after the clearance have no 
effect on the throughput (cf. section 3*2.6.), we define the 
total pressure loss without that portion: 

* Pe- Pa * 1 4 (1 " $ce) - | 4 ( 1 “ sS ) . 

( 3. IS ) 

The throughput equation then is of the form 
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^Pa - “ 1 4 + 1 4 




( 3.10 ) 


Hoi70T©r, sine© tli© angles ar© also a. f motion of the throughput 
» as ^7e shall see in th© follov;ing section » som© speoial 
considerations (se© section 3*2.5») Biust still h© applied to 
these calculations® Accordingly, Table 3ol provides the entire 
pressure distribution through all support points i,k in the 
seal clearance. 


3.2.4. Momentum equation 


The change in a roonentum entering and leaving a control surface 
is equal to the subi of all ©ssternal forces acting on the control 
space. If the itonentum change is to be calculated in this manner, 
the external forces have to be sufficiently known. On a control 
space as shown in Figure 3*9 there act ~ disregarding 
gravitational and centrifugal forces = compressive forces at the 
limits as well as supporting forces due to friction against 
channel walls. Applying the mcsnentum theorem perpendicularly to 
the entrance and exit surfaces is unfavorable for the control 


space selected, since it requires knov/ledge of the pressures 
acting on the end planes A^ and Aj^, froBi the energy equation. 

If the momentum theorem is supplied in the tangential direction, 
then these compressive forces are eliminated. Then only the much 
smaller forces acting on the lateral limits of the stream tube 
must be determined by iteration from a pressure distribution 
along the perimeter. 


If we take into consideration variable radii a.t the entrance and 
exit of a control space, as in Figure 3»9» then the momentum 



oauili'bjfiuiii for tii© for©©^ aatiiig at tXi© support point is 

m [ fj V/j C0S64.J •= pj Wtj COSOCj j + ' J ^ 


/M. 


3 £ Cos OCj “• C'OSOCjjj 1 


f 3.17 5 


She first term contains th© changes in torsional forces from 
entrance I to ©zit II, iTher© the differing radii have an effect 
only in clearances with a radial flow direction. She nesrc term 
refer® to the compressive and supporting forces effective at the 
control space surfaces, which due to v/all friction or kinetic 
energy vortezing act against the flow direction, in a first 
approximation. For th© nraientum equation above it must be assumed 
that no moments are exerted on the control space due to the 
forces acting on it, which is plausible if we take into 
consideration the impulse forces normal to th© peripheral 
direction. 


Prom considerations based on a imple control sps,ce, it can be 
reasoned tha.t frictional forces as well as those that balance 
the flo?; during vortexing, must be proportional to the 
throughput and the velocity defined as reference velocity in the 
energy equa.tion. !I!hus, by means of coefficients f we can in 
general encompass the support forces 

S “ tn I ( 3.18 ) 


that are active due to wall friction. The magnitude of these 
impulse loss coefficients f is determined in section 3»3 ij^ 
connection with the pressure loss coefficients q. 

The compressive forces P acting on the control space are 
determined from clearance geometry and the pressure 
distribution. Hence, the calculation of the compressive forces 


48 



is possiM© only wh&n th© ps^ossur© distplMtion ia all sts?©aa 

ha!3 15© on d©t©5?mined. ®h© variation of t!ie pressur© along 
tSi© p©risi©t©r ean l5© repree©nt©d at loeations i_ by aieans of th© 
S’oiJrier s©ri©[3 


/M 


- Qq 4 a^ODstf -S' Oj,cos 2 (p + ... 4 a^eosm<f 

u > I, • L ' t 5 

4 bj»a»tp + hj^sinZ^ 4,.. + D,„ Sinmsf . 

Frora 2 Isnown pressures p. , in a support plane 2 
corresponding peripheral angles a massimura of m - (j/2)-l 

coefficients can be deteraiined for the function p= ( » for 
instance by a least squares fit method. However, since due to 
their itera,tive determination (section 3*5*) the pressures are 
not given ejcactly, a curve-fitting with fewer coefficients can 
lead to bettor convergence of the calculation procedure, 
especially for labyrinths. 

According to Figure J.lOg the resulting tangential compressive 



Figure 3*10 Compressive forces acting on the control space 
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force P is eoBiposed of the forces ©M P, 


’rom "botli 


sides on the free portions of the stream tube. So this must he 

added the supporting fore© S^, ©sorted on tli© fluid as a 
consequence of the channel curirature. Por a siBiple calculation 
of these forces it is assumed that the pressure distribution and 
the variable clearance x'yidth of a support plane i are constant 
over the v/idth As. . Por a smooth clearance 5 it is possible to 

ci» 

attain an exact solution i7ith sufficiently narrov/ spacing of the 
support points in the flov/ direction; for a labyrinth, instead, 
because of the unlmown pressure distribution inside the chamber, 
no more precise considerations are possible. With a variable 
clearance v/idth 


I’m ■ t + s “ e cosip 


which takes into consideration the height of the labyrinth 
chamber, the individual forces can be determined across the 
width Az. 


1^ - Az plffbl , 

As J ^ p((^) d(f 

fo, ^ 


Ag ej ci<f 

*fa 


( 3.20 ) 


/^ 2 . 

By introducing the coefficients a and b of the pressure 
distribution (3«19)» these forces can be calculated within the 
prescribed interval to in the control space. The 
determined integrals occurring for the support force are of 
the form 
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eos X s5n cp 
bg slo j{ tp 5!n (j! dl<p 


» 


Kieir solution is gonorally possible, witli the usually tabulated 
foraiuiae, for any lumber of eoefficients and shall b© assumed 
Isnown, her®. She resulting tangential compressive force acting 
on the control space now is 

744 


® Aa^ { (li +s) [ - p(tfb) ) 

“ e [ eos (fa p^tpo) - C04 (fb p{<pb^ ] ( 3.21 ) 

T e j p(tp) sinfp dtp I 

ia ^ b'* 

It can be shown that if we assume p(il() = constant, the resulting 
compressive force vanishes. In a labyrinth the effect of the 
compressive forces increases with increa.sing chamber height and 
for variable radial clearance widths is dependent on the local 
gradient ds/dij;. Hov;ever, the sum of all the compressive forces 
along the perimeter must be zero. For the radial entrance S and 
escit A the calculation performed above for the compressive 
forces acting on the free portion of the control space is much 
simpler, since the supportive forces S. (3*20) vanish, because 
of the constant clearance v/idth. 

Once all eisternal forces in the impulse moment equation (3«18) 
ha.ve been calculated, that equation can be used to calculate the 
flow angles at a support point II, if that angle is known for 
the immediately preceding support point I, Since at the entrance 
the flov/ angle is given by the main flow in the meridian 
channel, in conformance to Figure 3»3» the angles at all other 
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suppoi'^'B p©in*B® eaii l3© ealeiilatods U^ing eontimiit^r ©^ua'feion 
( 3 » 9 )» ^v'© ©Main tli© giysteai of eq,y;ations of Saol© 3.2 foi* •feii© 
angulaE’ eliangofs along tJi© atroaui -fenb© k. Sli© irapuls© lose 
eoeffieient® f are indexed Imre in a manner iji'milar to that o 


TABLE 3«2 Stream tub© impulse equations 
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SaMo 3®^ @©ntaia %liQ eomli‘«5i©np (3«3.2) ami (3.1^f*)» iineo 
tliG Gn‘i5j?aii©Q aagl© fiepoMs on vliQ tlisfonglipyt ©f a □ts’Qaa tu‘bQ» 
tliip 2 ?©la'fei©ngliip E’Qaains applioaMo 15© all fm?‘t3li©n angloa. At 
thG olQamnoe tSio a1©¥7 anglG ©f tJiG Xq»s t snp»po5?t loaation 


lnoi“©asop, du© to th© siissefi f^’iotion fos*e© 


wIioGG Riaanitud© 


is fi©toriaia©d in soot ion 3 . 3 . 4* in eonnoetion with th© loss 


eoGffioients 0 


It can b@ inforned dinsotlj from th© ricrasntuin ©quations that /M 
for a r©diietion of tho eross^seotioiml areas A, fro® one support 
point to th© neztj th© flovj a,ngl© must inoreas©, Beoaus© of 
f riot ion or vortesingi f > 0 , of th© velooity energy a 
reduction in th© torsional force is obtained also for constant 


surfaces . 


3 a 2 0 5 e Throughput caleulat ion 

The flow ang 3 .es in Table 3.2 can be represented j for any 
arbitrary support point » as a f motion of th© stream tub© 
throughputs in the form 




\k 



Gg COS 6 ^^ 





( 3.22 ) 


Her© 9 
supper 


th© coefficients b indicate the angular changes due to th© 
t forces g whil© the coefficients c reflect the effect of 


esEtemal compressive forces. Based on the isipulse equations in 


Table 3«2g these coefficients were calculated for all the 
support points in a stream tube, as shown in Table 3.3« 


If equation (3«22) 


is replaced in the energy equation in Table 


3 . 1 * an equation is 
throughput 5 that in 
combination of the 


obtained for th© calculation of the 
comparison to equation ( 3 * 16 ) contains a 
coefficients mentioned above, instead of the 


angle j t 
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S as*G s5irtiila£*l3T 1 

bQ siaipl^ 2?Qp5?eg©at©d iisixig the fosrsiat in 


©aoii ©* 

© 3,k 


Qtm 


MBLE 3,4 Susig to ealeialat© tliroiighput v^itfc equation ( 3 « 23 ) 



By eom'bining the energy and iKxtentm equations, the throughput i 
no?/ gi¥en hy a single equation, readily solved by iteration. It 
is possible in principle , hov/ever, to first calculate the 
angles from Table 3«2 and then the throughput using equation 


( 3 * 3 . 6 ). Of oonrf 3 ©j ia tliat ease 'Gke ealeialations aequiaod 
inci’oase eon®i<l©£*a‘bl3r» depeadiag oa the magaitade of tlie 


torsional fore© 


til w 


til© entrane©. However s tlie underlying 


assuaiption in throughput oaloulations according to the ahov© 
equation is a Irnowledg© of th© local flow eross=sections and 
loss coefficients, which can h© ohtained hy iteration from th© 
course of the flov; lines. 


3.2,6. Discussion of the hasie equations using a sirapl© 
clearance shape 

A dimensionless description of the basic equations for a stream 
tube is possible only if the bomdary conditions ahead of and 
behind the clearance are eliminated, while retaining the effect 
of the seal geometry. Since a,s v/ill be shown later, the loss 
coefficients introduced depend on the local flow angle and 
Reynold's number, a complete similarity can not be attained for 
different initial pressures and velocities. 

The total pressure differential from equation ( 3 ® 15 ) is chosen 
as the reference raagnitude, xvhich besides the static pressure 
difference contains only the a^cial components of the velocities 
ahead of and behind the clearance. To describe the torsional 
forces, the dynamic pressure of the tangential component of the 
corresponding velocity is related to this total pressure 
differential 5 


( 3.24 ) 

( 3.2s ) 


L- 
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2 Cg Xe 

tcti 


c: - 


‘a 2 
f CDS 




While the magnitudes included above are independent of the 
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GlGamnoQ foraij 


s?GfGi»Qn0© sm'*faG© ismst feo defined 


fOi’ 


5T Tg 6. 


{ 3.20 ) 


Fon Q. eeatral sfotor pesitioni let 
ividth at a nadine r-p>« S’roai the rare 
deterriiin© a thnonghput, 


Sg h@ tliQ saiallest elearane© 
satire drop Apg iv@ ©an nov/ 


/so 


”>B 



( 3.27 ) 


v/hieii oorres 


e to a aiassiiflura value » if w© start from the 


prerais© that only the velocity energy at the sma.ll©st 
erosS“ sect ion Ag is completely turbulent. We thus obtain a 
reference magnitude to judge the quality of the seal. On the 


other hands with the introduction of the reference area 
mentioned s we can plot th© variable throughput mj^ of the 
stream pips, as a function of the seal clearance's perimeter. To 
obtain a simpl© description of a stream tube's ba,sic equations, 
w© shall neglect here the compressive forces at the surfaces of 
the control space, which appears permissible at least for smooth 
clearances (cf. section 3»6), If we furthermore consider a flow 


not affected by riiixed friction forces at the exit, then in 
equation (3*23) the subis and are eliminated. With the 
reference magnitudes introduced earlier, we now have for the 


throughput 



( 3.28 ) 


Since the sums S^ and S-^ in Table 3«^ sire quite large, we shall 
show the essential relationships v/ith the example of a labyrinth 
?/ith n sealing peaks and without radial entrance or exit. Here 
we need consider only the coefficients and of the 
velocity ahead of the clearance, as v/ell as the pressure loss 


5 ? 



and iapuls© lo 
p©alss. Ife ttos 


E3S ooGfficients 5i @r £“©gp©G‘(5iTQl]/-5 
Imv© the following eoMitiones 




Q.*c5 


f)E 0 , 


at 

iWD 



nA = 0 , 


SwA 
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)wA 
5 eA - 





n 5 


For the cosffioieniB of Sable 3»3 5i©v; have 

- { 1 |ce 1 j . 

y«1 


/il 


Prom this ive obtain the sums of Table 3»^> 


ri 



where for i = 1 the product above talces the value 

0 



If we take the dependence on the peripheral angle into account, 
described by the subscript of the stream tube k, then we ha.ve, 
using equation (3»28), the referred throughput 
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She stmoture of the equation is maintained even if w© take into 
consideration all the loss coefficients defined in the previous 
sections for more general clearance forms t In this oas©, just 
the sum- terms in the numerator and the denominator will become 


larger . 


If for a central rotor position we assuRie equally sized flow 
eross-seotions at all peaks, then the throughput » for complete 
vortexing of the velocity energy - v/ill he proportional to l//n 
in the Imown manner, where n is the nuaiher of seal peaks. /i2 
Assuming equal total pressure differentials, a decrease in 
throughput could occur due to an afflux affected hy torsional 
forces, if the expression in the rounded bracket of equation 
(3»29) becomes less than zero. For this, the relationship 
between the coefficients c and 1’ - which will be discussed in 
more detail in section 3»3«2. - is of primary importance. 

The local flov; angles can be determined, in the exaaiple chosen, 
at an arbitrary support point 1 >x^n in stream tube k, from 
Table 3»2 


° CgCOS«g 1 -^ (l-Jee) H' (1-5,,,') . J ) 


In a similar manner it is possible to calculate looa.1 pressures 
in the clearance from the differences for location x,k in 
Table 3 » 1 


“ Pa ^ 2 f* 


y I 


1 1 da\v 



/ 

i) 

'• ( 3.31 ) 

■ i,“X 

1 nj t 

1/ 

k 

and 

the flov/ 

angles above it is now 


possible to describe the pressure distribution only in terms of 
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loss) eoeffiaients and tlie looal surfaees. As a spooial ©a,s©i and 
as shall h© sliov/n y©t (ef. eo[v.ation {3*k9))t iv® oan assysi© a 
dependena© of th© irmentum loss eoeffioients 1’ on tlie psfessur© 
loss eoeffioients of th© form 

1-lL and (l”fcef - 


In this manner th© throughput would be independent of the 
magnitude of the relative afflux energy C^s and v;e obtain^ from 
equation ( 3 « 29 )» 



( 3.3 2 ) 



/H 

Under these assumptions, the flov/ angles calculated above (3. 30) 
can be described very simply as a function of the face areas A 
of the control spaces and the loss coefficients, for a !*ocation 

X 



( 3.33 ) 


According to this equation, for equal surfaces, with increasing 
torsional effects at the entrance and larger pressure loss 
coefficients, there will be smaller flow angles in the 
clearance Q Along the clearance (x>l), however, with complete 
vortexing (c = l) a purely axia.1 flow will be achieved, with 
etga=0. Pressure at the support points can be rendered in a 
dimensionless manner with equations ( 3«31 ) to (3.33) 
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px,k ” pA 


( 3.34 ) 




^ i r . 



o 


Hes*© the ©ff©ot of torsional influs^ oan no longer b© dirsetly 
r©eognia©(i, sinoe nov/ it is liaiited to changes in th© local 


surfaces (3«7)» which depending on the course of the flov/ lines 
(3»3) sir© a function of the angles (3*33) • It can he showni 
however, tha,t for stream pipes v/ith decreasing cross=seetions 
there is a pressure increase, compared to those in which the 
cross°sections increase in the flow direction. If, for an 


eccentrically positioned rotor, we start from the premise that 
the floviT through the clearance proceeds in the same direction 
anyv/here along the perimeter, then there will be a higher 
pressure ahead of the narrowest clearance than behind it, due to 
the surface ratios. For the rest, equation (3o3‘^0 applies to any 
clearance form, provided the local flov/ cross»sections and loss 
coefficients are Imown. This equation is hence particularly well 
suited to qualitative studies unaffected by a balancing flov/ 

(cf. section 2o3«2.) 


If one assumes that the floi7 through the clearance is purely 
axial, then for equal radial clearance widths, the 
erosS“Sections along a stream tube become constant, and we 
obtain the pressure distribution 



For constant vortexing coefficients there hence will be no 
pressure differences along the perimeter even for eccentric 
rotor positions. However, if the pressure loss coefficient is 
inversely proportional to the clearance width » as in the pla.in 
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el©a.rane0 ° tli©n a p2?©g0US“o iiiassimiiFii will ©@eiu’ at tli© jriai'^^’owQ^t 
el©a2’ano0» as In Lomakin's eas© (ef. section 2e3»2.). Fus^tkes’ 
pressure distribution tendencies in a seal elearanc© ivitii afflms 
affected by torsional forces can be derived from the sarsiple 
caleiilation in section 3.6. 


3.3® Coefficients to describe flov/ processes xvithin the 
clearance 

Table 3«5 i® s- corapilation of all the coefficients introduced to 
describe the clearance floxv. The first coluran lists the pressure 

TABLE 3»5 Characterization of loss coefficients 


Pressure > . 
loss coeff. 

1 

Impulse loss 
coeff. 

Refer .velea. 
geschwindighi. .t 

Location 
Bild 3 3 

ICE' ^CE 

Ice 


Entrance 

(we 

( tffi 

1 

(e X 

fe X 

'^E i 

Radial ' 
entrance 

3 £0 

( ES ( ES 

(es 

nB 

1st bend 

mm 

f.SE 

”s 1 

i 

fs i 

«s'x 

(8 Radial cle; 

(L (Ia 

f SA 

'% nS 

2nd bend 

ils (L 

(as 

«A 1 

(a. i 

fft i 

^'Ai 

^ radial eKit 

^ WA f WA 



Exit it t 

( CA ( CA 

... 

%A 
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10 E 3 S 5 eoQffieionts, vjliil© tli© eeeoncl displayfu th@ raagni tildes 
5 ?esp©nsi'ble for angular* elianges acoo^ding to tli© muentum 
tii© 02 ’©m. f 5 ? 0 B the eors’esponding rofQr&noo veloeities it is 
possible to readily recognise the eignifioane© of the individual 
subseripts, in eon junction ivith the location designations shown 
in Figure 3»3« 


Below vm shall determine the pressure loss eoefficients ? along 
the stream tube v/ith equations frora the literature, while the 
inonienturi' loss coefficients X can be determined only by means of 
analogy considerations. Since the relation betv/een these two 
'magnitudes is very essential to the throughput of the stream 
tube, an additional condition at the clearance entrance must be 
satisfied, which establishes the rela,tion between the incoming 
velocity energy and the entrance incmentuiri (of. section 3 . 3 * 2 . 3 .) 


3 . 3*1 Pressure loss coefficients 
3 . 3 « 3- • 1 « Plain clearance 


The loss coefficient already introduced for the energy 
equation (cf. equation 3 « 10 ), can be deteraiined » in the case of 
a. plain clearance - according to the laws of hydraulics s from 


the length I, 


the htdraulie diameter and the tube’s 


friction coefficient X. According to [16], the hydraulic 


diameter 


dhwdr - jv- - 2s 

' Ub«n 

is defined as the ratio of four times the flov; cross-section to 
the isretted perimeter, and corresponds here to twice the 
clearance width. The 'tube's friction coefficient depends on the 
tube wall roughness and on Reynold's number 

Re - 
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Fos* kydmiilioally ©mootii valuos as tlios© ia 

Pi^?£’© (©f. Sruelsenteodt [16]) eaa "b© us©cle At tiio seal 



Figiir© 3»11 Tu'b© friction cosffieients for 
hydraulically sBiooth walls 


clearances in 
restrictions 5 
flow. Thus we 


turhin® stators, we shall assme here, without 
the calculation procedure for turbulent /57 
have, for the friction coefficient in von Karman's 


form 


f4 - 2,0 

\fX ’ 


tg (RelfK) 


1 , 0 . 


( 3»36 ) 


Sine© the flow lines in th© clearance can be represented 
sectionally by straight lines, the flov? paths 1 can be 
represented (see Figure 3.12, below) as a function of the flov/ 
an.gle and the distance a to the reference line. We then 
obtain, for the loss coefficients in the case of pla.in 
clearances 


Entrance s 






i 


a * &r 


E 
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HacliSll ©l©Q,'L'’ail8© 8 «} ffl 5 <» e } .a » &Zj 

Sisi’i; 8 S’ « Sa 5 a • 



Figure 3A2 For the caleulation of loss coefficients 

in a plain clearance 


/5§. 

In the considerations above it was assumed that both channel 
v/allB were at resto In order to take into consideration the 
modified pressure loss for rotating channel walls, we start from 
the premise « as we did before •=■ that the velocities w and their 
directions a represent average values across the clearance. 

^Dhus, the relative motion will then affect only the direction of 
the stream tube and the pressure drop. The resistance forces ■=■ 
v/hich in a friction bearing, for instance, '^ause the 
characteristic pressure build-up ahead of the narrowest 
clearance = are disregarded here, because of ' ’le considerations 
in section 

According to Figure 3»13» below, at the wall rotating with a 
tangential velocity u, the relative velocity v will exist, as 
opposed to the absolute velocity at a fixed wall, w. If we 
consider only a wetted surface of the clearance, then the 
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hydraulic diameter 'becomes = 4 s. Accordingly ^ the 

friction coefficient moving along the fixed wall from I to II 
will 'be only half that of equation (3e36). At the rotating 
surfaces the fluid traverses the path froai I to II' at relative 
velocity v and suffers the pressure drop 

§ XI 

“ 2 ^ 4 1 5ln/i> ■ 

in the processs as a first approximation o 

If lire assume an average friction coefficient X for the pressure 
loss of "both clearance surfaces, then the total pressure loss = 
using the definition in equation (3»10) •= will he /51 


A , sing. V \ 

' sl«f2 ' 


C 5.37 ) 


The coefficient C to be used for the pressure loss at a fixed 
channel wall is that of equation (3«3^)» 5?he local relative 
velocities v and their angles 3 are easily calculated, from the 
geometry in Figure 3 •13s as a, function of the absolute velocity w, 
its angle a and the tangential velocity u. At least for flov/s with 
a strong torsional effect, the equaticjx above will provide an 
estimate of the modified pressure loss® 
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3e3ol.2e aleama©Q 


SliG ps’doossQs taking plaoG in a la^jra^inth aliaBljej? ?;itli in©oaipl@t;Q 
V05?t©2sing 0f thQ •FQloeities we2?e invoatigated in detail hj feoddeels 
tZ62» Using an ©sp 2 *©sslon fos* the aiissed f2‘ieti©a 'betweGn tho fluid 
ilov/ing awajr f5?orii tli© seal peak and that at s^Gst in th© ehamUer, 
wiiil© taking into ae count the friction against the fixed channel 
v/alls the impulse eauation ims satisfied vdth the stipulation that 
the pressur© x-amined constant during vortexing. By means of a 


friction coefficient A. 


which conta,ine "both the wall friction and 


the mixed friction “ it is possihl© to calculate the velocity 
reduced along the length ^ in the chaHiber» 


- 


w. 


1 


1 + X 



as a function oi the velocity v/^ at the preceding peak. ^Saking 
this residua.1 velocity into accoimt, the energy equation "becomes 


pi ^ fHt 


pJi^ t % 


With the definition of equation (3«11), w© obtain Groddeek's loss 
coefficient 


/§0 

By applying the iaipulse equation norrually to the flov/, Groddeck 
also o"btains a dependency on the pressure increase in the 
la"byrinth eham"ber, due to the "broadening of the stream upon 
leaving the peak. Shis causes an additional increase in the 
throughput} ivhich ©according to Groddeck, however, can be 
neglected in raost practical eases. She effect of imll rotation 




« 1 


( 77 ^) 

, I . 2s ' 


/ e 'Sio \ 

% ) 
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could l3G liaMled "b.v &i?©cMock's Qqiiaticn as in tliG casG ©f tko 
plain cleamncQo \ 1 q sliall disnogaM it Jiono, !i©woTQr, sine© a 
SQpamtion l)Qtw0Qa tliQ Iossgis ©ai^sod "b^f fnietion and tiios© di^© 
to V05?t©2sinjg is possilal© only with diificultjre 

In an analogous sonsid©mti©n» ffeiisiann [24] took into account = 
in a dioptjfie labyrinth = th© portion of flow ©n©rgy, (which in 
contrast to a tru© Isibyrinth is retained) for which th© velocity 
energy i© oosipletely turbulent. IJhen he applied it to a seal 
ivith a number of peaks 9 he obtained a so-called esseess 
pressure factor, which can be described by raeans of the loss 
coefficient defined in eq.ua.tion (3.10)5 



For a large number of eea.1 peaks g 
©access pressure factor is a linear 


-> a>, according to [24] this 
fun« felon of the ratio between 


clearance width and distance ^ of the peaks 


1 + m 


v/here the proportionality factor m can be determiner, from 
comparative studies between dioptric and true labyrinths. tJe 
then obtain as loss coefficient for a seal peak/ 



This pressure loss coefficient after Neumann is shown in Figure 
3.14, below, in co'mparison to Groddeck's equation (see equation 
(3.3s)). With th© ratio t/s, both equation? . a function only 
of the seal’s geometry, while at least for t hydraulically /6l 
smooth tube there is always the additional influence of the 
Reynolds nuHiber. For the resistance coefficient , Groddeck 
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Figure JA^l- Loss coefficient for a labyrinth seal 
__ after Neuraann [24]} after Groddacls [26] 

gives a value of Ooic Weuraann» using Egli's method (cf. also 
[ 15 ]) finds a proportionality fa.etor of m = 16 . 6 5 from his own 
ealciilations [ 25] 9 it is m = 8.9« At small clearance widths £ or 
respectively s large distances ^ between th© seal peaks 9 the 
limiting value ? = 1 is attained s corresponding to complete 
vortexing. Equations of this kind can be used alsos as 
approxiaiationss for flow affected by torsional forces 9 if 
instead of th© chamber v/idth v/e use the length of the local flow 
line* in accordance with Figure ^ol2. Using a local clearance 
v/idth § at a seal peak, according to [2k] the pressure loss 
coefficient for a labyrinth vdll be 




1 

s-ecosipsi,k 

Asi/sinoc^lf^J 


C S.4G ) 


Hsnes, a more marked vortexing of the velocity is achieved for 
small flov/ a,ngles. 


3.3»1»3« Entrance, bend and exit losses 



tiBlOES t 


AeeoMing to Flguro 3»3» flow ©Mn^ea eoiiE*©© fom'’ 
psrpeMiei^lar-lj t© tli© 5?©fe5?one© ©dgos. Wkile s,t th© ©ntimn©© 
and th© ©sit th© flov/ throuMi the turhin© stators is of 


isiportane© to tli© loss eoeffieionts 
around edges in th© mdial elearane© 
geometry ©an "be of signif ieane© . 


defined her®, at th© flov; 
(h©nd loss) only th© local 


If V/© assOTii© that at th© clearanee ©ntranoe th© ava-ilahl© 
velocity energy is fully used, then th© loss coefficients ivill 


he taken as 




0 (ef. eq.ua t ion (3«50)). In addition, we 


must tak© an ©ntrane® loss into account, relative to th© 


velocity at the first support plan©. Since th© tangential 
component does not change direction in space, w© can set the 
loss coefficient - 0, while the loss coefficient for the 
normal coriiponent is assuaied to he approxiaiately - 0,2, 
although much larger values are also possihle. In principle, 
loss coefficients could occur here that are dependent on, say, 
clearanee throughput or, according to section 2.3*3«» osi local 
overlapping. 


If two plain clearances follov/ one another, then we can start hy 
applying the bend losses from tube hydraulics. It must be 
remembered though, that the flow around such a bend is diagonal 
and additionally, that it can show very variable entrance and 
e3cit surfaces. However, with the energy eq,uation given in 
( 3 - 15 ) » the possibility eicists to consider separately th© losses 
in tangential and in normal direction; here the total pressure 
loss at a bend can be esspressed as a function of influss and 
outflov; velocities. The pressure losses caused by the peripheral 
or tangential component of the velocity can be disregarded as 
= 0 , follov/ing the same considerations as earlier. This means 
that if the control spa.ce's entrance and exit surfaces are 
eq.ual, the peripheral velocity is fully maintained. Considering 
only the normal eoriiponent, we can use conventional bend losses, 
although according to [16] these strongly depend on the 
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eoar/snietioa fom. S'oi’ application to aeal eleanancoai w© eouM 


nso oonstant ¥alii©£ 


.a „ 


0.2 to 1.0, at the pei'^iphe^-’j, where 


the effect of variable aiirface ehameterirtstiee can be deseribeci 
at leapt apprezimatelj by a corresponding subdivision into local 
influx and outflow velocities. 


If a labyrinth seal follo\7S a radial entrance, then the velocity 
vortezing at sunnort location E„r. be expressed in a manner 
similar to that used for a seal peak. If v/e assume the pessure 
loss coefficient® to be the same for the tangential and the 
azial components of that velocity. 



» 


then, in analogy to eq,uation (3*^0), we should have 





1 



( 3.41 ) 


In addition, it would be possible here to provide loss 
coefficients for the entrance loss at the first peak. 
However, with = 0 v/e shall have to a.ssurae that the 


’SS 


pressure loss occurring here is negligible coaipared to the 
vortezing of the flo?/=off velocity at support location An 

analogous procedure is possible if a plain seal (outlet) follov/s 
a labyrinth seal. In this case we describe the vortezing at the 
last peak by 



?/here the coefficient is to be calculated using equation 
(3«^0)« With C^g = C^g = 0, any loss at the entrance to the 
radial ©zit can be neglected, once again. 
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Behind the seal eloamno©, th© flow in th© m©5?idian ehannel 
transfers no velocity energy (Sg_^ ® 0) to the elearanee flow. 

Heno© the statio preesur© of th© main streari will he ir»®osed on 
the exiting flow. Under certain oonditions, whioh depend very 
strongly on th© geometry of the eles.ranc© exit, a portion of /6i 
the exiting flow energy can he transformed into pressure. Shis is 


taken into account hy means of exit losses 


V'J^ 




=*Tfl7A 


< 1, ^7hieh 


very significantly affect the pressure distribution i7ithin the 


clearance. In the normal ease, hov/ever, with 5., 


1, we have to 


assume that the velocity is completely turbulent, at the exit. 


3.3«2. laipul sexless coefficients as a function of pressure loss 
3»3>2.1 Analogy with a straight stream tube 

The coefficients (equation (3.18)) introduced in th® impulse 
equation for the support forces S, which act on th© control space 
due to friction, can be determined “ in a straight stream tube •=■ 
with the length 1, from Figure 3«12. From th© iaipulse equation for 
the direction of flow 

m (Wj- Wjj) + Pj “ Pb ^ ^ 

and the energy equation 

t i 1 1 a 1 a 1 

fpl" Pe) ^ “ |-W2 “ I Wj + P^r 


we obtain, considering a constant flov/ cross-section Aj = Aj-j- = 
A sin a , the frictional force 


s 


“Pv 


I 
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ae a fimotion ©f tli© pj?©^(3Ui?© lo©® 


Q 2 

p v“ 2 


Y 


A different derivation is possible » using the loss-energy 1 
obtained from the decrease in the specific kinetic energy (i/p) 

mul'&iplied by ths mass ixisidG tSi© stream tube under consideration, 




Sine© this loss^ivork amstg 


in turns be ©aual to tho product 


/n 


between the frictional force 
obtain the same result above 
(3.9) and the pressure loss 
the frictional fore© 


and the friction path lengthy xj® 

. Using the continuity equation 
from equation (3oi0)i, w© now obtain 


S O 

■= w w 



( 3.42 ) 


at any arbitrary support point o A compa.rison with the © 2 £pression 
chosen for equation (3.18) provides a simple relationship 
between the impulse loss coefficient ^ and the friction 
coefficient c after equation ( 3 » 36 ) 


5i,lt “ 2 • ( 3.43 ) 

It is assumed here^ ho¥ifever, that both channel imlls are at 
rest. The support force generated by the friction against the 
rotating xvall is determined similarly to the previous 
derivation, from a stream tube v/ith an inclination angle 3„ 


> ,.4. “ A «> I- V*" A sin /^ 

tuft r'»r-rofr 2 . ij. £ s{,n/3> 


C 3.44 ) 
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vjlievB tlm ealei^lated f2?@a ]?igm->© 3 . 1 3 . 

$h© frictional forces against th© direction of flow « doseri'bed 
"bj til© angle a •» ar© deserilaed “by means of equation (3.18), in 
the iBipulsQ equations for th© tangential direction. Sh© aho^© 
frictional force on the rotating pmrtf hov7ei"^®r, according to 
Figure 3*13 acts at the support point at an angle 3. In order to 
deteraiin© an isipulsa loss coefficient which corresponds to 
definition (3. IB) at "both channel walls, we siust hence introduce 
a correction factor. As a f motion of th© pressure loss 
coefficient 5 for Wo walls at rest (equation (3»36)) we then 
obtain 


Jr * I 5 


( 3.4S ) 


For large tangential velocities and sraall torsional effect, the 
iterative consideration of any rotation, using the iaipulse loss 
coefficient above, can fail for numerical reasons. In 
calculations for such eases, the frictional force frorti equation 
(3.44) would have to be introduced directly in the iaipulse 
equation. 

/66 


3. 3. 2. 2. Considerations based on the energy equation 

Especially in the application to flov/ around bends, we ma,y 
assume » certainly for a central rotor position «■ that in the 
tangential direction there is no pressure loss, but merely a 
decrease in that velocity eoaiponent. Consequently, the energy 
equation (3.I6) can be divided into two partial equations, for 
the pressure loss in normal direction and the velocity decrease 
in tangential direction. 

" "I ^2 "* f sin^o(j(l-^j), 

‘ ' ( 3.46 ) 

0 «• I w| oos^ocj (l + ) “I oos^a'i^ (1 - 1| ). 
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®li© ssoooM ©qimtlQn yield© bii© same 2 *©siilt as 
equation ( 3 . 17 ) if ©irbemal eompnessiv© fonees 


til© iaipuls© 
ana neglected 


9 


i.©.s if V/© set 


1 - I’j • i/T^ , 

■5 ■‘‘I’e ® ( ^ + {a • 


( 3.47 ) 


Shese relationships can also be applied to a plain clearance, 
and for small clearance v/idths = i.e., large pressure loss 
coefficients = yields more precise results, there. Increasing 
the number of support locations, for the same clearance length, 
the flow paths in a control space will be shorter, and hence the 
pressure loss coefficients ivill be smaller. Test calculations 
shov; that using equation (3»^7)» the number of support locations 
has no noticeable effect provided ^ < 1 is maintained. In 
contrast, this is not valid for equation (3*^-3)» obtained from a 
simplified, description, on a straight stream tube. Nevertheless, 
it will lead to the same result, provided the pressure loss 
coefficients are sufficiently small. 

3.3»2«3« Considerations on a labyrinth clearance /§1 


The support forces that cause a decrease in the torsional effect in 
a labyrinth have already been taken into consideration by impulse 
equation Since the afflusc from the labyrinth chamber to an 

immediately subsequent seal peak was asswiied to occur without 
losses, the support force acting on the entire control space (cf. 
Figure 3® 9) can be determined by a consideration of only the 
processes of velocity energy vortexing inside the chamber. If it is 


assumed that in a labyrinth chamber as that in Figure 3. 15 there is 
no pressure change due to the turbulence of velocity w.^, then the 
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Figure 3«15 Direction change in 
the labyrinth chamber 


residual velocity can be determined by the energy equation 

I wf (1-ii) - I v4 

independently of the pressure loss coefficient Eliminating the 
chamber velocity j the impulse equations in the tangential direction 
and perpendicularly to it will be ^ provided that for a concentric 
rotor position and constant chamber pressure we disregard the 
compressive forces on the control space's surfaces = 

til Wj [css Kj - cos OCg l/T^j ] - S cos 0's 0 , 

m Wj [sb sin sin ocg » 0 . 

The angle depends on the processes in the labyrinth chamber and 
can be determined precisely only from measured flow line courses o 
It is possible to determine -» depending on this angle “ the 
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®iipp@2*ting foK^a© W and 



I C0G 1-5 r 

an s5no<„frS|J. ' 

S » ihuj . f ^ = 2 tos(o<>j.-o(tj) +^'* "Ji) « 


( 3.48 5 


Sine© iaipulse equation (3»1?) already inoludes the assumption 
that the support fore© works against the velocity Wj at the seal 
peak, it is necessary to introduce a. correction factor y that 
v/ill take into account the redirection of flow inside the 


labyrinth chamber. With the definition of the support force by 
equation (3.18) and the direction established in aquation 
(3.17), using the impulse equation for the tangential direction, 
we obtain the iaipulse loss coefficients 


? 



v/here . 

0 Ceso<2 


( 3.49 ) 


According to Figure 3.1^, Y > 1 the stream tube would be 
redirected in a peripheral direction, for instance, through 
rotation of the channel walls for y < 1, larger afflux angles to 


the next seal peak would be attained. The latter is conceivable 
if the vortexing of the tangential coraponent were enhanced, for 
instance by building special devices, such as crossbars, into 
the labyrinth chambers. On the one hand, this would improve 
turbulence ivithin the labyrinth chamber. On the other hand, in 
this manner the velocity becomes only partially turbulent, which 
corresponds to redirection v/ithout pressure loss. This effect is 
more thoroughly s'uudxed m secGxcn 3*^, usXxig ^ caxouxa&xon 


example . 


For a total loss of the velocity energy C = 1, equation (3.49) 
always provides, with impulse equation (3*3.7) - independently of 
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the correction factor y - a flow timt is normal to th© nest seal 
peaks in contrast , without redirection and v/ithout velocity 


th© same flow ansi© as that at the 


turhulene© ( y ~ i » ? ” • 
previous peak will be repeated. Shus, the effect of th© 
redirection factor introduced increases with lower turbulence 
coefficients ?o For the usual seal construct ions » however, it 
should be Y « Is then the same impuls© loss coefficient is 
obtained as was found from energy considerat ions , in the 
previous section. 


Figure 3 • i ^ » below 


sshov/s the effect 


of th© redirection factor 


on a stream tube. For the same affluse angle and the same 
velocities at th© seal peaks, iv© calculated th© outflow angles 
and the magnitude and direction of support force S, from the 



Figure ^A6 Stream tube changes due to the 
constant y (c=0o3§ aj= 20°5 


above relationships. Due to redirection within th© chamber, the 
stream tube is changed considerably, however, as a first 
appro3simation we shall assume that the flow=line course =■ 
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speaking iralid only for y “ 3. = ©an b© ealawlatod as 
oiiulined in seetion 3.1. !• Du© •&© tho paF«?ial ireloeity 
'iJU 5 ?bulea©© , tli© out flow angle alviays beeoa©s las“g©r and depends 
on til© redireotion y« Shrougli tli® isipuls© equation* an angular 
eliang© is coupled to the stream tub© throtighput which, on the 
other hand, depends substantially on the flow angl® within th© 
elearane©, Shis relationship is described by th© disoharg© 

©qua t ion ( 3 • 29 ) • 


For a relatiTS affluss energy G„, th© round braolsets contain a 

iit 

eoiiibination of the pressure and impulse loss coefficients, which 
v/e simplified her© for the case of two seal peaks (n = 2, " 

5 ^= 1 )= 


D * 



Using th© impulse loss coefficient of equation (3«^J-9)» we obtain 
for the sas©, for instance, that the afflu 2 s energy and the 
torsional force can be used completely, for = 0, 


0 




Since < 1 must be, a redirection to the tangential direction 
Y > 1 xvith D < 0 iTOuld bring a throughput decrease, as a 
consequence, according to equation (3*29 ) 9 and viceversao The 
cause for this lies in the modification of the outflow angle 
\vhieh occurs (see Figure 3® 15)) due to the throughput -related 
redirection. For y = Is with D = 0 we obtain the special case 
already mentioned in section 3«2.6. Physically, an effect of the 
afflu3E torsion on the throughput is not meaningful, since it 
must be assumed that an additional suppo 2 ft force, acting on the 
stream tube du© to th© redirection, is balanced by 
correspondingly changed boundary conditions at the clearance 
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GSTsmaocD. a. si ans 
relationship 


gy ■'•5® equation C3ei!*9)» v;© assiM© the 


” |eE “ ^G6 ^ 0 

then for D - 0 (i»eo, no effect of the afflu 2 s torsion osi the 
throughput ) 5 w© can define the factor 




i / fc? 




for the labyrinth clearance considered in equation ( 3 “ 29 ). 

For the torsional effect supplied at the clearance is 

less than the flow before the entrance makes possible 5 hence the 
impulse coefficient 1 ms to be corrected by 


/H 



In contrast,, . the flow is redirected in an assial direction in 
the labyrinth chambers with Y^ls then it may be assumed tha.t at 
the entrance a smaller proportion of the velocity energy will be 
supplied to the tangential component o The equation 



takes this into account » The condition (3«50)i> above 9 must be 
satisfied for each of the stream tubes variable along the 


perimeters v/hich is readily possibles in connection with the 
iterative determination of the loss eoefficientsj also for the 
more general clearance form with radial entrance and exits 
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3® 3® 3' 


i5i22C 


iQS’eo at tSi© elmmiQl Qziit 


Past tko stat@5f ©gq, 1 ©leamneo, tliQ ta-sigoatial veloeit^r i'a tk© 
fii©s?idiasi eSmKiiiGl is relativQl^f largo o It is eonooi'valglQ tkat 
kos'o a iiizoi frietional for©© acts on tk© ol©arane© flox^Q Using 
Pln©g©l'e ©©nations {efo also [26]) tkis fore© ©an k© doseriked 


a Miz©d f riot ion 


5:3,©i©n’6 


4’ 


as r» •■•'■jjrio'Biovi 


tk© ©sit 


si^,rfa©© ar©a and tk© kinotie ©nsrgy of tk© difforsne© b©tiv©©n 
tk© loeal tangential irsloeiti©©!, 


Rtw “ \ ^A.a I ^ f. (3.51) 


Sk© diff©r©no©s ketwoen tk© tangential velooities is determined 
iteratively for ©aok stream tube, aeeording to section 3®^!-» 
simultansously ivitk tk© loss eoeffieientse To begin with# tke 
fore© s.bov® only aots to modify tk© ontflov/ angle, of tk© last 
eontrol space p sine© dn@ to tke fluid’s viscosity tk© frictional 
forces v&TQ disregarded in tke remaining control spaces. Ifitk /72 
a torsion-frss afflux and eccentric rotor positions an 
equalising flow ocourss v/kiek along the perimeter is symmetriea.1 
to tke position of tke narrowest clearance width. Taking tkis 
modified outflow angle into considerations larger mixed 
frictional forces in tke tangential direction are generated in 
front of tk© na.rroi'yest clearance than behind it. Due to the 
outflow angles corrected by these additional forces s tie 
pressure losses at the exit noiv changes as does the pressure 
distribution within tke radial clearance. Due to tkisp in spite 
of a torsion=free afflrac and without considering rotation, a 
transverse fore© now also acts perpendicularly to tke deviation 
(cfo section 3.6. ) 


Tkis simplified description at least perruits a qualitative 
evaluation of tk© mixed friction beyond tke cleara,nee§ according 
to [26] s for tke mixed friction coefficient we can use ^ 0.1. 

According to Table 3® 2s the outflovj angle caused by tlie mixed 
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f5?ieti©n 

depend El 


is a iime'bion o 


strei: 


^nigly ©n the ©zil5 


f 'fell© tlis^owglipirb 

loss imid hen© ©8 


SiaoQ it 
also 


in 

©El tilQ 


outflow angle )= at 
satisfied as tliat , 
eeans of whioli it 


tliQ essit a similas'^ eoEidition 


pl*e5 ‘tJO 1^0 


for the entrane©s in tiie psfevious seotioEi, "by 
is assumed that the th5?oughnut remains 


constant under the 


effect 


of missed friction. 


3 . 3 . 


h. 


Deteriiination of loss coefficieEits 
pressure distribution 


from a measured 


If the pressure 


distributions 


are Imovm from measurements s 


for a 


fsv/ essential reference planes, 
determine the loss coefficient 


then it should be possible to 
s from these in such a mamier that 


there is agreesient beUveen measurements and theory. Plowever, 
Imoivledg© of the course of the floiv lines is assumed, here, 
since smali changes in the torsional forces can decisively 
affect the pressure drop in the stream tube. If one makes 


assuBiptions rega.rding the impulse loss coefficients or their 
relationship to the pressure loss coefficients, then it Biust be 
possible to determine them from the measured pressure 


distributions, by means of the existing theory. So this end i?e 
can apply the same iterative calculation procedure ife use in 
section 3 '5s in which v/e calculate the loss coefficients frora /73 
Table 30 I 9 means of the pressure differences obtained from 
measure'ments. In this manner, we obtain the pressure loss 
coefficients for each measurement plane or, in conjimction with 
Table 3*2, also the impulse loss coefficients. The tivo 
quantities can be determined simultaneously only if, for 
instance, the flow lines are rendered visible so that the angle 
change becomes knoivn, and the losses become knovm by pressure 
measuremeiits. In a, labyrinth, however, the assiEiiiptions aiade on 
the composition of the flo?/ lines from flow angles are 
additionally of iEiport;ance, as are twose made about stream 


contraction. 
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^ 0 0 


Galoi^ilation ©f oho 




:©r*?e ©0 


a©'6a.ag 


©a. Q. E*Q'd01'* 


la addition to the ts‘angv©rs© fore©© diseussod below, there ar© 
also momenta that aet on th© eecentrioally located rotor, duo to 
th© olearanee flow. ®iey ar© generated especially by the 
coB^ressiv© forces on th© seal's front side, while moments due 
to differences in th© pressure drop are relatively small, within 
th© z^adial cl©arane©» Bseaus© of th© mode of construction of 
turbine steps, such moments can occur especially in the case of 
th© so-called recess steps, which have large frontal areas at 
the clearance entrance and exit. Since at these surfaces the 
pressure distribution still strongly depends on th© estimated 
values for th© b®nd=loss©s, and because th© effect of such 
moments presumably is only small, w© shall forego their 
numerical evaluation, in the context of this study. 


3.4,1. Transverse forces from the pressure distribution 


The general integrals for the forces due to the variable 
pressure p( 2 , ), acting along the perimeter and the length of a 
cylinder ivere indicated in equation (2,24). Since the calculated 
pressure distribution is given only at support points, an 
appropriate interpolation procedure must be found to solve these 
integrals. For th© longitudinal direction Zs as an approxiamtion 
the pressure calculated for the support point should be con- /?4 
stant across the control space's width, Az. Hence the forces 
acting on the rotor can be expressed as sums that result from 
the integra,tion of the pressure distributions p(^)j^ for each 
support plane, 



a 


zir 

cos f df) 
0 ' 

.1 sinip d?) 


C 3.S2 } 
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‘She pressure dist riM’S ions esa “be desorlbed by mesas of a 
Fourier series, to be dote^aiaed for the support points given in 
analogy to equation (3.19). If we select a dimensionless 
description of the pressure distribution, according to equation 
(3«3‘^0» then iv© have 


« i^n t A, 


(j T ri| C 0 S<P +A2,C.0S 2<j>+... A|j|C0Smj{) 
<■ s«n ^ Sm 2»p + ... sin vutf 

PiM-pA 


( 3.53 ) 


For j. support points, the maaber of coefficients is m - 
to raainta,in the functions pCil;)^ exactly as possible o 
here\'7ith integrate equation ( 3 . 52 ) alonp* the perimeter. 


(3/2)-!*1, 

If vm 
the 


Fourier terms of higher order cancel or . as is readily shown. 
Thus the compressive forces a.eting on a rotor are obtained as a 


sum of the first Fourier coefficients or respectively, of 
each support plane i of a radial cleranee, multiplied by the 
widths Az. of their control spaces, 


nS»ft 

Qi0 « -ApB % ir £ (Aii Aisj + ] 

Qzd * “Api t J { B|J^ t&i ¥ ) 

l"1 


( 3.S4 


Because of the different definitions for the control space in a 
labyrith and a plain clearance, xve again introduced the control 
quantity f^^, as in equation (3.1^). Because of the dimensionless 
description of the compressive forces, the above transverse 
forces are a, linear function of the total pressure drop and also 
of the shroud band radius, ivhile the effect of the total ividth, 
b = SAz. , is coupled to the amplitudes of the pressure 
distributions . 


For full rotor eccentricity e = s, in general the pressure 



dlfferene©© at th© ol©amaee ©nts’ane© ©an amount only to APg» 
Slaking a double amplitude 2A < 1 possible. Assuming a linoan 
pressure d©OK*©as© in axial direction we then obtain a fore© 


i; tr 


b e_ 

I ft 


{ 3.S5 ) 


Shis maximum attainable eoaipressive foroe = in which the 
relative ©eeentricity has been assumed to be linear - can be 

taken as a reference magnitude for the above transverse forces i 
from th© pressure distribution. 


3.4.2. Transverse forces from the clearance loss 


The local tangential force changes with the variable clearance 
throughput at the entrance to the rotor blading. The integral of 
the coraponents of this force along the perimeter yields the 
transverse forces acting on the rotorj as indicated in general 
in equation (2.10). This requires a knowledge of the local 
clearance lossj which can also be calculated from the ratio 





C 3.S6 3 


of the clearance throughput to the total mass flowo if we 
introduce a correction factor 

C 3-57 3 

The function g(y) takes into account that the decrease in the 
t^gential force included in ? is not directly proportional 
to the throughput ratio u* Any indication on this “ also /76 

available in [14] - can be made only based on efficiency 
measurements for the entire turbine step. It must be assumed » 
hov/ever® that the function g(y) is different for stator and 
rotor clearance losses, since the tangential foroe at the 
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blading is aff©ot©d in on© oas© by th© apps^oaohing mass flow and 
in til© other I by th© outgoing on©. 

With th© abov© definition of th© el©arano© loss, v;© obtain for 
th© integral of equation (2,10), bearing in mind th© quantities 
that ar© constant along th© perimeter, 




1^0 

Ms 

Itio 


2 « 

J 9(^) “ 




Up 

Qzs j cos(p dmsp . 


C 3.S8 ) 


The local throughput per are length can be described by means of 
a Fourier series 


*■ Co t C^ eos»f 't €>i eos2^ + - 
■ •# D, + a, S\n2.tf + .. 
A 

AtfH 


. Ctncostn 

Q^sin mq) 


( 3.59 ) 


whose coefficients are determined by the calculated throughputs 
per are length of the stream tubes. While at the rotor 
clearance we must use the channel ividths Aijj and the tangential 
angles if; of the first support location, at the stator it 

is appropriate to use values for the last support location, 

^nA k* equation (3»58) as a first approximation we set 

g(p5 = 1» then its integration ivill yield the transverse forces 
acting on the rotor, ^vith the higher-order Fourier coefficients 
canceling out, again? 


Q^s 

Qas 


Us 


Us 

ibo 


D,if, 

C,lr. 


( 3,60 ) 
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®kG eleamne© tlirouglipirb for tk© ©ntir© seal ean also b© 
obuained by integration of equation (3.59) » witli tk® sum of tk© 
tkrongkputs of all stream tubes of course yielding the same 
result 5 


21t ^ 

** j ®* 0^ *“ 1. *^ k 

0 k«i 


C 3.61 ) 


In principle, tke above transverse forces can also be determined 
as tke sum of tke individual forces at each support location, 
altkough this could lead to larger errors, depending on their 
number. For tke description of tke transverse forces due to 
clearance loss referred to, the ratio of the isentropic 
tangential force to the turbine step's throughput must be imowni 
for test calculations it may be assumed to be proportional to 
the tangential component of the afflux velocity. 


3.^-.3« Forces due to friction on tke rotor surface 


She frictional forces at tke rotor surface also cause a 
resulting transverse forces hov/ever, for the usual cases of 
calculation, it is t?/o orders of magnitude smaller than those 
mentioned previously. Equation (3*18) shov/s an expression for 
the frictional force acting on any arbitrary support point on 
the rotor surface, defined as positive in the direction contrary 
to that of rotation. If the components of this force are a.dded 
for all support points in the entrance, the radial clearance and 
the exit, then we obtain the resulting transverse forces 


Qir “ 

Qifl ® 




I (I S,. , .inf,. 4s, 


( 3«{53 } •* 
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iSiie f: 
■fchese 


fictional moment to i 7 hieh ths rotor is 
forces es.n "be calculated in a, simila,r 


sulsjected 'because of 
aianner , 
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rnC 

. X 

k*i ' l»J 



lift 





( 3.63 ) 


The efficiency decrease for the turbine step resulting thereby 
is negligibly small, at least for the rotor » 


3.5° Iterative solution to the problem 


The solution of the equations in the preceding sections requires 
the extensive iteration outlined in Figure which also 

follov/s the steps of the calculation programo The first inputs 
are the number of subdivisions n in the flow direction and in 
tangentis,! directions as well as the seal geosietrys including 
all radii 0 clearance widths 5 length and recess or cha.mber 
heights. Next come the pressures, velocities and angles for the 
main floiv ahead of and behind the seal, as well as the average 
density and kinematical viscosity. In addition, the input must 
include the bend loss coefficients and the constants for the 
calculation of the turbulence or vortexing coefficients in the 
labyrinth . 


The caleiils,tion itself starts with the estimated average flov/ 
angles a, established for each stream tube with its support 
points at peripheral locations channel v/idths Aip and surfaces 
Ao In a first iteration step, the loss coefficients along each 
flow line are determined. According to section 3.3., they depend 
on the flow paths and hence, both on the angles and the local 
Reynold's number. Once the loss coefficients have been 
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Figure 3.1? Iteration plan 
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d@t©min@d “ initially for oonstant thron^ipnt •=• with suffioiont 
precision, th© throughput of th© @tr@@M tubes is calculated and 
iteratively improved beyond the defined level 11, already taking 
into account th© changes in the flow angle, in accordance with 
section 3,2»5> Shis process is performed independently for all 
stream tubes k, as indicated in Figure 3»1?» above, by the 
subscript field. 


/I 


In a further block, the geometry of the stream tube is 
calculated in accordance with section 3»1 and subseq,uently the 
throughput is corrected (to level B82). Since through the channel 
width the stream tubes affect each other, the calculation must 
be performed ea,ch time for all stream tubes. If th© difference 
between the new throughput and old throughput of a stream 
tube is too large, then the loss coefficients must be 
recalculated, beyond level Mr 


Once the throughputs, angles, loss coefficients and surfaces of 
all stream tubes are precisely knoiTO, the pressure distribution 
is calculated. Once the Fourier coefficients have been obtained, 
this allows the determination of the compressive forces acting 
at each support location on the control space surfaces. These 
primarily change the geometry of the stream tubes? for this 
reason the iteration need not start from the beginning, but from 
index P. Subdividing throughput calculations into two partial 
iterations, starting at index Ml or respectively, M2, was ehosen 
only to save calculation time, since the loss coefficients 
appparently only change very little. 


For normal application eases the convergence is very good. 
However, for large torsional effects at the entrance and large 
eccentrioities, the flow angles vary considerably along the 
perimeter, causing significant changes in the channel widths. It 
is recommended not to perform further iterations directly with 
the new value, but with a weighted average of the old and the 
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nov! '^allies}. ^Jhil© tk© eM xvill 13 © attain©d 

soaewiiat more ®lov/ly, a div© 2 *goao© eau@©d “by ©koqssIv© ohaag©© 
in 'fek© iadi'vidiml values vdll b© avoided uheK'e'by, A siaiilar 
procedia^© aay also b© necessary in the iterative calculation of 
th© compressive forces, especially for large labyrinth chambers. 


At the end of the overall iteration, th© compressive forces Q. 
acting on the rotor can be calculated from the pressure 

distribution, and the forces Qg caused by th® variable 
taiigential force, from th© variations in the clearance 
throughput o Wo nexv initial values need be set to calculate 
further ©ccentricities . 


D 


/§k 

Based on the siniilarity of conditions, favorable initial values 
are obtained rather by increasing the eccentricity gradua-lly. 

Since th© cours© of the forces in general is only weakly 
non-linear, a few individual forces can cause an increase beyond 
the deviation as required for the vibration calculations in 
section Zol^ 


3, So Testing the calculation procedure using simple clearance 
f orms 

The essential parameters describing the clearance flow affected 
by torsional forces can be studied in Figure 3»18, using simple 
clearance forms o The clearance's diameter and length were chosen 
for their similarity to a test turbine, but without radial 
entrance or esiito Besides the radial clearance width and the 
chamber height h, the geometry of the seal remains unchanged and 
is supplemented only in Figure 3»37 with a radial entrance and 
eaeit . 

Table 3®8 clarifies the most important equations for the 
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smd tlie l3oimclai?3r ©oMitions 


oaleialation of loss eo©ffa.oiQa’fes, 

Qivon hy a aiaia flov/ maffootod fey tJie Gloamnes flow. A 
coast aat voXoeity eoapoa©at e is assuaed ahead of the 
cleajfaneOf in nozriial directioai which is fully used with the 


assumed eoeffieients Behind th© olearanee, only a normal 

velocity ©ssists; hov/ever, fee cause ~ 1» it has no effect on 
th© olsaranc© flow. By ehanging th© afflux angle j it is nov7 
possifele to desoi'dhe conditions such as they may occur ahead of 
th© rotor clearanc® of a turbine step. Under thee© conditions} 
for a constant static pressure drop, the available tots,l 
pressure drop Ap™ according to equation (3«i5) also rertiains 
constant. To calculate the transverse forces that act on the 
turbine rotor due to the clearance loss, we set the ratio 
between the • isentropic tangential force in equation (2.12) /83 

and the throughput equal to the velocity difference Ac^. It 
corresponds to the tangential component of the afflux velocity, 
since tiia tangential or periferal velocity disappears, beyond 
th© clearance 0 The relative afflujc energy can be modified by 


either th© angles cxv. or the pressure decrease APr 




For dimensionless descriptions the transverse forces acting 
on the x'otor are expressed in terms of the force Qg as the unit, 
according to eq_uation (3 '>55)'' This also reflects in a 
qualitativsly correct ma.nner the ratio between th© compressive 
forces Qg and. the forces Qg from the clearance loss. To clarify 
the results, we have in pa.rt also provided tl'ie difliensionless 
pressure distributions along the perimeter, at support 
locations 1, according to Figure 3. 18. Next to these pressu;:- 
distributions, the throughput of the stream tubes, 15^ ■= m^ym^ 
has also been recorded, here referred to the total throughput 
according to equation (3o2?). 


The examples of calculations, belmv, are intended to illustrate 
the effect ©f the afflux conditions as v/ell as that of the 
pressure and iBipulse loss coefficients on the transverse forces 
acting on the rotor. Initially, in a dimensioned description, we 



iavdstlgat© til© <l©poM©aoo of tlie f 05 ?a©s on tli© or. ?*snts"ieity aad 
til© pj?©s!3iir© cliff©r©ne©j sinae tiiGjr ar© inelMoc? in tli© 
s*©feK*ene© quantity;. For a piirel^r axial afflux ^ Figur© 3»i9 ©liov/© 
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Figure 3=19 Plain olearanosp 
s-lramp a-p .==90 s Ap-^^SOaibarp ^^=0 

c„-l 



Figure 3»20 Labyrintlis S"irrMp 
H=Akibis ag -90 p APg "50 aibarp 
=0p C^”1 d ?|(ai3_==9)9 yd) ~1 


the variation of the restoring force 
Figure 2o5) across the eccentricity. 


"°^1D (ef. definition in 
Since the eosipressive force 


Qpjj vanishes perpendicularly to the deviation due to a 


syaimetrical pressure distribution, 


the forces Qg 


ca,used by the 


variable 
anulled , 


tangential force at the turbine blading must be 
b©ea,use of the assuEiption regarding the velocity 


increase Ac^ in fable 3,60 $he approximate integration of the 
pressure variation iii equation (2 .,25) yields a liae®,r dependence 
on the ec'tentricity 0 which in comparison to complete 


calc t ions 


is valid only for small deviations. 


Due to the 


re V'; 1 v 


clearance length the effect 


of a balancing flow 
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is n®gligi'ble 5 in comparison to th© calculation of th© pressure 
variation from equation (3 = 35 ) o for which a purely ascial flov/ 
was assumed 0 Taking into account a variable tube friction 
coefficient X(R©)p larger restoring forces are obtained „ whose /84 
dependence on the eccentricity is also non-linear o The same 
calculations were performed for the labyrinth clearance „ in 
Figure 3 = 20i, above o Due to th© incomplete velocity vortaxing in 
the labyrinth chambers a restoring force “Q-j^ acts on the rotors 
which varies nearly linearly across the eccentricity o Here too 5 
the effect of the equalizing flov/ is only slightly larger than 
for the plain clearance « In accordance with equation (3»35)8 
taking into consideration a contraction coefficient y that 
varies with the local clearance v/idth does have en effects for a 
purely axial flow 5 since the surfaces at all seal peaks are 
reduced in the sarae manner o 
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Figure 3 ”21 Plain clearances s = 1 mm 
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For a flov; aff®ct©d by torsional foreesj Figure 3 = 21 ^ above^ 

shows that there are transverse forces acting on the rotor o in 

the direction of the deviation!, and perpendicularly to itp 

Qgo The cause for these forces can be seen froni the pressure 

variationsp next to the figure « For a central positioningp e/s - 

0 and the throughput and the pressure along the perimeter are 

constant 0 The strongest pressure decrease occurs betiveen the 

conditions ahead of the clearance, p'*' = 1, and the first support 

locations p^ (cf, also Figure 3»l8)s v/hil? pressure losses along 

the clearance and to the end states p^ = 0, are relatively 

small 0 For an eccentric rotor positions this pressure decrease 

varies in an inversely proportional manner with the local 

clearance v/idth, which is smallest for ip = Oo Due to the 

torsional effect, there is a spiral flov/ in the clearance, which 

results in variable cross-sections along the stream tubes, and 

causes the characteristic pressure maximum ahead of the 

narrov;est clearance, as seen in the direction of torsion. As a /86 

consequence, compressive forces arise, Q 2 J)P perpendicular to the 

deviation, v/hich for a plain clearance are of the same magnitude 

as the restoring forces The local throughput (broken line) 

modifies the tangential forces at the rotoi blading, whose 

resultant generates a force , perpendicular to the 

s 

deviation, and the much smaller force Q- , v/hich acts in the 

is 

direction of the deviation. If the variations in the throughput 
were symmetrical around ip = 0, then according to equation ( 2 , 15 ) 
that force would have to vanish. Across the eccentricity, the 
transverse forces are linear, ivith a small departure at 
This is also true of the flow with torsional effect, in a 
labyrinth clearance. Because cf this it is possible to build 
dimensionless increases of the forces across the deviaition 
using the reference quantity of equation (3-55)i) and in 
accordance with Table 3®^“ all further calculations, this 
increase is determined from this force only for a relative 
eccentricity e/ s = 0,6 , We thereby also obtain an approximately 
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average value for a minimally linearised increase o 


The effect of the radial clearance v/idth is shown •= with a view 
towards a vibration calculation ■= on the basis of an increase 
Q/e of the force Qj assumed linear across the deviation £o 
Qualitatively this also corresponds to the representation of 
transverse forces for equal rotor eccentricity o As Figure 3 a 22 



Figure 3»22 Plain clearance 


shows 0 for a plain clearance the forces from the pressure 
distribution increase steeply for smaller radial clearance 
v/idthso The amplitudes of the pressure variations depicted are 
approximately equal;, since they are valid, in each case, for 
half the maximum eccentricity, e/s = 0o5" For smaller clearance 
widths, the pressure loss increases along the clearance, thus 
increasing the pressure level at the entrance, 
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As already sjcplained in section 3»3«2o2op in a plain clearance 
and luith decreases in the radial clearance v;idtho under some 
circumstances inadmissibly large values can arise for the 
pressure loss coefficient 5 in a control space 0 Since the 
pressure losses depend on the flow paths within the control 
space p if necessary the number of support locations in the flov/ 
direction may have to be increased,. Using equation ( 3 «^ 7 )!> 'the 
calculations here were performed for all clearance v/idths with 
nS = 11 support points 0 However p no deviations worth mentioning /87 
v/ere observed for radial clearance widths s>Oo 3 mmp even for nS 
= 3. The determination of impulse loss coefficients from 
sectionally straight stream tubes in accordance with equation 
(3-^3) is also valid only for correspondingly short flow paths » 
However p the use of this relationship yields the same result - 
for a sufficiently large number of support locations - and was 
applicable 5 in particular and without restrictionsp to nS = 3 
and s = 1 mmo 

For labyrinth clearances the number of support locations is 
determined by the number of seal peaks 0 In additionp only 
pressure loss coefficients ?<1 can occur. As Figure 3»23 showsp 
the increase in the forces due to the pressure distributionp 
with small radial clearance widths is very much less than with a 
plain clearance. While the exciting force Q2^ is independent of the /88 

radial clearance width, in accordance with equation 2 . 15 , this results 
in a change based on the variable pressure loss coeficient. By disregarding 
the flow contraction ff = I, the exciting force O25 is proportionally increased 
to the throughput. However, the pressure distribution forces are not 
similarly altered, as the cross-sections along the flow tube become more 
equal because of the flow contraction (s) , which is dependent on tlie 
local clearance width (c.p. fig. 3 . 8 ). A flow contraction which is constant 
along the perimeter would therefore cause an increase in the exciting force 
from the pressure distributor. 
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Radia"' clearace 
Fig ure 5,g3 lRByrinth h=4 hm 

ag=20° Apg=50 mbor (e/s=0,5) 

^E=0 ;*=' Llm,=9),v,=l 
:^nr:}u=i i:~~)uis> 

For a flow affected by torsional forces, the influence of a pressure drop 
on a plain clearance is shown in fig. 3.24, while fig. 3.25 shows the in- 
fluence on a labyrinth. The lateral forces from the asymmetrical tan- 
gential force at the turbine blading indicate a parabola -shaped curve, as 
they are solely dependent on the throughput (which represents the square 
root of the pressure drop) , based on the assumptions of tangential force 
according to table 3.6 (c.p. equation 3.27). Calculation of the exciting 
force Q 2 g with a calculation of clearance loss according to Traupel 
(equation 2 . 15) , results in reduced forces , as complete vorticity at both 
seal pealcs is assumed. 


99a 





Ee=0 Ea=' 




Assijming a purely axial flowj the restoring force s ea '.w a nearly /89 

according to the pressure distribution shown in eqmtion '5. “v .vi th a pressure 
drop of ^ Pg- This dependence also seems to be valid for a flow affected by 
torsional forces, at least for a labyrinth. However, the exciting force Q 2 p 
at the two clearance dividers is non-linearly dependent on the pressure 
loss, as the relative afflux energy changes simultaneously, for which the 
significant influence is shown in the following figures. The pressure 
dropi! pg = 50 mbar, which could be chosen from the weak reaction stage at 
the rotor clearance, was chosen as the initial value for further sample 
computations . 

As all parameters contained in the reference values were changed, an 
additional dimensionless representation according to table 3.6 may 
be selected, in which the forces (f and the throughput m^ is plotted on the 
same scale. The relative afflux energy was varied (see Figure 3.26) with 
the angle ©<'g ahead of the clearance, and with the pressure dropjpg, starting 
from the average value marked. The pressure distributions depicted next to 
the figure show the effect of the 
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Figure 3.26 Plain clearance „ s = 1 mm 


afflujc angle j due to which the location of the pressure maximum 
is displaced towards the direction of torsion^ with a 
simultaneous increase in amplitude « 

The local throughput changes in the same manner p with increasing 
torsional effects at the entrance o Because of the spiral flow 
within the clearance p the flow paths become longer p which leads 
to a slight decrease in the throughput of reference p m^^p as a 
function of the relative afflux energyp C^o Variations in the /gl 
angle or the pressure gradient do not yield the same results 
here either p for the transverse forcesp since Reynold's number 
and hence the friction coefficient X changes somewhat with the 
throughput's absolute quantity o However p within the range 
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depicted,, the relative afflux energy Cg is characteristic for 
the generation of the corapressive forces that act on the rotor o 
On the other hand 5 this also follows from the calculation of 
floViT angles using equation (3"33)p which for a given geometry 
and pressure loss coefficients are only a function of Cgo It 
follows that exciting forces from the pressure distribution 
are generated only in turbine steps with small reactiono Thus 5 
besides a high tangential velocity ahead of the clearances a 
relatively low pressure gradient must exist simultaneously 0 In 
contrasts the restoring forces due to the changed friction 

values depend only on the relative afflux energy 0 

The exciting force Q, , acting in the same direction as the 
compressive force Q2p» also increases with the relative afflux 
energy s on the basis of our assumptions regarding the tangential 
force at the turbine blading, in accordance with Table 3 » 6 o As 
the angle a^, varies, this increase is caused by the velocity 
difference, Ac , for an approximately constant clearance losso 
In contrast, as the pressure difference Ap^ changes, the 
velocity difference Ac^ remains constant, while the clearance 
throughput changes 0 Obviously, both effects cause the force 
- which acts in the direction of the deviation - to follow the 
same course 0 

In a similar representation, Figure 3»27 shows the effect of the 
afflux energy in a labyrinth clearance o Here, variations in the 
afflux angle and the pressure gradient yield exactly the same 
results, since the loss coefficients depend only on the 
geometryo Thus, at the clearance, for different throughputs but 
the same relative afflux energy, the flow angles will be the 
sameo Due to stream contraction, the throughput (for the same 

e ^ 

reference quantity m_) is smaller than it is for a plai,;* 

JJ 

clearance, and hence the transverse forces generated, 
smaller 0 Since according to equation ( 3 »^ 0 ) the turbulence / 
coefficients C depends on the flow angle, the thi’oughput becomes 
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Constant flow angle ApB=50mbar 



Figure 3»27 Labyrinth, s=l mm 
H=4mm, 5 ^= 0 , 5^=lp ?;^(m^=9)o 
Yi = l 



without pressure 
equalizing flow 

Figure 3® 28 Labyrinth ^ s=l mm 
H=4mm„ Cg=0p c^=l , (^^^=9) » 

y(s) 


smaller with C^o Compared to the restoring force -Q^qp the 
excitation force from the pressure distribution is 
relatively small.. However, Figure 3 <>28 shows that this changes 
considerably, if it is assumed that no compressive forces are 
acting on the surfaces of the control spaces (cfo Figure 3<>10)o 
Because of the chamber height, these forces are much larger for 
a labyrinth than for a plain clearance » They act at the 
perimeter, in the direction of the lower pressure and modify the 
stream tubes in such a way that not only does a, pressure 
equalization take place, but there also is a phase displacement o 
This effect shall be further elucidated with the variation of 
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chamber height o Because of the unordered floiv in the labyrinth 
chamber p we could assume that the compressive forces acting on 
the stream tubes are much smaller than indicated by equation /£2 
( 3 » 21 )o On the other handp it is conceivable that the effect is 
balanced by a corresponding change in the local loss co- 
officients. If we neglect the pressure equalisation flov/ in the 
labyrinth chamber p as above j the transverse forces acting on the 
rotor are of approximately the same magnitude as in the plain 
clearance p in accordance with Figure 3=260 

Figure 3»28 shov/s a strongly simplified calculations under the 
assumption that all flow angles in the clearance are constant and 
correspond to the afflux velocity angle = This allows a closed 
solution for the pressure variation with equation (3"3^0o since 
the local surfaces for all stream tubes are given = The 
excitation force from the pressure distributions increases 

considerably as the angle decreases = In factp the maximum of the 
pressure distribution could be displaced so farp by prescribing 
the flow angles that the force will act in the direction of 
the deviation 0 Such a calculation could be very illuminatingp 
qualitatively (cf= Rosenberg [ 28 ])= floweverp the excitation 
forces it yields from the pressure distributions 'too 

large unless one assumes an average flow angle that is 
considerably larger than the afflux velocity angle 0 It is 
possible to perform an approximated calculation for this average 
flow angles as a furiction of the relative afflux energy p in 
accordance with equation (3»33)p for instance assuming constant 
flow cross-sections and loss coefficients 0 



Flow-line course in a rolled out (to scale) 
seal clearance (Figure 3 » 18 ) 

without pressure equalization flow 

with pressure equalization flow for h=12mm 
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In Figure 3“29 p above 5 the height of the labyrinth chamber was 
varied (cfo Figure 3°l8)o While the forces related to the 
clearance lossp were only marginally affect edp especially 
the excitation force p Q 2 B'’ ‘‘'^® p 2 *essure distributions 

decreased very markedly with increasing chamber height 0 An 
explanation for this arises out of the course of the flow-lines 
in a rolled-out 5 true to scale representation of the clearance 
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as that shovm at the bottora of page 102 5 above 0 for a large 
chamber heights in comparison to a calculation v/ithout pressure 
equalizing flov/o Due to the lateral forces acting on the free 
portion of the control spaces s the flov/ lines are deflected 
someivhat more strongly just ahead of the narrowest clearance 
location ( = 0) 0 Based on this minimal variation in the outflow 
angle and the channel width of the stream tubes p the pressure /£5 
variations in the labyrinth chamber are substantially affected p 
since the cross=sections normal to the direction of flow appear 
squared in the energy equationo Viewed in the tangential 
direct ion p ahead of the narrowest clearance are generated 
smaller exit velocities p and hence smaller pressure losses p 
whereby the characteristic pressure maximum in the labyrinth 
chamber is weakened with increasing chamber height 0 For very 
large chamber heights p even negative excitation forces may 
occuro The calculation with a chamber height h = 0 corresponds 
to the same assumptions made for a plain clearances taking the 
pressure equalizing flow into consideration changes the 
transverse forces by a maximum of 


In subsequent test calculationsp we shall now investigate the 
effect of changes in the loss coefficients » Starting from the 
standard casep Figures 3 <>30 3»31s belowp show the variation 

in entrance and exit losses for both clearance forms o With 
increasing entrance loss Cgs "the throughput and the transverse /£6 
forces acting on the rotor become smaller » Maintaining the exit 
loss < 1 p the static pressure at the end of the clearance is 
lowered p since a part of the velocity energy is used as pressure 
againo This increases the throughput p which may become larger 
than the reference throughput (m^>l)p since the latter is 
calculated only for a loss = lo With increasing throughput p 
the average flow angles in the clearance also become larger (cf» 
equation (300))° ThiSp in turnp displaces the maximum in the 
pressure distributionp peripherally to the location of the 
narrowest clearances causing an increase in the restoration 
force “Q. 


‘ID 


5 while the excitation force becomes smaller.) 
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Figure 3»30 Plain clearance „ 
s = 1 mnij ^2 ~ A (Re) 



Figure 3»31 Labyri.nth clearance 
s=linms H=4mms c|=lo6p 

= Ip y(i) 


Assuming similar effects as in Figure 3->29 p this could become 
negative for the labyrinth <> 

/2Z 

Figure 3°32p belowp shows the effect of the friction coefficient 
X, which in the normal case is calculated from the local 
Reynold's number ^ For the restoring force a maximum is obtained » 
which also occurs with a purely axial flowo At low friction, the 
throughput increases, and at A = 0 it takes the value of the 
reference quantity, m^ = 1, in v/hich only the exit loss was 
taken into considerationo The pressure distributions choseii show 
very clearly that at lower friction coefficients o the decrease 
in the torsional effect in the clearance also is smaller, 
whereby the compressive force Q2j^ increases □ In Figure 3 0311 
below, the pressure loss or, respectively, the impulse loss 
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3o32 Plain Llearance 

Cp“Op 



Figure 3«33 Plain clearance ^ s=lmm 







coeffieients ar® ealculatad from ©quationii (3o3?)o 
respectively 0 (3'>'*^’5)i> to investigate “*h© effect of a rotating 
channel mllo The throughput has a fiat majdrnum at a tangential 
velocity of 60 ni/sj which can he attributed to changes in the 
loss coefficients caused by the friction path (cfo Figures JolZ 
and 3olJ)o Correspondingly ,> the compressive forces acting on the 
rotor are affected in a manner similar to the effect caused by a 
change in the friction coefficient 



Figure 3»3^ Labyrinths s=l mmo H=4mm 


21 «o. 0 


1«o 

p° 



In Figure 3"3^o above „ a loss coefficient Cj constant along the 
perimeter was assumed t at the first seal pealCo In contrast to 
purely axial flowj because of the variable flov/ cross-sections « 
different pressure distributions occur within the clearance o For 
a flov/ v/ithout loss ~ 0) the throughput attains its limiting 
value j = I5 since the stream contraction was disregarded 5 
hereo While the restoration force disappears completely for 

= 1 b the excitation energy retains a finite value 0 
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Hov/©v©rp this value can inereas© no further v/ith larger numbers 
of ssal peaks p sine® du® to oomplste turbulence p the floxv with 
an impuls® loss coefficient of f ~ 1 runs normal to the second 
ssal peak ■= io®op purely axial flov/o The pressure levsl in the 
labyrinth chamber is significantly affected by the loss 
coefficient and attains the values p approximately p that would be 
expected for purely axial flow (cf. equation (3o35))” Assuming a 
pressure loss coefficient constant along the perimeters the 
excitation force from the pressure d: strijution is larger, here, 
than the restoring force, which also substantially changes the 
relations in Figure 3-27 or 3*28, respectively o 

For a variable pressure loss coefficient (from equation 
(3»^*0))p Figure 3 = 35 shows the effect of the coefficient m^ o 


ca 



Figure 3»35 ijabyrinth, s=lramp H=4mm 


It is inveresely proportional to the loss coefficient, as can 
already be seen from the behavior of the throughput o In a manner 
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similar to that of th© friction coefficient for the plain / 

clearance (Figure 3 = 32)0 the restoration force has a 

maximum that can be displaced v;ith changing afflux torsi on » Due 
to the pressure equalizing flov/p the excitation force is 
considerably smaller than for the plain clearance 0 


In Figure 3 = 3^5 belov/p the factor y v/as modified 5 which in 



equation ( 3 =^ 9 ) describes the relationship between the pressure 
loss and the impulse loss in the labyrinth chamber 0 Starting 
from the standard value y = I 5 it was assumed that the flov/p for 
instance due to the shaft’s rotations suffers a small decrease 
in the torsion (y> 1 )= The pressure level in the chamber is 
thereby raised p since the pressure loss at the second seal peak 
increases because of the increased velocity (cfe Figure 3ol6), 


109 



Beoauss of th© changed pressur® distributionp somev/hat larger 
excitation forces Qgj^ are generated o If the velocity is 
redirected 5 in the labyrinth chamber j v/ithout further pressure 
losses (y< 1 )o then only a little of the velocity energy can be 


destroyed at the clearance exit. This causes a lowering of the 
pressure level in the labyrinth charaberp v/ith a simultaneous 
displacement in the location of the pressure maximum. A strong 
reduction in the torsional effect on the flov/p caused for 
instance by crossbars within the labyrinth charaberp leads to a 
considerable reduction in the excitation force Qp^o taking into 
consideration the lateral corapressive forces acting on the free 
portions of the control spaces. This force can even adopt 
negative values p as has been observed during tests v/ith these 
fittings [30]” The calculations above were perforraedp as a first 
approximationp with the composition of flow-linos frora the /lOO 
tangents to the flow angles » v/ith a weighting factor g = 0 o 5 n 
described in equation Changes in this factorp even in 

extreme cases p caused no deviation v/orth mentioning in the 
remaining calculation examples. Taking into account the pressure 
equalizing flow in the application case described (y< 1 ) is also 
only approximately valid p since by adding features in the 
chamber p additional corapressive forces act on the stream tube. 



Figure 3 "37 Plain clearance with radial 
entrance and exit (cf. Figure 3 , 18 ) 
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Figure 3o37o above j, shov/s two further clearance forms j, in which 
the flow's radius remained as in Figure 3® IS- Form A 
appro3cimately corresponds to the shroud band on a rotors while /lOl 
form B could b© used as the stator clearance of a turbine step in 
chamber constructiono The effect of the radial entrance in form 
A is shovm in Figure 3o3^> belowo With a modification to the 
entrance clearance = to correspond to the ascial clearance 
betv/een stator and rotor blading = the forces due to the 
pressure distributions Qj^o remain nearly constant j, even though 
very large pressure differences occur in the entrance clearance „ 
due to the higher velocities » In contrasty the forces Q 
decrease p due to the lower clearance loss in a smaller a^cial 
clearance p something confirmed also by measurements to be 
discussed later. 



Figure 3® 38 Clearance as in Fig. 3.37(A) 
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Figure 3=39 Clearance as in Pigo 3 “37(B) 


/102 

In the case of clearance form the clearance of the radial 
exit “ which in a turbine corresponds to the axial clearance - 
was changed simultaneously with the entrance clearance „ as shown 
in Figure 3 “39c above o The velocity difference which 

according to Table 3“^ determines the transverse forces arising 
out of the rotor blading’s variable tangential forces is caused 
here by the tangential velocity behind the clearance o If the 
entrance or exit clearance is reduced g the clearance throughput 
will be lowered and this will lead to a decrease in the 


excitation force 
axisg the forces Q. 


Because the afflux is parallel to the 
and Qgp vanish g since the pressure 


equalizing flow occurs in peripheral or tangential directions 


symmetrically to the deviation. The transverse force acts in 
the direction of the deviation for small exit clearances » since 


here the pressure distribution depends very strongly on the exit 


velocity. Assuming even larger pressure loss coefficients at 
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the elbov/Sj a decrease in the restoration force occurs j v/hich 
eventually disappears completely „ if only pressure losses occur 
that are independent of the local clearance width o 

While we have assumed,, so far^ that the tangential velocity of 
the main flov/ behind the clearance has no effect on the 
clearance flowp in Figure 3»^05 belowp we took into account a 



Parameters as in Figure 3»39 
s-E = s^ = 4 mm 

Figure 3-40 Clearance as in Figo 3° 37(B) 


mixed friction force in the sense of section 3o3»4o It increases 
proportionally to the coefficient and causes an asymmetry in 
the clearance flowo 


/103 


Because of itp the force and the excitation force are 

additionally generated from the pressure distributiono Since 
under the assumptions stated the mixed friction force affects 
the clearance flow significantly only at the exitp its effect is 
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relatively small = From this we may conclude that the esscitation 
forces related to the pressure distribution are mainly generated 
at the rotor clearance of a turbine stepo At the stator 
clearance 5 besides the excitation forces caused by the clearance 
lossp transverse forces are generated by the pressure 
distribution that may act against or even in the direction of 
the deviationp if we take a radial entrance or exit into 
account » 
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^ EXPERIMEML DETERHI .NATION OP THE TRANSVERSE FORCES ACTING ON 
THE ROTOR 


In agreement v/ith their definition„ the individual coefficients 
of the 'matrix proportional to the deviation (see equation 
2o3) can b© determined by means of the measurement of the 
trans^verse forces acting at the turbi'ne rotor 5 as a function of 
its static displacement 0 In contrasts the forces proportional to 
the velocity p in matrix can be determineu only from a 
predetermined vibratory rotor movement.. Kinetic tests (cfo 
Wohlrabp [30]) make it possible to study the simultaneous effect 
of all the coefficients of the matrices proportional to the 
deviation and to the velocity p from the fading of the aisturbing 
vibration for known system darapingo This study focuses only on 
the static methods which has the advantage - besides the 
separate determination of the forces proportional to the 
deviation = that the effect of a pressure distribution inside 
the seal clearance is relatively easy to investigate 0 


4 olo Test assembly 


4 ololo Installation construction 

The transverse forces caused by thfe flov/ medium at a turbine 
stop are many times smaller than the rotor's inertiap 
representing at most 10?^ of the tangential force acting on the 
rotor o Since the clearances between rotor and housing are very 
smallp special problems arise in the exact positioning of the 
rotoro The force measuring apparatus is arranged around the 
rotor in its stationary positions an eccentric positioning can 
be simulated by moving the housingo 




Figure 4ol Test assembly schematic 

Figure aboves shows a schematic representation of the 

arrangement chosen o Prom the fixed piping a^ the fluid (air) 
streams into the housing b« This is placed on a twofold support 
c in such a manner that it can be displaced both in axial and in 
horizontal direction in relation to the stationary rotoro This 
makes it possible to vary the axial clearance between the /105 
bladings 5 and the eccentricity of the rotor with respect to the 
housingo In an expansion of this arrangement p the housing b can 
also be rotated in the horizontal plane,, making it possible to 
adjust the rotor's inclination with respect to the housingo The 
rotor d is held in the bearings whose retaining rings are 
fixed to the bar structure f . The latter connects the two 
bearing brackets and exits the housing without contacts 

By means of pressure gauges and Qgo which in principle 
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^ Tubing ^ Movable hou^sing Moving mechanism: ca Axial clearance cb Eccentricity 

cc Inclination cd Clamping arrangement ^ Rotor £ Bearings Rod arrangement £ Guide Ball 

rollers _i Axial prestressing j_ Coupling k Housing ^ Rubber grommet m Sealing for rod arrangemnt 
a Turbine step £ Sheetmetal guide £ Flov/ guide £ Spent air £ Dial gauge (housing position) 
s Longitudinal guide (Ball boxes) _t Balance beam (Q2) £ Dynamometer adjustment 3 C Tachometer 

^ Deviating roller y Position plane of housing £ Measuring surfaces on rotor 


Figure 4.2 Test installation 
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(to Figure 4o2j above) 


are very rigid spring dynamometersj the bar structure is 
horizontally and vertically fixed. Moraents acting through the 
bearing on the structure „ on the one hand are picked up by the 
guide gs which is freely movable over ball rollers h. On the 
other hands cables it, prestressed at location prevent 
rotation in the horizontal plane. In the coordinate system 
referred to the test apparatus t for a turbine rotating to the 
rights the excitation force is indicated at gauge Qg and the 
restoration force at Through a bending coupling the rotor 
d is connected to a vibration generators at which the tangential 
force U,, is measured. Force measurements at Q, and Qo are not 

a Id. 

affected by the couplings since its spring rigidity is very lov/. 

In additions such side effects are eliminated by calibration of 

the pressure gauges under actual conditions. The frictional 
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moment of the turbine bearinjgs is transferred to the guide £ via 
the rod structrre and hence can affect only the tangential fore© 


U 


B 


5 


but not the force measurements at 


and Qgo 


Figure k-oZs abovsp shows the test assembly built p with the 
construction elements a through corresponding to the schematic 
on page 117 o In order to relieve the vertical pressure gauge 
from a part of the rotor weight,, a balance beam t provides a 
counterweights through two ball=bearings 5 at the rod arrangement 
f 0 The pressure gauges are fixed to the rod arrangement through 
cross guides of two ball boxes eacho The distance and the angle 
of the pressure gauges can be adjusted - by means of a 
tightening screw and three lifting screws u - in such a manner 
that all guide elements are perpendicular to each other., The 
longitudinal guides used permit a nearly frictionless rotor 
suspension 5 which is also free of play due to the prestressing 
applied o The central rotor d positioning with respect to the 
housing and that parallel to the axis» is established prior to 
assembly by means of four dial gauges and can be reproduced 
after assembly 0 
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4 olo 2 o Installation operation 

The test turbine was operated on an open circuit 0 with the air 
supply being provided through the compressed air pipeline of the 
Technical Iniversity of Munich's heating power plant 0 Figure 
4 o 3 s belowp shows a diagram of the plant, approximately to 
scale, with a schematic of the air supply lines o 

A two-stage cell compressor a transfers the air through filter b 
into pressure tank From there, a pipeline (nominal diameter 
80 ) leads through a restrictor d to the quick-acting gate valve 
e» After a long diffusor £ with subsequent straight ener, the air 
reaches the afflux tube (nominal diameter 230) to the turbine ho 
Together with the vibration generator k, they are mounted on a 
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Figure Arrangement of the installation 

baseplate., The spent air is blown into the machine room via a /l 10 
silencer o On one side of the control panel ra are mounted the 
essential control elements for the test turbine ho The pressures 
are recorded at a U-tube ivall n^ and can be fotographed with a 
camera from the measuring siteo The measurement amplifying 
and indicating instruments are located approximately at £o 
Shielded cables lead from there to a recording facility r 
equipped v/ith selector switch^ digital voltmeter and printero 


4olo3» Measurement instruments and test method 

The forces v/ere determined by means, of Hottinger pressure 
gauges 8 whose spring elements are outfitted with wire strain 
gauges o The table belov/ provides a review of the gauges used and 
their spring rigidity j the possible measurement ranges » v/hich 
are much larger than the forces actually measured at the 
facility. The precision is not affected hereby s since a 
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Porca Maasuramant range Spring rigidity 
20 kp 800 kp/cra 

Qg 10 kp 350 kp/cm 

Ug 50 kp 1600 kp/cm 

calibration is alv;ays performed in the range indicated by the 
prsstressingo Because of the relatively high pressure gauge 
spring rigidity a path change due to force interference can be 
disregarded 0 The wire strain gauges were balanced with a 
Wheatstone bridge and recorded on the indicator of a carrier 
frequency amplifier = Since due to the eccentric rotor vibration 
high frequency oscillations v/ith large amplitudes are 
superimposed on the signals to be measured j, the working range 
chosen for the amplifiers had to be very larger With a low=pass 
filter in line behind them and digital readout „ it was possible 
to attain the required measurement precision o 


P2,tot P2 «t 


POjtOt 

Meas. plane ya p^^gj-Meas .plane yb 

Pd A ^ ^ 



Figure UoU Teat turbine longitudinal section 
with measurement points 


Figure above 5 shows the pressure measurement points on a 

longitudinal section of the turbine 0 The static pressure p . 

O p S u 

and the total pressure p . .. ahead of the step are measured 
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with two Prandtl=tubss» arranged rotated by 90° o The pressures 
Po <,.<. and p« . behind the step are measured in the same /l 11 

manner 0 For control purposes, several drill holes pg are made 
on the inside of the outflov/ channel wallo For standing blades, 
large underpressures were observed at this flov/ guide, due to 
the high eicit velocity, which depending on the rotor's 
eccentricity, showed an uneven distribution along the perimeter » 

In order not to affect the measurement of the transverse forces 
with this additional effect, the flow guide v/as bolted to the 
housingo During some of the measurement series, an attempt was 
made to measure the pressure betiveen the stator and the 
rotor, in a ring chamber above the rotor's coveringo However, 
this v/as in part very imprecise, since the corresponding housing 
portions could not be sealed off in the manner desired o During 
tests with shroud band, the pressure distribution in the rotor 
clearance was measured using a maximum of 24 v/all holes, 
p^ 1 24' exact location and size is described in section 


All pressure measurements are performed with U-tubes, which are 
fotographed for fast and reliable recordingo Figure 4o5s below, 
shoivs such a fotograph, for instance one taken in measurement 
sequence 115 (© = 0"7 mm) o At the extreme left are the mercury 
columns for the restrictor pressures p Pg 2 pre-step /112 

pressure of the turbine o Next to these, to the right, 

appear the pressure differentials in the two Prandtl tubes ahead 
and behind the turbine step» The pressure (center) had to be 

measured with a mercury column, hereo To the right are collected 
the pressures p^ along the shroud band, whose single =arm 
manometer shows different levels for the labyrinth chambers K1 , 

K2 and K3o The barometric pressure is recorded from an aneroid 
barometer bo Asides from this method of evaluation, the most 
important pressures were also recorded by means of the pressure 
sensors of the data acquisition system used, whose code number 
(at bottom) appears on every fotograph » However, because the 
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tihi^ were wSMsmreA ty aeafts; of the- osrresfJoMiag: 

iroa-oonstaiifaji.: tli'#ri»e elements at fft# measiiFeffient orifice tp^ » 
as wall as ale.a^l Md: beJiied th«- torlsiae stea » %,■ » ~ foaci 

ori tiae: iigitaX f oltias fcer sf tlia: roeoriing Bfyipwsnl:* Siaae only 
siatleiiary afs- slaaisd wi-th- th® t-sst- iasXmliatioKt: 

all paffes o-oald: fee aeasared with orainany ilai gauges at: fixed 
liO'Using par-tsi Iieaes- the ipraoislon of aAJnsMjmrit -was- tisiMlly" d-f 
lass thmi Oi'Ol iiffit file- arrangement of tfte ae-asareXiSHf: aiits^ -tS' 
ieteriiine tlit rotor fositiOtt: as; eitlier coatral 'o;r pamlle'l- t'O- 
tbs; axisj:' witlt rsifect tbe housing in pairs:» displasai hy 
fO'*^ with raspsol: to eaoh other « W.b shown in figure 
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At the beginning of a test sequence, the static rotor was placed 
in the desired position with respect to the fisced housing parts, 
and the arrangement checked for freedom from friction. For a 
precise adjustment of the central position v/ith respect to the 
movable housing, airfloiv was established across the step without 
rotations hereby it was possible to reduce measurement errors 
caused by arrangement elasticity to a minimum , Once the dial 
gauges were removed, the rotor v/as brought to operating speed 
and the pressure gauges were calibrated using weights. It was 
then possible to change the housing position without modifying 
the operating conditions and record all forces, temperatures and 
pressures for each position. Thus a measurement sequence yields 
the variations in the transverse forces across the deviation and 
on the basis of repetitive measurements provides relatively 
precise average values for all turbine data. Within the 
precision of the measurements, these were independent of the 
eccentricity = or respectively, inclination » of the rotor. 


4,2 Test program 

The constant-pressure step represented in Figure 4,6, below, was 
used primarily to investigate the effect of the seal form on the 
clearance excitation forces. With the rotor suspension chosen, 
it was possible to measure simultaneously the forces acting in 
the direction of the deviation and perpendicularly to itf 
because of the test arrangement, the moments acting on the rotor 
were not included. Because of the non-variable compressor, the 
turbine throughput was fixed at approximately 0,4 kg/s. For an 
operating speed of 8,000 rpm, the pressure coefficient \|j was 5,2 
to 6 , 2 , The Reynold's numbers calculated from the afflux 
velocity and the blades' chord length was approximately 8 * 10 '', 

Prom measurements performed with standing blades - listed 
completely in one report [ 5 ] - it was possible to derive 
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Figure 4 q 6 Constant pressure Figure 4o7 Variations in the 
step geometry form of the seal 


essential trends in the clearance excitation forces p which also 
served as orientation for the remainder of the test prograrHo 
Because they are more applicable „ we shall mention here only 
tests with rotating bladingp which are in agreement with the / 
measurements above p once differences in operating conditions are 
taken into account o 


The seal geometry was varied by placing certain inserts in the 
turbine housings shown in Figure 4o7» above, for blading with 
shroud bando Table 4c 1, below, provides an overview of the 
sequentially numbered insert conditions, indicating the 
characteristic magnitudes o The variations in the radial 



TABLE 4ol Test program 


Insertion condition 

Geometry 


EB 

s* 

s** 

Form** 

21 


1.1 


1,0 

A 

22 w/o shroud-band 

1 ,S5 


1.5 

A 

23 1 

0,55 


0,5 

A 

] 





plain shroud 

1,1 


1,0 

B1 






B2 

26, 27 

two 

seal peaks 

1.1 

1,0 

C 

29 

three 

1,1 

1,0 

D 

seal peaks 


clearance widths v/ere performed only for blading without shroud 
band (form A) 5 for reasons of tolerances,, the stator clearance 
was slightly larger than that for the rotor» Since distinct 
pressure differences could already be observed along the rotor 
shroud band for standing bladingp special attention was paid to 
banded bladingo First „ the plain clearance was investigated with 
a simple gauge ring B1 5 for different operating speeds^ However, 
the pressure drop in axial direction could be measured more 
accurately only after adding a second gauge ring, B2» To 
conclude the test program, two labyrinth seals commonly used 
with turbines (C and D) were installed at the rotor 0 In 
addition, this made it possible to study the effect of rotor 
inclination with respect to the housing,. 


4o3« Measurement evaluation 
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4o3«l<> Turbine data and force measurements 

The turbine's mass flow is determined from the temperature at 
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the measurement orifices 
according to DIN--’ 1952 



and the pressures and Pg 2 ‘’ 


m - ot e m £ o"' f2^ (pg,--pB2). 


The pipe diameter is gi’ven as D = 80 mm and the opening ratio of 
the orifice as m = = 0o366. The density p can be 

calculated from the gas equation^ as a function of temperature 
at the orifice 0 The expansion coefficient 

throughput coefficient a(Rep m) were determined from DIN 1952 » 
taking into consideration the temperature dependence of the 
kinematical viscosity for Reynold's number « Because of the 
dependence of Reynold’s number on the mass flov/p iteration is 
necessary 0 


The turbine’s isentropic gradient is obtained from the pressures 
Po j >2 ■'^he temperature T^p by 



( 4.2 ) 


The force U-n at the brake's lever arm 1, is first recalculated 
B b 

using the calibration factor of the gauge p then the turbine 
step's tangential force 


u 


QJ 


JjPL 

dfa/2. 


( 4.3 ) 


is determined 0 With the angular velocity w at the rotor p we 
obtain the effective efficiency 


*|e 


U ^a/2. 
Ah'! ra 


( 4.4 ) 


DIN = Deutsche Industrie -Norm = German Industrial Standards 
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Correspondingly ^ th® i eternal efficiency can be calculated from 
the temperatirro, .nd (ahead of and behind the step) u 

p). a SslklSl ( 4.S ) 


From the tangential velocity u = w'r we can nov; calculate the 
reference quantity = from equation (2A2) = for the measured 
excitation forces o To determine the theoretical excitation 
forces requires knov/ledge of the step’s triangle of velocities « 
Prom the mass flow ra^t, the surface and the density we 
obtain the afflux velocity at the stator bladings 


m 




( 4.6 ) 


According to Traupel [14 ]d the blading efficiencies can be taken 
as Ti' = O 089 and n” = 0o82 for the existing stepo From the 
energy equation we obtain the stator exit velocity 


I 



( 4.7 J 


as a function of the stator blading's isentropic gradient 


Ali- en.[l -(!)*“ 1- 


I 4.8 ) 


An iterative solution is possible for the pressure p^ 5 if the 
stator's exit surface A^ sin has been very precisely 
determined 8 for instance through measurements of the individual 
channel widths o From the continuity equation we new calculate 
the velocity 



and "bhe density from th© gas ©quationj as a function of the 
turbin© condition 


i « fi. Jo ^ ^ 

§0 96 • 


The temperature can he calculated from the energy equation 
for th© stator bladingp 


“ fihg 1)' * I 


( 4.11 1 


From (4ol0) we now obtain the density ratio 


a . S i 


( 4.12 ■> 
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If herefrom we now calculate the velocity c^ from equation (4,9) 
and replaces this value in the energy equation (4,7) p then we 
obtainc with (4o8)„ the initial equation for the pressure p^^s 


1 

I 




[1 -(I) “ J 


+■ 


C 4.13 ) 


If we solve the right-hand side for p^^ t an estimated value p^ 
can be improved by iteration. We then obtain the step's 
percentage eaction 


t - 



( 4.14 ) 


However, the result of this calculation strongly depends on the 
surface sin and the stator's exit angle contained in it. 
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Disregarding the stator clearance loss in equation ho9 yields a 
much lower percentage reaction for the existing test steps while 
for loss-free flow (n’ - 1) there is an increase againo Since 
the two effects approximately cancel each others all further 
measurements “ considering the possibility of further sources of 
error = v;ere performed assuming n’ = 1 and 1 * 1 ^^ = Oo It was 
observed s in addition^ that even small changes in the throughput 
or the turbine's pressure gradient had a large effect on the 
percentage reactions following the above calculations » However ^ 
the latter is needed in the theoretical determination of the 
pressure distribution across the shroud band and hence also 
impairs the comparison between theory and measurements 0 

With the loiovm pressure pj^ we obtain the velocity c^ from 
equation (4o7) and the complete triangle of velocities 


=> Cj + - 2 mC, CotoC^ , 

14 >= rj* ) , 

C| • Wj + w '■^2 » 

o(, “ afc+oH 

* Wjeo*4j"U ' 


( 4.15 ) 


/m 

as a function of the tangential velocity and the angles 
indicated in Figure 4 60 We thus obtain the peripheral 
efficiency 

0 a CaCOtOta) ( 4.16 ) 


for the stepp assuming c^ = C 2 » Furthermore p the abscissa values 
for Figure 2=7 can be givens to facilitate calculation of the 
theoretical excitation constants from equation (2el5)o Since no 
equations are given in [14] for the clearance loss for plain 
clearances 9 as an approximation to it we used the value for a 
single seal peako 
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In each measurement sequence o the calibration factors for the 
pressure gauges were determined first o by linear regression on 
the electrical readings obtained for different calibration 
v/eightso After recalculation of the measured transverse forces „ 
their increases £ across the deviation or respectively 5 a 
were calculated by linear rogressiono For better applicability 
of these resultsj in conformance to the theory (see equation 
(2ol6))p these increases arc related to the ratio between the 
isentropic tangential force and the rotor "bucket" length 0 We 
thereby obtain the so-called clearance excitation and 
restoration coefficients through K2^p which describe the 
step's transverse forces 9 in a dimensionless formo The most 
important operating parameters and measurement results here 
discussed for the test turbines are compiled for selected 
measurement sequences in tables in the Appendixo 


4 o 3 “ 2 o Relationship between excitation forces and efficiency 
measurements 

The measurement of the clearance loss can be accomplished 
indirectly p via the efficiency where the 

peripheral efficiency is determined neglecting the additional 
losses due to extrapolation to radial clearance width s = Oo 
Since the excitation forces are very sensitively affected by 
efficiency variations p exact measurements are necessary p for 
many radial clearance widths 0 It hence seems appropriate to 
proceed inverselyp ioSop to establish comparisons to efficiency 
variations based on the many excitation force measurements 
available 0 It can be shown at least for blading without shroud 
band that the efficiency varies non-linearly with the radial 
clearance width 0 This can be taken into account by means of a 
polynomial of the clearance loss 


131 



( 4.17 ) 


fsp ® 4!^ + Qg S® 


as a function of the clearance width o If we assume ^ disregarding 
equalizing flov/Sp that this clearance loss is valid also for 
local radial clearance widths £ according to equation (2ol4)8 
then the integration of equation (2.10) yields the transverse 
forces acting on the rotors 

^<6 * I . ( 4.18 ) 

Q^8 “ f [(q<^2q 2S ^ S^) £ 4- I j. 

Since the clearance loss was assumed proportional to the local 
clearance widths the force in the direction of deviation 
vanishes. While for a linear equation a^ = = 0 - as we 

already established - the force is independent of the radial 
clearance width £p this effect can already be noticed for a 
parabolic equation = 0. However p it is only for higher-order 
clearance loss equations that an '"s““Shaped course of the 
excitation force across the eccentricity is obtained. 

If we assume that the excitation forces caused by pressure 
distributions can be disregarded for bladings without shroud 
bandp then the coefficients of equation (4.18) can be determined 
by the least squares method applied to the measured excitation 
forces Q« = Q .5 . Here we obtain a family of curves of the 
excitation forces as a function of radial clearance width £p 
across the eccentricity p whose coefficients determine the course 

of the clearance loss that would result as a consequence of 

sp 

the excitation forces. By further adapting this course to the 
measured efficiency 






( 4.19 ) 
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/121 

v;e can determine the peripheral efficiency and a factor 

v/hich the clearance loss from the excitation fore© 

sp 

measurements vrould have to be corrected. For Yq - 1 o there would 
be complete agreement (see Figure 4,8, belovir) between excitation 



Figure 4,8 Variation of the effciency across 
the radial clearance width 


force and efficiency measurements. For Yq < i » 'the measured 
excitation forces are larger than the forces that can be 
explained based on the clearance lossp which leads to the 
conclusion that there are additional effects p for instance due 
to the pressure distributionp at work. In general the factor 
will be Yq > 1 0 since the local clearance loss can be reduced by 
means of equalizing flows along the perimeter. In this case the 
measured excitation forces are smaller than those one would 
calculate from a precise course of the efficiency n^(s). 


It was assumed I in the evaluation of the results p that the 
effective efficiency determined from output measurements at the 
electric dynamometer are to be set equal to the internal 
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effi©i©noyo 


n., = n_o Th© m©asu5f©msnt of th© internal sffici©ney 

X © 


from the temperatures ahead of and behind the turbine is 
affected by measurement errors that are too large s for these 
considerationso But it is to be esspeeted tl-at the friction at 
the test turbine's ball bearings v/ill affect the efficiency 
in the same manner for all parameter changes., 


Especially important to the above correlation is the assumption 
that viewed across the perimeterj the smallest clearance 1< ss 
will also occur at the narrov/est clearance width o Because of the 
torsional affect on the flov/j according to section Jo6 this will 
not always be trueo As Figure belowj shows s a small phase 

displacement occurs also because the clearance loss occurring /122 



Figure Ao9 Flov/=line course 


at the stator reduces the tangential force only after an arc 

length bo Because of thiso in contrast to equation (Aol8) there 
s 

are transverse forces Q. active here that must be taken into 

Is 

account p in the correlation v/ith the efficiency variations o The 

measured restoration forces however point towards the fact 

that in spite of the relatively large flov/-line advance (ib = 

/ ^ 
2b^/d^^ =: 1 7 ° ) this effect is small. 

Ao3o3" Pressure distribution at the shroud band 

The purpose of the pressure measurements at the rotor seal 
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eleamnc© is th© dstarmination of the forces arising there o 
Integration of a prsssur® curve along the perimeter ip and in 
axial direction 2 , presupposes a tv/o^diraensional function p(so«|/) 
at the rotor surface 0 Assuming that the pressures measured at 
the housing v;all are constant across the local clearance widths 
this function can be determined from a limited number of 
measurement points 0 For the axial direction it is assumed = for 
a plain clearance « that in conformance to Figure 4ol0s beloiVj 




Figure 4olO Evaluation of measured pressure curves 


the measured pressures p^ correspond to the average values 
across the widths Azj for a labyrinth j instead p constant 
pressure may be assumed along the chamber width 0 On this base it 
is possible to calculate the compressive forces acting on the 
rotor as was done in section Since in the peripheral d’.rec- 

tion a maximum of only eight measurement points is available 
for a pressure curve p(4^)p it is described by means of three 
Fourier coefficients from equation (3»55)- This corresponds /123 
to the fitting of a sinusoidal function of variable amplitude „ 
with phase displacement with respect to the course of the local 
clearance widths it will be recorded v/ith the measurement values 
submitted later 0 

Based on the assumptions regarding the course of the pressure 
curves the forces determined from it could be subject to 
systematic error 0 Despite some deviations p the prescription of a 
sinusoidal pressure curve in the tangential or peripheral 



direction appears to lead to only sraall errorst sine© higher 
order Fourier coefficients (equation (3”5^!-)) drop out during 
integrationo In contrasts the result is severely affected by hov/ 
the i";idths Ag in Figure 4olO are establisheds since the pressure 
curve is non=linear in axial direction. It v/ould be possible in 
principle to fit a two-dimensional function p(zn 4 ') to the 
measurement points ^ but it v/ould lead to substantial calculation 
efforts 9 since the parameters v/ould no longer be independent of 
each other. This would make restrictions necessary such that no 
increased precision in the results could be expected a in 
comparison to the above evaluation. Qualitatively j moreover a 
measurements with only one measurement plane perpendicular to 
the carrying axis (cf. Figure 4.10b) yield useful results. 

Within a measurement sequence „ the pressure variations are 
recorded for each eccentricity and the forces Qj^ are computed o 
from whose course across the eccentricity the increases q^ are 
calculated by linear regression. Whether the rotor is displaced 
in a positive or a negative direction (cf. Figure 2 . 5)9 for a 
corresponding definition of the tangential angles the same 
pressure curves should be obtained. Hov/everi, due to minute 
differences in the seal's microgeometry 9 in part substantial 
deviations can occur. Therefore 9 in order to obtain 
representative descriptions of the pressure course 5 prf sures 
for equally large positive and negative eccentricities re 
averaged . 


4.4. Transverse forces from an eccentric rotor position and 
comparison to the efficiency curve 


/124 


Starting from the central position of the turbine step within 
the housings the eccentricity could be adjusted to up to 75 ?^ of 
the radial clearance width, in both directions. The transverse 
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forces acting on the rotor v/ere measured at 15 points o within 
this ranges the course of the forces was approssimately linear 
across the deviation (cfo[5])s for all parameter combinations 
here investigated o The pressure gradient at the turbine occurred 
as a function of the throughputs which was fixed because the 
compressor could not be varied o Correspondingly t similar test 
conditions are best differentiated for equal throughputs which 
is recorded in the Appendix^ together with other measured 
turbine step parameters o Hov/everj the necessarily somewhat 
different turbine operating conditions had only little effect on 
the consideration of the dimensionless slopes through of 
the forces across the deviationo 


4o4olo Blading v/ithout shroud band 

Figure 4ollj belowt shows the clearance excitation coefficients 
determined from force measurements for rotor blading 
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Figure 4oll Rotor without 
shroud “band 


Figxire 4.12 Rotor without 
shroud “band 
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i7ithout shroud=bando For the turbine step investigated ^ the 
clearance loss increases as the axial clearance increases 5 
thereby causing an increase in the clearance excitation 
coefficients Kgo In contrast to the linear theory of equation 
(2 o17)d they depend strongly on the radial clearance widtho The 
restoration coefficients -=Kj shovm in Figure 4ol2t above 5 
increase from an axial clearance = 1 mra on 9 which can be 
attributed to the flow line displacement illustrated in Figure 
or to differences in the pressure distribution in the 
stator clearance 0 We can not state here 9 with any certainty 9 
whether these restoration forces ■=• which according to equation 
(2ol5) should vanish = are caused by stator or rotor clearance 
losses 0 Such a statement is possible only if the housing portion 
of the rotor clearance can be displaced independently from that 
of the stators eccentrically to the rotor (cfo [30])» For a /125 
very small axial clearance there may be a finite value for the 
restoration coefficient 9 which could be measured somewhat more 
clearly during tests with standing blading [5]= For this extreme 
position of the rotor with respect to the housings due to the 
step’s construction there v/ill be a clearance loss only at the 
rotor bladings which according to equation (2ol5) should not 
give rise to restoration forces o It is thus not impossible for 
transverse forces to be caused by a non-uniform pressure 
distribution along the perimeters even for blading without 
shroud -band 0 

In order to explain the relatively large deviations of the 
measured clearance excitation coefficients in Figure 4oll from 
those expected in theorys we performed a correlation between the 
excitation forces and the measured efficiencies 9 in accordance 
with section ^o3»2o Figure ^ol3 belov/o shows the course of the 
thus referenced excitation forces for various radial clear- /126 
ance widthSs across the eccentricity. The family of curves 
plotted corresponds to the function (4. 18) 9 whose parameters 
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Figure 4ol3 Blading without shroud-band 
Curve fitted w/2nd 

order polynomial for clearance loss 
from clearance loss , [ 14 ] ,eq . (2 . 15) 
were determined from the measurement points ^ by linear 

regressions as a function of the two variables £ and £« Since 
the course of the excitation forces is linear across the 
eccentricity 5 a second order polinoraial (4»17) will be 
sufficient for the clearance losso However ^ a higher-order 
polinoraial was able to explain measurements for a standing rotor 
[ 5 ] very wellp even though there the course of the excitation 
forces is non-linear across the eccentricity j, characterized by a 
steeper slope for increasing eccentricity » 

It is possible to establish a comparison between the parameters 
of the functions {4ol8)o determined from excitation force 
measurements 0 and the measured efficiencies s for all axial 
clearances invest igatedo Starting from the calculated peripheral 
efficiency (for s = 0 )i, the broken straight line records the 
course of the clearance losses from [143 d which for this turbine 
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Figure 4»l4 Blading without shroud-band 


step agree only very little with the measurements o The slope and 
curvature of the curves plotted were determined for each axial /127 
clearance from the measured excitation forces o If together with 
equation (4,19) we were to introduce a correction factor y=^lo2p 
the agreement between these curves and the measured efficiency 
would improve 0 In other words, the excitation forces are only 
83?^ of the forces one should calculate, from the actual course 
of the efficiencies 0 Taking into consideration the certainty 
required for vibration calculations, the excitation forces in 
bladings without shroud-band can be calculated from a known 
efficiency curve, using equations ( 2 A 6 ) or 


Blading with shroud -band 

In order to reduce clearance losses, the blading is fitted with 
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© Plain shroud-band (EB 25, 28) 

* Two seal-peaks (EB 26, 27) 
o Three seal-peaks (EB 29) 
w/o shroud— band (EB 21) 

Figure 4,15 Effect of the clearance form 


shroud-bands 0 Figure 4 o15d above, shows the course of the 
measured efficiencies for equal radial clearance widths, as a 
function of the axial clearance. The efficiencies are 
significantly improved with respect to blading without shroud 
band. The recessed labyrinth with three peaks is the most 
favorable? here the seal effect of the radial clearance is 
little affected by variations in the axial clearance. Due to the 
smaller clearance losses, smaller clearance excitation forces 
are to be expected for banded blading. Figure 4,16, below, /128 
shov/s hov^ever that the clearance excitation coefficients are 
much larger than for blading without shroud -band. As shall be 
explained in section because of a pressure distribution 

that varies along the perimeter, in the rotor clearance, the 
measured forces are approximately twice as large as those 
calculated from equation (2.17) from the clearance loss only. 

With increasing axial clearance the clearance excitation /129 

coefficients become larger} there is a qualitative relationship 
to the efficiencies of Figure 4.15* Because for increasing 
clearance loss - i.e,, sma,ller efficiency - the clearance 
excitation coefficients become larger in the same manner. 

However, an exact correlation between the excitation forces and 
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• Plain shroud-band (EB 25+ 28) 

* Two seal-peaks (EB 26+27) 

■ Three seal-peaks (EB 29) 

*■" From clear. loss, eq. (2.17) 



o Plain shroud-band (EB 25+28) 
* Tiifo seal-peaks (EB 26+ 27) 

® Three seal-peaks (EB 29) 


““ w/o shroud-band (EB 2 1 ) 


Figure 4,16 Effect of the 
clearance form 


Figure 4ol7 Effect of the 
clearance form 


the efficiencies is possible onlyn when the dependence of the 
curve on the radial clearance width is loiowno Figure 4 o17p 
belowp shows the restoration coefficients across the axial 
clearances which depend strongly on the form of the rotor 
clearance 0 For a plain rotor clearances the restoring forces are 
largest and can be primarily attributed to a pressure 
distribution over the rotor shroud-band » as the comparison to 
blading without shroud-band shows o 


4o5« Pressure distribution over the rotor shroud-band and 
comparison with theory 
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4o5ol» Shroud“band v;ith plain clearance 


Figure 4ol8 shov/s a two-dimensional representation of the 



Figure 4ol8 Pressure distribution across a plain 
shroud-band 5 e = 0»7 mm 


pressure course p measured for an eccentricity of 0=7 mm at a 
plain rotor clearances both along the perimeter and in axial 
directiono The static pressure drop at the rotor blading's 
external section was approximately 60 mbarp in this measurement 
sequence (cfo Appendix) p with an afflux velocity of 140 m/so 
Starting from the entrance edge (z = O 5 cfo Figure 2o8)p the 
pressure courses for five measurement planes normal to the axisp 
in accordance v/ith Figure 4»7 (B2)p are shown along the shroud 
bands the axial direction was considerably magnified in 
comparison to the perimeters The clearance at the entrance side 
of the shroud-band here is only half as large as the average 
radial clearance widths In agreement with the definition in 
Figure 2o3s the smallest local radial clearance width is found 
at the tangential angle = Oo Because the flow is subjected to 
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torsional offsets j thers is a pressure maxiaium just before the 
narrowest clearance p which for individual measurement planes 
along the shroud=band travels almost 80 ° against the direction 
of rotation,, Due to the characteristic pressure maximum in 1 3^0 
front of the narrov/est clearance ^ a forc{; acts on the rotor that 
has a component in the direction of the ceviationo and another 
one perpendicular to ito In tangential direction the pressure 
curve does not correspond to a pure sine fimctionj since the 
maxima are steeper than the minimao 

Figure 4ol9s belov/ shov/s^ for three selected measurement 
sequences s the effect of the axial clearance on the measured 
pressure courses o with the fitted sine function plotted in each 
caseo Measurement sequences 90 and 93 are only qualitatively 
comparable to sequence 89 j already in axonometric 
representation 0 Because of a defective sealp the turbine's 
throughput and hence also the pressure and the velocity in front 
of the clearance were much smaller,, as can be seen from the 
recorded pressure PI of the central sectiono 

/132 

To calculate the theoretical pressure courses ~ also plotted on 
Figure 4ol9 - the pressure and the velocity were recalculated s 
from the turbine step's central section to the entrance radius, 
using the potential vortex law (cfo section 3<>2olo)o The loss 
coefficients introduced in section 3 'bo describe the clearance 
flow affected by torsional forces, were used 5 their effect has 
already been thoroughly investigated by means of test 
calculations!. The radial clearance was subdivided in the flow 
direction in such a manner that the pressure courses represented 
are valid for the position of the measurement planes 0 In 
addition, we plotted the course of the throughput - variable 
along the perimeter - and for small axial clearances, the course 
of the pressure at the end of the radial entranceo Due to the 
manner of construction of the seals studied, the radial exit .can 
be ignored 0 
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MEASUREMENT: For housing insert B2, fitted to sine function 



CALCULATION: Ece-"0 Ewe=Ees=&= 0^ MRe) 

ScE =0 £we=£e5'5se=0 



SQ,sO,5mm Sagsl mm So»s4mm 

fi I. O 2. 0 3. a 4. + Meas . plane 
® radial entrance (calculation) 
local throughput (calculation) 


Figure 4ol9 Pressure distribution for a plain rotor 
clearances for various axial clearances 


When the axial clearance is changed 5 it must be remembered that 
from Figure 4 o7p simultaneously the position of the measurement 
drill-holes with respect to the shroud-band is displaced = For 
the gauge ring B2 considered here^ with five holes perpendicular 
to the axisp one can obtain for s =4 mmp for instances the 

S.X 

interesting special case in which the drill-holes of the first 
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plane are already over the radial entrance » In agreement with 
theory,) v/e obtain here a nearly constant pressure course along 
the perimeter 0 By decreasing the a^cial clearances j the 
velocities at the shroud-band* s front side increase and cause 
very large pressure differences along the perimeter „ as shown by 
the calculations for s = 0o5 i’or instance o A direct 
comparison with measurements is not possible,, since the 
drill-holes for the first measurement plane are s,lready over the 
shroud-band (cf= Figure 4o7)= This measurement plane’s low 
pressure must be attributedj at the tangential angle = 180°, 
for instance, to a flow separation due to cross-section 
enlargement, which passes from a clearance width s^^^ = 0 = 5 ™ 
for the radial entrance, to the locally large radial clearance s 
= lo7 Kimo During calculations, this would be taken into account 
by means of contraction coefficients « In addition, the bend-loss 
coefficients could not be considered constant, as here, but 
dependent on the local geometry » 


/m 


Experiments with modified afflux conditions are particularly 
suited to test the calculations procedure = In a turbine, these 
conditions can be affected by the percentage reaction, which 
essentially can be modified with the rate of rotation and the 
pressure differential for the entire step= Starting from a 
central operating condition, the pressure curves for both 
possibilities of variation are shown in Figure 4 =20 5 the simple 
gauge ring B1 (see Figure 4=7) was used» It has only two 
drill-holes in the frontal and back measurement planes, which 
however are so arranged along the perimeter that the course of 
the curve can be determined from positive and negative 
eccentricity. With increasing pressure differentials the 
amplitudes of the pressure curves become larger and at lower 
rotation rates the position of the pressure maxima shifts. 


Figure 4,21, below, shows the slopes of the transverse forces 
acting on the rotor, as a function of its pressure differential. 

Due to the reduced number of measurement points, these forces /134 
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ng = 4C00 1/min ofitXIO l/min ng=6600 l/min 

Pj-PjotSmbor p, -Pj,» I7mbar . p, -p^sSImbor 

Figure 4o20 Pressure distribution for a plain rotor 

clearance,, e = 0„6 ranis s = 1 o 0 mm 

ax 

measurement: for housing Insert B1 with 
matched sine function 



pressure difference imbbbj 

Figure 4o21 Plain shroud-band (Bl) 
measurement calculation 

linear course for , o ® » 

were determined based only on the pressure curves of the central 
measurement plane o Since according to Figure 4ol8 the course of 
the pressure is non-linear in axial directions the transverse 
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foress thus dateraiinsd must to© considered only qualitativelys 
n©v©rth©l©sS(, the sara© assumptions v/ere mad© in th© coraparativ© 
ealculations shovmo Tha clearanes ©sccitation and restoration 
eo©ffioi©nts ra©asurad increas© with increasing rotor pressure 
gradients 0 However,, different tendencies are noticed for th® two 
pressure gradient variation possitoilities » The causa for this is 
not so much the change in th© loss coefficients torought atoout toy 
the rotating channel v?all (cfo Figure JoZk), tout must toe 
attritouted mainly to th© effect of the afflux velocity c^ » 

During the tests performed at constant rate of rotation n = 
8000/rainp the relative afflux energy v/as of approximately the 
same magnitude (cfo Appendix) and hencep according to Figure 
3o26 and equation (3o55)s compressive forces are a linear 

function of the pressure differentialo In contrasts during the 
variation of the rotation rate the afflux velocity was nearly 
constants due to which according to Figure 3»24 the excitation 
forces increase paratoolically v/ith the pressure differential,, 


Both tendencies were well descritoed by means of the calculations 
performed using the data in Figure 4o20o If more precise results 
can not be expected o it is only because due to a measurement 
uncertainty of Ifo in the throughput and a reading error of 1 mm 
Hg in the pressure differential for the stepp the pressure 
difference for the rotor may already be affected toy an 

error of 5 mbaro Systematic errors that could toe included in the 
simplified calculation of p^^ from equation (4ol3) have not been 
taken into consideration here. 
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Figure 4o22p toelov/g shows the clearance excitation constants 
over the axial clearance p for measurement sequences with 
approximately equal throughput o On the one handp these constants 
were determined from the transverse forces acting on the entire 
step {q_ 2 ) 0 In comparison to themp the portion q£j^ has been 
plotted which on the other hand is obtained by integration of 
the measured pressure distribution over the rotor shroud-toando 
Starting from the measurement point at s =0,5 ro^p a constant 
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feree meaBurement ©iri entire stag© 
fr©fii pressure distribution (eo'«'©r band) 
theory for rot or gap 




Figure 4»22 Plain shroud=band 
(B2) 


Figure 4»23 Plain shroud=band 
(B2) 


course across the axial clearance may be assumed for the slope 
of these forces p which is qualitatively confirmed by 
measurements with the housing insert B1 (cfo Appendixp MR 62 to 
64) 0 These slopes q^'j^ are also confirmed very clearly by 
calculationo Besides the compressive forces p the calculated 
slope q^ contains the portion q^'gii caused by the variable rotor 
clearance losso The transverse forces arising from the stator 
clearance lossp represented by the difference ^2 ~ ^2 ” *^2° 
not calculated here (however p see in Figure 3»39 for a 
qualitative comparison) o 


/1 36 

Figure 4»23 p abovep shows the restoration coefficients on a 
similar plots the portion due to the pressure distribution was 
again plotted as a qualitative course from comparable 
measurements. The restoration coefficients calculated from the 
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rotor clsarancs pressure distribution agree with the 
measurements only for large aisial elearancesc due to the 
deviation explained on the basis of Figure 4ol9= The transverse 
forces (q^^) for the entire step increase v/ith increasing axial 
clearances evidently the portion qj - q^ ■=■ q^ can not je 
neglected herso For the small axial clearance s - 0.5 mmo a 

^ EX 

transverse force acting in the direction of deviation is 
generated in the stator clearance „ which based on a sample 
calculation for a similar clearance form, had already been shown 
in Figure 3o39 (s^^^-SA)» The restoration constant q^ calculated 
only from the rotor clearance flow is alv/ays somewhat lov/er than 
the portion q^'^ from the pressure distributionj for torsional 
f loiv 0 


4o5o2o Labyrinths with tv/o seal=peaks 

The pressure curves measured in the chambers of the dioptric 
labyrinth (cfc Figure 4o7 C) for an eccentric rotor position are 
shown in Figure 4»24 for three axial clearance widths o Also 
drawn was the pressure PI in front of the rotor 5 calculated for 
the central section from measured values? it increases with 
decreasing axial clearance p because a better seal effect is then 
achieved.. Recalculated for the radius of the clearance entrance „ 
a pressure is obtained that for a large axial clearance 
approximately corresponds to that measured in front of the first 
seal peako Noticeable pressure differences are observed in this 
chamber only for a very small axial clearance » For s^^^^ = 

2 o 0 mmp the measured value deviates a little frora the fitted 
sj.ne function, which can be attributed to an untight seal at 4 » = 

0 and ip = ±l80o Due to it, locally higher velocities and hence, 
somewhat lov/er pressures, could occur » In the chamber between 
the two peaks the pressure curve’s amplitudes depend only little 
on the magnitude of the axial clearance o Behind the last peak / 138 
the pressure is nearly constant along the perimeter and 
corresponds approximately to the pressure P2 measured in the 
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MEASUREMENT: for housing insert C with fitted sine function 
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Figure 4.24 Pressure distribution for a rotor 
clearance with two seal peaks 5 e = 0,7 ™ 


central sectionj, behind the rotor 0 

For the labyrinth 0 a calculation of the pressure curves depends 
very much ^ as shown in section 3«6 ■= on the assumptions made 
for the pressure and impulse loss coefficients at the 
seal-peaks 0 Assuming a constant contraction coefficient y = 0 o 7 n 
the calculations were performed with the pressure loss 
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coefficients from equations (3o^0) and (3o^l)o Through the 
constants applicable to the seal geometry s these coefficients 
depend on the local clearance v;idths and flov; angles <, The 
impulse loss coefficient was determined as a factor y = I from 
equation (4o49), as a function of the corresponding pressure 
loss coefficients 0 Based on the considerations in relation to 
Figure 3” 29s a pressure equalizing flov/ was not taken into 
consideration 0 In additions the channel widths for the 
individual stream tubes v/ere assmned to be constants based on 
the same assumptions ioSoj that this effect is balanced by 
corresponding local variations in the loss coefficients o 

In agreement with the measurements, a relatively low pressure 
level is observed in the central chamber, caused by the first 
peaks's low loss coefficient and the changes in the flow angles 
within the clearance » The amplitudes for these pressure curves 
are approximately as high as those measured o The very small 
axial clearance s =0=5 mm caused considerable pressure 
differences in front of the first peak (cfo Figure 4ol9)p due to 
the variable velocities in the radial entrance o With the 
existing empirical loss coefficients, no better agreement 
between measurement and theory can be attained in this chamber » 

It remains unclear, in addition, v/hether the assumption is 
warranted that the pressure in the chamber before the first 
seal-peak remains constant, with turbulence and simultaneous 
redirection of the velocity. It would be conceivable that due to 
cross-currents within the chamber, the pressures determined 
through existing measurement drill-holes (cf. Figure 4.7) are 
not representative for the entire chamber. 

/ 123 . 

By means of section 3o3o^»i) the loss coefficients can be 
determined in such a way from the available measurements, that 
complete agreement exists between the measured and the 
calculated pressure curves. But this approach also depends on 
assumptions that limit the general validity of the result. 

Agreement will especially not be achieved, with the contraction 
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coefficients chosenj if we start from a relationship between 
pressure and impulse loss coefficients described by means of a 
factor Y assumed constant along the perimeter (cfo equation 
(3»^9))o Because according to Figure 3<’36p the pressure curves 
are considerably affected by even small changes in the flow 
angle o 


Assurtiing constant impulse loss coefficients 5 for instance^ the 
pressure loss coefficients shown in Figure 4o25 (A) below 


A 



-leo -»o 0 Qo 180 


Perimeter 
















*\ 









»1Bo **9 q 0 Oo IQ 0 


Perimeter 



C' Ces chamber before 1 st peak, o Cg 1 after 1 st peak 
A from measured pressure curves (Fig. 4 .24 ) , '5 = 0.2 
B used for calculation in Fig. 4. 24 with 


Figure 4o25 Loss coefficients for rotor clearance 

with two seal peaks 
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are obtained p which are confirmed qualitatively also for other 
assumptions regarding the impulse loss coefficients o 
Disregarding small deviations obtained for the limiting values ^ 
= 0 and C = Ip the measured pressure curves are described 
exactly by means of these coefficients.. However p no unequivocal 
dependence on the local geometry can be formulated p described 
for instance in terms of clearance widths p flow angles and 
distances between peakSp or also by the seal peaks' edge 
sharpness a For a large axial clearance p the loss coefficients 
C„c. the radial entrance become very small o For this reason 
the equations taken from the literature (Figure 4o25 (B)) =• as 
shovm in Figure 4o24 = yield useful results that are sufficient 
for turbine steps p since for the usual axial clearance 
construction approaches „ the axial clearance at the radial 
entrance is much larger than the radial clearance width., 


Figure 4o26 shows the measured clearance excitation constants as 
a function of the axial cleaiancep in comparison to calculations 



— Force meas. (entire stage) 

— from press .distrib . (shroud-band) 
= Theo.for rotor clearance 

Figure 4o26 Two seal peaks 



Force meas. (entire stage) 
from press .distrib. (shroud-band) 
Theo.for rotor clearance 

Figure 4«27 Two seal peaks 


154 




psrfomed v/ith the data from Figure While the transverse 

force out of the pressure distribution is nearly constant, the 
total force increases, because the clearance loss increases with 
increasing axial clearance v/idth« This course is confirmed also 
by calculations containing only the portion of clearance flov; 
coming from the rotor o The excitation force from the stator seal 
can be determined as the difference (qg - qg) between the /ikO 

measured total force and the forces determined from the rotor's 
clearance flowo While the stator clearance loss of the turbine 
steps studied is just as large as the rotor’s (cfo Appendix), 
due to the equalizing flov/s only v/eak excitation forces are 
generated here , 

/3.41 


The restoration constants are shown over the axial clearance, in 
Figure ^<.27n above. The transverse force out of the rotor 
clearance’s pressure distribution is larger, at s^^^ = 0,5 mm, 
than the total force measured at the runner, since the stator 
clearance has the same effect as in Figure 4, 23, The 
theoretically determined restoration force increases very 
steeply at small axial clearance v/idths, which may be attributed 
to the course of the pressure ahead of the first seal peak (see 
Figure 4, 24), Disregarding a small axial clearance, the 
calculation procedure provides good agreement with measurements, 
despite the uncertain loss coefficients for the transverse 
forces rotor clearance's pressure 

distributioiiD Taking into consideration the transverse forces 
generated in the stator clearance (cf, qualitatively and Qg 
in Figure 3o39)s the forces measured for the entire turbine step 
(q.-i» <lp) can also be completely explained, 

/142 

4,5»3» Off-set shroud band with three sef. . peaks 


While the dioptric labyrinth was somewhat more favorable as to 
its seal effect than the plain shroud band, it was possible to 
considerably increase the step’s effiency by means of the 
off-set labyrinth common in turbine construction (cf. Figure 
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4.15)o The pressure distributions along the periraeter o shovm in 
Figure 4o28o belov/p were measured in the chambers before the 


MEASUREMENT: for housing insert D with fitted sine function 
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before 1st peak © 1st 0 2nd labyrinth chamber 
local throughput (calculated) 


Figure 


4o28 


Pressure distribution for a rotor clearance 
with three seal peaks p e = 0»7 mm 


first seal peak and between the first and second peaks (see 
Figure k„7 (D))» From the results obtained for the dioptric 
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labyrinth p it can bs anticipated that the pressure is constantp 
beyond the last peak and corresponds to the pressure P2 measured 
beyond the rotor o 

Along the perimeter the pressure distributions are similar to 
those for the dioptric labyrinth? the measurements in the 
chamber before the first peak again lead us to suspect an 
untight sealo The pressure course for an axial clearance = 3 
mra constitute an exceptiono Here the clearance between the last 
peak and the front of the shroud-band shoulder becomes very 
smallp and the pressure in the chamber increases p due to the 
better seal effect o In this arrangement of the seal peaks with 
respect to shroud-band shoulder p the position of the pressure 
maximum is also displaced considerably against the direction of 
torsion« 

The theoretical pressure curves were calculated under 
assumptions similar to those made for the dioptric labyrinthp 
with the constants for the variable loss coefficients given 
according to the seal peak positions « Here toop because of the 
small clearance at the radial entrance p excessive pressure 
differences occur p while the agreement is better for the 
amplitudes and phase angles of the pressure curves corresponding 
to the central chambers » For the large axial clearance s =3 
mmp the pressure levels in the individual chambers can be 
reflected by corresponding loss coefficients o The pressure 
distributions along the perimeter p however p agree only little 
with the measurements p since the flow cross-sections changed due 
to the peak positions were not taken into account » In additionp 
from the position of the measured pressure maxima it must be / 144 
concluded that due to the oblique oncoming flow at the last 
peakj a redirection of the flow occurs in tangential direction » 

The slopes of the measured excitation constants are shown in 
Figure 4 o29p in comparison to the calculated values o The portion 
due to the pressure distribution is practically independent of 
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“"Force raeas. (entire stage) 

from press odistrib . (shroud-band) 
= Theo.for rotor clearance 


Figure 4«29 Three seal peaks 



Force meas. (entire stage) 

from press, distrib. (shroud-band) 

= Theo.for rotor clear. 

Figure Three seal peaks 


the axial clearance i, as was true of other clearance forms o The 
slope of the excitation force at the stator - formed as the 
difference qg - between the measured value for the entire 
step and the value calculated for the rotor = agrees with the 
results for plain rotor clearances and for the dioptric 
labyrinth 0 

The restoration constants decrease with increasing axial 
clearance,, under the effects of compressive forces t as shown in 
Figure 4o30o The restoration constants calculated frora the rotor 
pressure distributions agree only qualitatively with the 
measured values o due to deviations in the pressure curves 
already explained in connection with Figure 4o28o The difference 
q^ - q”s between the measured total force and and the / 145 

compressive forces corresponds to the compressive forces 
generated in the stator clearances if we disregard the small 
restoration forces out of the rotor clearance losso This courses 
which can also be recognized in Figure 4o27s is qualitatively 
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confirraed by the complete calculation in accordance v/ith Figure 
3 o39 tor a stator clearance v;ith radial entrance and exit.. 

It was also possible for the set=off shroud band to determine 

the loss coefficients in such a manner that the measurements are 

reproduced exactly o However, it is not possible to provide a 

generally valid dependence on the local seal geometry, which is 

here expanded with the shroud=band shoulder » It should 

nevertheless be expected for the calculated transverse forces to 

provide better agreement, for ti*e axial clearances that are 

usually much larger than the radial clearance widths o Special 

cases, though, resulting from the position of the seal peaks 

with respect to the shroud band shoulder - as occur in Figure 

koZ8 for s =3 MM “ Sire exceptedo 
ax 


4 060 Forces due to a rotor-to-houslng inclination 


Normally, according to Figure 2o5o an inclination of the rotor 
with respect to the housing occurs only coupled to rotor 
eccentricity, with the ratio of the axial clearance change a to 
the eccentricity e being a function of the deflection bending 
line, the rotor diameter and the arrangement of the turbine step 
between the bearings 0 In order to test whether the forces due to 
both kinds of deviation are additive, a separate and a common 
displacement were performed on the test turbine 0 Figure 4o31, 
below, shows a plot of the relative excitation force Q,/U over 
the eccentricity e (MR 76 ) and the force Q^/U over the 
inclination a = 0»535e (MR 80 )= Furthermore, the displacement 
was performed in such a manner that under the action of both 
effects the forces became larger in one case (MR 84), but 
smaller in the other (MR 83)0 Within the measurement precision, 
they can be readily composed from individual measurements, with 
the same conclusions also being applicable to the restoration 
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Figure kojl Relative excitation forces over the 
eccentricity and the inclination of the rotor with 
respect to the housing 
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Figure 4o32 Two seal peaks 


Figure 4c 33 Two seal peaks 
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forces and respectively,, Q^, o The ratio of the tv/o deviations 
selected here corresponds approximately to the maximum to be 
expected for a high“pov/er turbine o The excitation forces for the 
steps at the turbine entrance v/ere increased here by 
approximately lOfoi v/ith those at the exit being correspondingly 
reduced (cfo Figure 2o5)= 

Figure 4»32s above 5 shows a plot of the dimensionless slope of 
the excitation force generated due to the rotor's inclinationo 
The clearance excitation coefficient decrease^ .•tron^.Iy with 
increasing axial clearance and then attains . . tant limiting 
value? the portion originating in the pressure aistribution acts 
similarly 0 In the calculation of the theoretical coefficient 
according to equation ( 2 , 22 ), it was assumed that the entire 
clearance loss at the rotor was due to the axial seal effect » 

For this reasons for the large axial clearance the s;im of the 
theoretical forces out of the clearance loss and out of the 
pressure distribution yields excessively high values 0 Due to the 
linear clearance loss equation^ the theoretical coefficient is 
independent of the axial clearance „ The restoration coeffi- / 14? 

cients are shown in Figure 4o33s aboveo The forces originating 
in the pressure distribution depend particularly strongly on the 
axial clearance width „ hereo 

The pressure differences generated due to an inclination of the 
rotor with respect to the housing are shown in Figure 
belowp for various axial clearances o In accordance with its 
definitions the axial clearance - now variable along the 
perimeter because of the inclination - is smallest at = Oo 
Correspondingly s this is also where the highest velocities 
occurs which cause a minimum in the pressure curve in the 
chamber before the first seal pealCo For a large axial clearance 
this effect no longer exists s due to the smaller pressure 
differences. The position of the now only weakly marked pressure 
maximum in the central chamber depends on the axial clearance 
width. The pressure differences are small in comparison to those 
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MEASUREMENT? for housing insert Cp with fitted sine function 
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Figure Pressure distribution for a rotor 

clearance with two seal peaks 


generated for an eccentric ro. r position j even though very 
large axial clearance changes were tested „ such as would hardly 
be possible with a vibrating turbine shafts 

/JM 

Assuming a linear clearance excitation theory j, it is possible to 
formulate a relationship between the excitation force caused by 
the inclination, of the rotor, and efficiency curve o From 
equations (2ol9) to (2o22) we obtain the slope of the local 
clearance loss over the axial clearance, 

ilLi 7 7 kls 

sa, Ui/r “ L 


as a function of the clearance excitation coefficient K„ , 

caused by the changes in the local tangential force o It can be 

determined as the difference K_ = = K„„ between the total 

3s 3 3D 

clearance excitation coefficient and the portion due to the 
pressure distribution, from the available measurements o Figure 
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slope from escit. force meas . 


Figure 4o35 Two seal peaks 


4o35d above j shows the efficiency curve, as a function of the 
axial clearance o Starting from the calculated tangential 
efficiency for the step, the linear efficiency decrease which 
results from the clearance loss due to the shroud-hand entrance 
shoulder, is plotted, according to equation (2., 19)0 The slopes / 149 
determined from the excitation force measurements are 
qualitatively in agreement with the non-linear efficiency curve » 
However, there are deviations for large axial clearances, caused 
either by compressive forces at the stator clearance, or due to 
the fact that an inclination of the rotor can also cause small 
radial clearance changes a If we talie these side effects into 
consideration, the relationship between the excitation forces 
and the efficiency course seems sufficiently confirmed also for 
the case of an inclination between rotor and housingo 
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5 SUMMRY 


The Totovz of thermal turboraachines can become endangered by 
self-=excited vibrations that may impose severe output 
restrictions on the affected power plants 5 if these vibrations 
are caused by clearance ejccitationo For a theoretical 
description of the systeraj it is necessary to know the 
transverse forces that act in dependance on a deviation of the 
rotor from its centered position v/ith respect to the housingo 
The component perpendicular to the direction of deviation 
acquires special significance 5 since in circumpolar vibrational 
movement it act 5 by setting up vibrations in the turbine shaft » 

The excitation force can be described as the resultant of the 
turbine stage's tangential forcej which becomes variable for an 
eccentric rotor positiono Therefore,, for a given rotation rate„ 
it is a function of the stage’s output and of the course of the 
clearance loss over the radial clearance width c In the case of 
blading with shroud-band 9 the excitation forces can be 
considerably enhanced due to a variable pressure distribution in 
the rotor’s seal clearance » This effect was studied in detail by 
means of a calculation procedure in which the torsionally 
affected flow at the clearance - caused by the main flowp in 
turbine stages - was especially taken into account 0 

Sample calculations show that the excitation forces due to a 
pressure distribution in the seal clearance increase with 
increasing pressure gradient and especially, with increases in 
the tangential velocity before the clearance,, However, the 
characteristic quantity is the relative afflux energy, to be 
obtained from the dynamic pressure of the tangential velocity 
before the clearance, related to the pressure gradient operating 
at the seal clearance. Correspondingly, because of the large 
relative afflux energy at turbine stages in impulse 



construction!, larger excitation forces become active due to the 
pressure distribution over the rotor’s shroud=band 5 than in /15I 
reaction stages = However!, the form of the seal clearance and the 
empirical loss coefficients!, which describe the clearance flov/ 
as a function of the geometry!, also have a significant effect on 
the magnitude of the excitation forces. 

The experimental studies were performed on an impulse turbine 
stages varying the rotor’s clearance form. The rotor was mounted 
in a kind of two-component balance,, connected to a 
dynamometrical brakes the transverse forces were measured as a 
function of its deviation with respect to the housing. By means 
of efficiency measurements at various radial clearance widths,, 
it was possible to show that for blading without shroud-band the 
excitation forces could be calculated from the clearance loss 
alone. In contrasts for the bladings with shroud-bard 
investigated 0 in spite of lov/er clearance losses 5 larger 
excitation forces were observed 5 these could be explained by 
means of measurements of the pressure distribution at the rotor 
clearance . 

Considering the complexity of the flow processes in a seal 
clearances the agreement between the measured pressure curves 
and the theory is relatively goods minor deviations - especially 
for small axial clearance widths » can be explained in terms of 
the simplifying assumptions made regarding the loss 
coefficients. For the usual constructions 5 the transverse forces 
acting on a turbine stage - caused by the clearance losss on the 
one hands and the pressure distribution along the rotor 
perimeters on the other - can be determined with sufficient 
accuracy. 
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6 DESIGNATIONS USED 
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surfaces 8 flow cross=sections 

specific works axial clearance change due to inclination 

seal clearance length 

relative afflux energy 

absolute velocity at blading 

specific heat capacity 

diameter 

insert condition to describe stage geometry 
eccentricity of the rotor with respect to the housing 
control magnitude (plain clearance = labyrinth) 
chamber height 
isentropic gradient 
clearance loss coefficients 

coefficient for the transverse forces of the flow medium 
blade lengthy flow paths in the control space 
measurement sequence number 

massj constant for seal peak loss coefficients 

throughput 

rotation rate 

outputs compressive force on control space 
pressure 

transverse force of the flow medium 

coefficient for the transverse force depending on the 

deviation 

Reynold’s number 

radii 

support force 
clearance width 
temperature 

distance of the seal peaks 

tangential or peripheral forces isentropic tangential 
force 



u tangential or peripheral velocity 

Vp w velocities within the clearance 

w relative velocity at the blading 

XpypZ spatial coordinates 

z' p z" number of seal peaks 
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ot flow angles angle of rotor inclination 

oip 3 flovi angles at turbine blading 
Cp C pressiire loss and impulse loss coefficients 
n efficiency 

ae isentropy exponent 

A friction coefficient 

y contraction coefficient 

P density 

tangential or peripheral angle p bending angle 
<1) pressure coefficient of the turbine stage 

0 ) angular velocity 

Jjp « load vector 

load matrix proportional to the deviation 
load matrix proportional to the velocity 
percentage reaction 
Subscript and superscripts s 

^ dimensionless representation 

clearance width variable along perimeter 


t 


’ stator 

” rotor 

1 in the direction of deviation 

2 perpendicular to direction of deviation (preceding 
the direction of rotation) 

A radial exit 

B reference magnitude 

D from pressure distribution 

E radial entrance 

ipn support point in flow direction 

kpj support point in tangential direction 
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S radial clearance 

s originating in clearance loss 

Stage control surfaces? 

0 before the stator 

1 between bladings 

2 behind rotor 
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8 APPENDIX? TABLE OP mSURED VALUES 


The Table belov/ lists the mid-section magnitudes of the test 
turbine and the slopes of the forces over the eccentricity for 
selected measurement sequences® The stage geometry is described 
by means of the insert as defined in section ^o2o The individual 
magnitudes are explained in the sequence of the computer 
print-out 5 with a listing for the equation number® 


MR Measurement sequence number 
SA]C axial clearance 
NB operating speed 
M throughput „ (4 ® 1 ) 

DHS isotropic gradient p (4«2) 

PSI pressure coefficient p (j) = 2Ah^/u 

US isentropic tangential force 

ETAE effective efficiencyp (4o4) 

ETAU tangential efficiencyp (4o5) 

REAK percentage react ionp (4ol4) 

REMA percentage reaction from measurement p^^ (cf® Figure 4®4) 

KSI ratio corresponds to q2s'^^2s 

Q2S theoretical excitation force „ (2® 15) 

US/L reference magnitude U /l for following constants 
Q2 excitation coefficient (entire stage) 

Q1 restoration coefficient (entire stage) 


Q2D excitation coefficient from pressure distribution 

QID restoration coefficient from pressure distribution 

CES relative afflux energy for rotor clearancep (3«24)p 
recalculated for external tadius using (3“8) 

PO pressure before turbine stage p^ 

DP02 pressure difference p^ - p^ 

DP12 pressure difference p^ - P 2 from (4®13) 

TO temperature before stage [°C] 

CO afflux velocityp (4®6) 

Cl stator exit velocity, (4® 15) 

C2 velocity behind rotor, (4.15) 

AL2 outflow angle, (4.15) 
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MR SAX N0 M QMS PSJ US ETAE feTAJ ETAU KEAK REAM KS2 QE§ 

HH yPH K6/S RJ/K6 N »»—••« N/MH 

Insert cond. 231 w/o shroud band (s=0.5 mm) 'Form A 

19 0*5 8000 *399 iA,9 5.90 76,9 ,720 ,70A ,768 ,222 ,068 ,57 7o2 

60 1,0 BO0O ,399 i6,7 9,60 78,2 ,72X ,69* ,790 ,2;3 ,069 ,97 7,2 

61 2,0 BOOO ,389 16, e 9,91 70,7 ,7X0 ,7X9 ,769 ,219 ,068 ,97 7,2 

62 3,0 80eo ,399 16,7 g.68 78,2 ,7X6 ,720 ,79u ,212 ,068 ,97 7.2 

Insert cond. 211 x^i/o shroud-band (s=1.0 mm), Form A 

X7 0,9 8000 ,393 X6,6 5,37 77,3 ,628 ,666 ,766 ,165 » ,57 7.X 

X6 1,0 8000 ,393 X6,X 9,26 79,8 ,996 ,65X ,T6V ,i2l n ,57 7.0 

X6 2,0 8000 ,392 X6,0 5,20 76,7 ,979 ,666 ,77X ,117 o ,57 6.9 

19 3,© 8000 ,39?. 16,1 5,25 79,6 ,578 ,639 ,769 ,12X " ,S? 7.0 

Insert cond. 221 w,/o shroud-band ,(s=l .5 mm). Form A 


io 

l.O 

8000 ,398 

16,0 5.22 

76.3 

,921 

,357 ,772 

,1U5 

,071 

,5? 

7,0 

n z,o 

BOOO ,398 

16,0 9,22 

76.3 

,459 

,397 ,772 

,105 

.071 

.57 

7,0 

27 


aoOQ ,399 

16.0 5,21 

76il 

,456 

,331 .772 

,105 

,072 

,57 

7,0 

32 

3f 0 

8000 ,398 

16,0 9,22 

76.3 

,440 

,357 ,772 

,105 

,071 

,57 

7,0 

Insert 

cond. 251 

1 Plain 

shrouu-band, Form 

B1 





1,0 

6000 ,366 

12,316,30122,8 

,405 

,466 ,96-. 

,091 

n 

,38 

7,9 

00 

1.0 

6000 ,368 

12,6 8,36 

64.9 

,617 

,617 ,709 

,117 

n 

,47 

6,B 

66 

l.@ 

8000 ,369 

13,6 9.06 

71.9 

,67? 

,923 ,769 

• 145 

n 

,49 

5o5 


1.0 

eoug ,622 

13.9 9,93 

91,7 

,659 

,560 ,753 

,139 

n 

,50 

7,3 


0,5 

BUOO ,629 

16,7 6,22 

97,7 

,707 

,699 ,742 

,169 

Of 

,49 

7,6 

U 

1,0 

GOOQ ,629 

16,6 6,12 

96,1 

• 666 

,939 ,747 

,157 

n 

,50 

7,7 


2.0 

S§0@ ,629 

16,2 6.02 

94,9 

,659 

,922 ,793 

,132 

n 

.50 

7.4 

Insert 

cond. 281 

Plain : 

shroud- iband, Form B2 




99 

0,9 

8000 ,623 

16,9 6,29 

97,7 

,707 

,770 ,741 

.171 

.041 

,49 

7o6 

— 

•Throu^^hput corrected 

MR ! 

30 to 93— 



90 

l.O 

600Q o326 

10,0 3.72 

44,6 

,671 

,748 ,788 

,098 

,046 

,50 

3,5 

91 

^ e 0 

6Q0Q 

lOoO 

44,1 

,660 

,722 ,786 

,131 

,0S4 

,50 

3,5 

92 

3,0 

800© 0^26 

10,1 

45,4 

,650 

,744 ,708 

,100 

«0f? 

,50 

3,6 
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10 ,3 3 g o ^ 

46 , 5 

,649 

,740 ,709 

,091 

,059 

,49 

3 o6 

Insert 

cond. 271 

Two seal-peaks 3 

Form C 





06 

Oof 

6000 9^00 

16,2 6,03 

80,4 

• 6Q8 

,731 .725 

,286 

U7S 

,43 

6o4 

7S 

uo 

8000 

16, e 6,26 

95,3 

,664 

,727 ,734 

,217 

aofl 

,43 

6 e 6 

76 

a«0 

6000 0^17 

16,9 6,26 

95,4 

,652 

,706 ,735 

,210 

ol^O 

,42 

6,6 

77 

3g0 

60UQ 

16,7 6,21 

95,1 

,642 

,701 .737 

g2Ul 


,42 

6*5 

78 

^ 0 0 

ByO0 

16,7 6,20 

95,0 

,635 

,704 ,737 

b200 

bW 

,42 

6o 5 

Insert 

condo 291 

Three seal- 

■peaks 3 Form D 




113 

0,5 

6000 8^06 

i<»,8 

82.3 

,700 

,774 ,756 

.175 

eg 

,37 

5,^ 


1,0 

8000 e^O? 

1^,8 5,52 

82,4 

, 608 

,762 ,756 

si'' 

81 

,37 

5g2 

US 

2,0 

6000 ,<>07 

14,9 5eS4 

62,6 

®6Q7 

,764 ,756 

,1’fS 

«3 

,37 

5©a 

lit 

3,0 

8000 

14,9 5,56 

82.8 

,696 

,756 ,754 

.18' 

(n 

,37 

5,2 
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MR US a 

02 

°Qi 

02D « 


ees po 

0PO2 OPU T0 

eo 

ei 

C2 

4^2 


N/HM 


N/MM 

N/HW 1 

fVHM 

OOC3 HBAR 

MOAR MBAR g 

M/s 

M/s 

M/S 

OR0 

Insert ^ 

cond . 

231 

x-j/o 

shroud=“band 

(s“0 .5 Iran) 

5 Form 

A 


39 

4.1 

6.1 


m 


1.4 U43 

167 39.3 23 

?4 

145 

41 

64 

^0 

4.S, 

o.a 

4.^ 

ts 

» 

1.4 U42 

109 37,4 24 

a4 

145 

41 

65 

H 

4.1 

10.7 

4|l 

m 

ft 

1.4 U43 

186 38, U 24 

a4 

145 

41 

64 

42 

4.1 

{>2.0 

Si6 

(0 

fl 

U4 1142 

109 37,3 26 

^4 

145 

41 

65 

Insert 

cond. 

211 

w/c 

1 shroud“band 

(s= 1 .0 Iran) 

5 Form 

A 


V 

4.1 

4.6 

Oi? 


V 

l.B 1129 

179 24,9 29 


15© 

37 

71 

It 

4«0 

6.1 

Qt^ 

• 

f9 

2,0 1125 

179 20,0 25 

25 

151 

35 

74 

14 

3.9 

e.e 

Xii 

m 

n 

2.1 1123 

173 19,2 24 

25 

150 

34 

75 


4.0 

. 8.4 

ii'f 


n 

2,0 U25 

179 19,9 24 

25 

15© 

33 

74 

Insert 

cond. 

221 

Vj/c 

1 shroud-band 

3 (S“ 1 .5 mm) 5 

Form A 


do 

4o0 

5.0 

“08 2 


n 

3,2 XX33 

175 17,3 25 


151 

34 

76 

31 

4.0 

0.5 

08^ 

e> 


2,S li3. 

179 17,4 25 

25 

151 

34 

76 

27 

4.0 

0.3 

0 8 5 

«3 


2,1 U34 

179 17,4 23 

25 

151 

34 

76 

32 

4.0 

7.1 

le^ 

m 

* 

2,2 113. 

175 17,3 25 

(^9 

151 

34 

76 

Insert < 

cond . 

251 

Plain shroud-band 5 Form B1 





59 

6.5 

9.5 

3. a 

4.5 

2,4 

2,6 1092 

149 19.0 £4 

24 

143 

72 

33 

6U 

4.4 

10.1 

5,8 

4.6 

3.x 

2,1 1100 

154 17.2 24 

^4 

142 

47 

40 

6a 

3.6 

JO. 7 

7,5 

5{)Q 

3,6 

l,e 1117. 

170 g3.3 21 

24 

146 

34 

75 

63 

4o0. 

J2.4 

8 j6 

S.9 

4, a 

1,9 1150 

292 Sa,4 21 

26 

158 

41 

64 

64 

5.1 

10,3 

6,8 

5.1 

3,3 

1,7 1161 

214 33.V 21 

?6 

159 

45 

59 

62 

3.1 

J4.0 

8i^ 

6.2 

4,6 

l.a 115S 

211 3J.1 19 

$6 

159 

44 

61 

63 

3.0 

|0.2 

xoii 

6.2 

4,1 

1,9 1154 

207 25,6 $0 

26 

160 

42 

63 

Insert 

cond. 

281 

Plain 

shroud-band 5 Form B2 





t)9 

5.1 

9.9 

7,5 

7,5 

9.x 

1.7 U60 

213 34.2 25 

26 

160 

46 

5G 

— 

•Throughput corrected 

MR 90 to 

93 





90 

2.^ 

7.5 

6.2 

4,6 

6.1 

2,2 1067 

119 U.3 25 

22 

126 

27 

U4 

n 

2.3 

0.7 

7.8 

3.9 

3,0 

2.0 1067 

120 15,0 26 

aa 

126 

a? 

109 

92 

2.4 

9.7 

8.'^ 

3.9 

5.4 

2.2 1069 

121 U.^ 26 

22 

129 

27 

iia 

93 

2.4 

9.5 

9.2 

3,6 

4,5 

2,3 1070 

123 10,7 26 

22 

lil 

27 

iia 

Insert - 

cond. 

271 

Txiro 

seal-peaks 

Form C 





60 

4.7 

7gC 

7.2 

5,6 

9,7 

1.1 U75 

20B 56.2 26 

24 

145 

49 

56 

75 

5.0 

JOoC 

6.6 

6.3 

e,o 

1.4 U67 

213 43,6 25 

26 

155 

40 

57 

76 

5.0 

ja.3 

6. r 

5,e 

4,4 

1.4 U66 

213 42.1 25 

26 

155 

40 

57 

77 

5.0 

ja.o 

t % f 

3.4 

3,0 

1.5 U65 

212 40.0 25 

26 

156 

47 

58 

70 

5.0 

13,2 

8.3 

5,3 

4,X 

lo5 U65 

211 39, § 25 

26 

156 

47 

58 

Insert cond. 

291 

Three seal-peaks 

, Form D 





U3 

4.3 

6,7 

5.2 

4,0 

6.2 

1.6 1164 

191 31.5 23 


149 

39 

66 


4.3 

0,i 

4.9 

3.7 

4.4 

1.7 U64 

191 31,3 23 

25 

150 

39 

66 

U3 

4.3 

9,<i 

3eQ 

3.5 

0,3 

1.6 1164 

191 31,6 23 

25 

150 

39 

66 

lit 

4.4 

B,6 

U4 

3.4 

3.1 

1,6 1169 

192 32.7 23 

(59 

149 

40 

66 
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